For  Reference 


NOT  to  be  taken  from  this  room 


©X  llBUS 
1MMI 


The  University  of  Alberta 
Printing  Department 
Edmonton,  Alberta 


UNIVERSITY  OF  ALBERTA 
LIBRARY 


Regulations  Regarding  Theses  and  Dissertations 


Typescript  copies  of  theses  and  dissertations  for  Master’s  and  Doctor’s 
degrees  deposited  in  the  University  of  Alberta  Library,  as  the  official  Copy  of 
the  Faculty  of  Graduate  Studies,  may  be  consulted  in  the  Reference  Reading  Room 
only. 


A  second  copy  is  on  deposit  in  the  Department  under  whose  supervision  the 
work  was  done.  Some  Departments  are  willing  to  loan  their  copy  to  libraries, 
through  the  inter-library  loan  service  of  the  University  of  Alberta  Library. 

These  theses  and  dissertations  are  to  be  used  only  with  due  regard  to  the 
rights  of  the  author.  Written  permission  of  the  author  and  of  the  Department 
must  be  obtained  through  the  University  of  Alberta  Library  when  extended  passages 
are  copied.  When  permission  has  been  granted,  acknowledgement  must  appear  in  the 
published  work. 

This  thesis  or  dissertation  has  been  used  in  accordance  with  the  above 
regulations  by  the  persons  listed  below.  The  borrowing  library  is  obligated  to 
secure  the  signature  of  each  user. 


Please  sign  below: 


THE  UNIVERSITY  OF  ALBERTA 


COMPUTER  CALCULATION  OF  MULTICOMPONENT 
VAPOR-LIQUID  EQUILIBRIA  AND  ITS 
APPLICATION  TO  THE  DESIGN  OF 
DISTILLATION  COLUMNS 

BY 

ANTOINE  M.  HAKIM 


A  THESIS 

SUBMITTED  TO  THE  FACULTY  OF  GRADUATE  STUDIES 
IN  PARTIAL  FULFILLMENT  OF  THE  REQUIREMENTS  FOR  THE 
DEGREE  OF  MASTER  OF  SCIENCE 


DEPARTMENT  OF  CHEMICAL  AND  PETROLEUM  ENGINEERING 


EDMONTON ,  ALBERTA 
July,  1967 


ATflMJA  **0  YTI25I3VIWU  3HT 


8TI  <3 HA  AIHSIJIUQ3  CIIUQIJ--H0SAV 
*50  KOI  33a  3HT  OT  K0ITA0IJ33A 
3KMLM0D  W0ITAJJIT3IQ 


MIKAH  .M  3KI0TKA 


3I33HT  A 

23IGUT8  3TAUCIA30  30  YTJU0A3  3HT  OT  G3TT I M8U2 


3HT  303  3TK3M3HIUQ33  3HT  30  TK3MJJI3J  3  :  ^AIP>IA3  41 


3DK3I03  30  33T3AM  30  33H03C 


ATH3SJA  , H0TK0MG3 


Tdei 


UNIVERSITY  OF  ALBERTA 


FACULTY  OF  GRADUATE  STUDIES 


The  undersigned  certify  that  they  have  read,  and 


recommend  to  the  Faculty  of  Graduate  Studies  for  acceptanc 
a  thesis  entitled  'Computer  Calculation  of  Multicomponent 
Vapor-Liquid  Equilibria  and  its  Application  to  the  Design 
of  Distillation  Columns"  submitted  by  Antoine  M.  Hakim  in 
partial  fulfillment  of  the  requirements  for  the  degree  of 
Master  of  Science  in  Chemical  Engineering. 


Date  . 


npisaQ  o3  cioiJ&oi IqqA  bne  & iidl  1  iup3  bxLfpi.I  •  icqs^ 


ACKNOWLEDGEMENTS 


The  author  expresses  his  appreciation  to  the 
following  persons  and  organizations. 

Dr.  F.D.  Otto  for  his  guidance  and  advice; 

Syncrude  Canada  Ltd.  for  its  encouraging  and 
understanding  attitude; 

Chemcell  (1963)  Ltd.  for  providing  the  necessary 
information  for  the  study  of  the  acetone-methanol  tower. 

The  staff  of  the  University  of  Alberta  Computing 
Center  for  the  computer  time  required  during  this  work. 


ssbyitx^s  pnibne^s^bni/ 


■ 


ABSTRACT 


The  semi-empirical  model  for  the  calculation  of 
vapor-liquid  equilibrium  presented  by  Prausnitz  is  evalu¬ 
ated.  The  Prausnitz  model  uses  the  Wilson  equation  as  the 
source  for  the  multicomponent  liquid  activity  coefficients. 

The  properties  of  the  Wilson  parameters  and  the  interaction 
energies  are  investigated,  and  two  optimization  techniques  are 
used  to  locate  the  Wilson  or  interaction  energy  parameters 
that  will  best  predict  the  vapor-liquid  equilibria  in  a  system. 

The  Prausnitz  equilibrium  model  is  incorporated 
in  a  distillation  routine  which  uses  a  tridiagonal  matrix 
method  to  solve  the  material  balance  equations  and  the  Muller 
method  to  calculate  bubble-point  temperatures.  This  distil¬ 
lation  program  is  used  to  study  and  optimize  the  extractive 
distillation  of  the  system  MCH-toluene  with  phenol.  The 
same  program  is  used  to  simulate  the  Chemcell  Ltd.  distilla¬ 
tion  of  acetone-methanol  with  water  as  the  extractive  agent. 

The  simulation  succeeds  in  duplicating  the  product  purities. 

The  program  is  then  used  to  determine  more  optimal  operating 
conditions  for  this  industrial  column. 
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I .  INTRODUCTION 

The  sound  design  of  equipment  for  separation  opera¬ 
tions  necessitates  the  availability  of  reliable  vapor-liquid 
equilibria . 

Experimental  data  for  binary  vapor-liquid  equilibria 
are  often  available,  but  experimental  ternary  equilibria  are 
rare,  and  for  systems  with  more  than  three  components,  labora¬ 
tory  data  are  practically  non-existent.  Thus,  the  value  of  a 
semi-theoretical  model  that  could  predict  multicomponent  vapor- 
liquid  equilibria  using  information  from  binary  equilibria 
only  is  obvious. 

Several  methods  for  the  prediction  of  multicomponent 
equilibria  are  available.  Carlson  and  Colburn  and  Wohl 
have  presented  excellent  reviews  of  these  methods.  All  of 
these,  however,  involve  long  and  tedious  calculations  in 
their  rigorous  forms  and  the  short  cuts  that  have  been  used 
to  speed  their  solution  have  led  to  a  lack  of  accuracy. 

The  development  of  large-memory  computers  has  en¬ 
couraged  the  formulation  and  use  of  several  semi-theoretical 
vapor-liquid  equilibrium  models  which  are  applicable  to  wide 
ranges  of  components.  Two  such  models  are  the  Chao-Seader 
and  the  Prausnitz  models,  the  latter  based  on  the  equation 
of  Wilson  .  The  Chao-Seader  model  is  primarily  designed  to 
handle  hydrocarbon  mixtures,  but  it  can  also  predict  the 
equilibria  for  several  non-condensible  gases  and  boiling 
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point  cuts  of  petroleum.  The  Prausnitz  model  predicts  multi- 
component  equilibria  using  liquid  activity  coefficients  ob¬ 
tained  from  Wilson's  equation. 

Orye  and  Prausnitz  were  the  first  to  conclude  that 
the  thermodynamic  model  presented  by  Wilson  did  actually 
provide  a  sound  model  for  prediction  of  multicomponent  equil¬ 
ibria.  Since  this  model  appeared  to  handle  polar  components 
well,  it  was  decided  to  study  the  possibilities  it  offered  in 
the  prediction  of  multicomponent  vapor-liquid  equilibria. 

One  objective  of  this  work  then  is  to  achieve  a 
better  understanding  of  the  Prausnitz  model  and  test  its 
ability  to  predict  multicomponent  equilibria  for  a  wide  range 
of  systems.  Another  objective  is  to  develop  a  computer  rou¬ 
tine  which  can  be  used  for  the  reliable  design  of  equipment 
for  multicomponent  distillation.  The  incorporation  of  the 
Prausnitz  and  the  Chao-Seader  models  in  such  a  routine 
would  have  special  applications  in  the  design  of  extractive 
distillation  equipment. 
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II  -  LITERATURE  SURVEY 

2.1  Determination  of  Vaporization  Constants 

A  comprehensive  review  of  the  methods  available 

for  multicomponent  equilibria  prediction  prior  to  1946  is 

1  2 

presented  by  Carlson  and  Colburn  and  Wohl  .  Since  the 

development  of  high-speed  digital  computers,  several  models 

and  data  correlations  designed  for  computer  calculations 

are  now  available  t^e  most  notable  being  the  Chao- 

3  .  4 

Seader  and  Prausnitz  correlations. 

3 

The  Chao-Seader  correlation  provides  equilibrium 
ratios  for  hydrocarbons  of  various  molecular  structures, 

7 

some  common  non-condensibles  (CC^,  I^S,  ^ CO,  O^)  , 
and  heavier  hydrocarbon  fractions,  all  over  wide  pressure 
and  temperature  ranges. 

4 

A  recent  publication  by  Orye  and  Prausnitz 

proposes  a  semi-empirical  model  for  the  calculation  of 

vapor-liquid  equilibria.  The  model  is  studied  by  Lyckman 
8  9 

et  al  ,  O'Connell  et  al  ,  and  more  comprehensively  in  a 
recently  published  nomograph"^. 

The  Prausnitz  model  calculates  independently  the 
fugacity  of  the  pure  liquid  and  the  fugacity  coefficient  of 
the  vapor  phase.  A  calculation  of  the  latter  at  low  system 

pressures  is  possible"^  through  the  use  of  the  Pitzer  and 

12  13 

curl  model  presented  by  Lewis  and  Randall 
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The  liquid  phase  activity  coefficients  are 

14 

calculated  using  the  Wilson  model  .  This  model  uses  the 
15  16 

Flory  and  Huggins  concept  of  athermal  solutions  whereby 
excess  enthalpy  is  set  to  zero. 

2.2  Determination  of  Enthalpy  Data 

Several  correlations  have  been  developed  to  provide 

3  6  17  18 

enthalpy  data  using  computer  calculations  '  '  '  .  The  Chao- 

3 

Seader  correlation  provides  enthalpy  data  that  are  thermo¬ 
dynamically  consistent  with  the  equilibrium  data  calculated 
by  the  correlation.  This  is  due  to  the  fact  that  the  departures 
from  ideal  enthalpies  are  determined  using  parameters  obtained 
in  the  calculation  of  the  vaporization  constants. 

The  Chao-Seader  correlation  does  not  cover  many  of 
the  systems  commonly  separated  by  extractive  or  azeotropic 
distillation.  In  these  systems  vapor-phase  ideality  is 
assumed,  but  partial  molar  heats  of  solution  have  to  be  cal¬ 
culated  for  the  liquid  phase.  These  partial  molar  heats  of 

19 

solution  are  calculated  by  the  Redlich-Kister  equation  know¬ 
ing  the  liquid  phase  composition  and  the  infinite  heats  of 

20 

dilution.  The  latter  are  listed  in  the  Critical  Tables 

for  many  common  systems  or  can  be  obtained  by  a  combination 

21 

of  the  Pierotti  correlation  and  the  relationship  of  the 

19 

activity  coefficients  to  the  infinite  heats  of  dilution 
Ideal  enthalpies  can  be  obtained  on  most  hydrocarbon  systems 
from  the  literature^  '  ^  . 
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2.3  Rigorous  Distillation  Calculations 

The  development  of  high-speed  digital  computers 

has  made  the  rigorous  solution  of  a  distillation  column 

possible.  The  most  common  computational  techniques  for 

the  solution  and  design  of  distillation  columns  are  varia- 

23 

tions  of  the  methods  of  Thiele  and  Geddes  and  Amundson 

Both  of  these  methods  rely  on  the  iterative  technique  and 

24  19 

are  extensively  covered  by  Holland  and  Smith  respectively. 

25 

Lyster  has  made  an  extensive  study  of  the 

method  of  Thiele  and  Geddes  as  applied  to  simple  and  complex 

23 

distillation  columns.  The  model  developed  by  Amundson  was 

the  first  attempt  at  setting  the  column  material  and  enthalpy 

23 

balance  equations  in  a  matrix  form.  Amundson  solved  the 

matrix  by  inversion,  however,  which  introduced  instabilities. 

The  simultaneous  solution  of  the  column  matrix 

using  a  modified  Gauss  elimination  technique  was  lately  pro- 

2  6 

posed  by  Wang  and  Henke  .  This  method  of  solution  is 
2  6 

claimed  to  be  numerically  more  stable  and  more  accurate 

than  the  inversion  method  since  it  totally  avoids  the  sub- 

27 

traction  of  nearly  equal  quantities.  The  Muller  method  ' 

which  is  claimed  to  cut  solution  times  by  at  least  a  factor 

2  8 

of  three  over  the  Erbar  technique  ,  was  used  in  this  work 

to  determine  the  bubble-point  stage  temperatures.  The 

29 

method  is  studied  and  recommended  by  Lapidus 
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2.4  Rigorous  Design  of  Extractive  Distillation 

Relatively  few  studies  have  been  reported  on 

the  rigorous  design  of  extractive  distillation.  The  topic 

19 

is  dealt  with  by  Smith  who  presents  the  first  iteration 
in  the  design  of  an  extractive  distillation  column  for 

30 

the  reparation  of  methylcyclohexane  and  toluene.  Coates 

describes  the  phenomena  that  take  place  inside  an  extractive 

distillation  column  and  lists  the  general  properties  that 

a  solvent  must  possess.  The  choice  of  the  solvent  can  be 

31 

made  using  the  molecular  similarity  technique  where  the 

solvent  is  chosen  according  to  the  hydrogen-bond  capability 

that  both  the  system  and  the  solvent  possess,  or  by  use 

21  32 

of  the  Pierotti  correlation.  Dunn  et  al  provide  classical 
data  as  to  the  effect  of  solvent  concentration  on  the  vapor- 
liquid  equilibria  of  a  system  to  be  separated.  They  also 
describe  the  result  of  their  work  on  the  optimization  of 
an  extractive  distillation  column.  As  for  the  simulation 
of  multicomponent  distillation  columns,  a  research  in 
the  literature  could  not  locate  any  work  dealing  with 
that  topic. 
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III.  FUNDAMENTAL  CONSIDERATIONS  IN 
EXTRACTIVE  DISTILLATION 

3 . 1  General 

Whenever  the  components  in  a  system  to  be  fraction¬ 
ated  boil  too  closely  together  to  be  economically  separated 
by  simple  distillation,  azeotropic  or  extractive  distillation 
may  be  resorted  to.  In  both  these  cases  a  third  component, 
called  '!entrainerr  or  "solvent'  respectively  is  added  to  the 

system.  This  addition  causes  a  modification  of  the  deviations 

19 

from  Raoult's  law  that  each  component  exhibits 
This  modification  is  obtained  by: 

a)  The  choice  of  the  third  component 

b)  The  concentration  of  the  third  component  in  the  system. 

In  the  case  of  azeotropic  distillation  the  "entrain- 
er"  forms  an  azeotrope  with  one  or  both  of  the  components  in 
the  system.  This  facilitates  the  separation  since  most 
binary  azeotropes  are  lower  boiling  than  any  of  the  pure 
components  that  constitute  the  azeotrope.  In  extractive  dis¬ 
tillation,  the  "solvent"  ideally  enhances  the  volatility  of 
the  naturally  more  volatile  component  in  the  system.  The 
solvent  itself  boils  so  far  above  any  of  the  pure-feed  compo¬ 
nents  that  the  possibility  of  forming  azeotropes  is  usually 
quite  remote.  This  also  causes  it  to  always  leave  a  distil¬ 
lation  column  at  the  bottom. 
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3 . 2  Desirable  Solvent  Properties 

30 

Coates  lists  the  general  properties  that  a  solvent 
must  possess  for  the  successful  extractive  distillation  of  a 
system.  These  are: 

1)  The  solvent  should  be  non-corrosive,  readily  available, 
non-reactive  with  the  components  to  be  separated,  and  of 
much  lower  volatility  with  respect  to  them. 

2)  It  should  cause  a  greater  increase  in  the  volatility  of 
one  component  than  in  that  of  the  other. 

3)  The  solvent  should  avoid  the  formation  of  immiscible 
liquid  phases  and/or  azeotropes  with  either  of  the  com¬ 
ponents  to  be  separated. 

3 . 3  Selection  of  the  Solvent 

One  general  approach  in  the  selection  of  a  solvent 
is  to  find  a  compound  that  is  chemically  similar  to  the  high 
boiling  key.  A  high-boiling  homolog  to  this  compound  is  then 
used  as  the  separating  agent  to  avoid  azeotropic  formation 
between  the  chemically-similar  compounds. 

The  fact  that  the  solvent  is  similar  to  the  higher 
boiling  key  causes  it  to  behave  more  ideally  with  that  key 
than  with  the  light  key.  This  causes  the  light  key  to  exhibit 
greater  positive  deviations  from  Raoult's  law. 
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Table  3.1  illustrates  this  method  for  acetone  and 

o  — 1  Q 

methanol  with  respective  boiling  points  of  56.4  and  64.7  Cr  . 
If  acetone  is  the  heavy  key  in  a  mixture  to  be  separated, 
then  any  of  its  four  listed  homologs  could  be  used  as  a  sol¬ 
vent  provided  it  does  not  form  an  azeotrope  with  the  light- 
key.  If  acetone  is  the  light  key,  such  as  in  an  acetone- 
methanol  mixture,  water  would  probably  be  the  most  suit¬ 

able  solvent  due  to  its  ready  availability  and  its  easy 
handling  characteristics.  According  to  this  method,  then, 
the  lowest-boiling  homolog  of  the  heavy  key  that  does  not  form 
an  azeotrope  with  the  light  key  is  the  most  suitable  solvent. 


Table  3.1 

Homologs  of  Acetone  and  Methanol  Respectively 


Solvent 

B.P. 

°C 

Solvent 

B.P. 

°C 

Methyl 

Ethyl  Ketone 

79.6 

Ethanol 

78.3 

Methyl 

n-Propyl  Ketone 

102.0 

Propanol 

97.2 

Methyl 

i-Butyl  Ketone 

115.9 

Water 

100.0 

Methyl 

n-Amyl  Ketone 

150.6 

Butanol 

117.8 

Amyl  Alcohol 

137.8 

Ethylene  Glycol 

197.2 
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If  the  concept  of  molecular  similarity  is  insuf¬ 
ficient  in  the  choice  of  a  solvent,  two  semi-quantitative 

31 

methods  are  available.  Ewell  et  al  divided  all  liquids, 
depending  on  their  hydrogen-bond  forming  capabilities,  into 
five  classes.  Since  any  reduction  of  the  intermolecular  bonds 
in  a  system  will  lead  to  more  mobility  and  a  higher  volatility 
of  its  components,  the  direction  and  approximate  extent  of 
the  deviation  from  Raoult's  law  caused  by  mixing  any  two  of 
the  classes  of  liquids  could  be  predicted  from  the  classifica¬ 
tion.  However,  the  hydrogen  bond  is  only  one  of  the  inter¬ 
molecular  forces  that  could  exist  in  a  system.  Consequently, 
this  classification  fails  to  predict,  for  example,  the  in¬ 
creased  volatility  in  aromatic  hydrocarbons  when  a  paraffin 

is  added,  even  though  both  fall  in  the  class  of  liquids  hav- 

.  19 

mg  no  hydrogen  bond-forming  capabilities 

2  i 

Pierotti  et  al  offer  a  method  for  predicting  devi¬ 
ations  from  Raoult's  law  more  quantitatively  than  is  possible 
with  Ewell's  classification.  The  method  relates  the  change 
in  the  excess  partial  molar  free  energy  of  a  solute  at  an 
infinite  dilution  in  a  solvent  to  a  series  of  interaction 
constants.  Since  excess  partial  molar  free  energies  are 
directly  relatable  to  activity  coefficients,  it  becomes  pos¬ 
sible,  through  Pierotti 's  correlation,  to  predict  the  in¬ 
fluence  of  the  molecular  structure  of  solvent  or  solute  on 
the  volatilities  of  the  components  in  a  system. 
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3 . 4  Solvent  Concentration 

An  extension  of  both  the  molecular  similarity  con¬ 
cept  and  the  hydrogen-bond  concept  would  be  that  increasing 
the  solvent  concentration  in  a  liquid  mixture  could  enhance 
the  volatility  of  the  solute.  This  is  found  experimentally 

to  be  true.  Classical  data  in  this  respect  are  those  of 

32 

Souders  et  al  .  Figure  3.1  shows  the  effect  of  varying  the 
solvent  concentration  on  the  ratio  of  non-toluene  to  toluene 
in  the  vapor-phase  of  such  mixtures  (on  a  solvent-free  basis) . 
The  non-toluene  mixture  is  a  narrow  cut  from  a  straight-run 
naphtha.  These  data  show  that  as  the  concentration  of  the 
extractive  agent  increases,  the  column  required  for  a  given 
separation  will  be  shorter  when  all  other  specifications  are 
kept  constant.  On  the  other  hand,  the  increase  in  the  sol¬ 
vent  will  require  increased  heat  loads  and  column  diameters, 
so  that  an  economic  balance  can  be  struck  between  the  advan¬ 
tages  and  disadvantages  of  increased  solvent  concentration. 

3 . 5  Description  of  an  Extractive  Distillation  Column 

A  typical  extractive  distillation  column  will 

consist  of  three  sections.  The  function  of  each  section  can 

30 

be  thought  of  as  follows  :  the  section  between  the  top  of 
the  column  and  the  solvent  introduction  point  is  the  ''sol¬ 
vent  recovery  section'5 .  It  removes  the  solvent  from  the 
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Figure  3.1  Vapor-Liquid  Equilibrium  for  Mixtures 
of  Toluene-Non-Toluene 
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vapors  going  up  the  column.  The  section  between  the  solvent 
and  the  feed  introduction  points  removes  from  the  vapor 
phase  the  heavier  components  and  is  known  as  the  'absorbing 
section".  The  "stripping  section",  which  is  below  the  feed, 
removes  the  light  component  from  the  liquid.  The  bottom  pro¬ 
duct  from  a  distillation  is  essentially  a  binary  containing 
the  heavy  key  and  the  solvent.  Thus,  a  solvent  recovery 
column  is  usually  needed  to  purify  the  solvent  for  recycle 
into  the  extractive  distillation  column.  A  typical  extrac¬ 
tive  distillation  flowsheet  is  shown  in  Figure  3.2 


3 . 6  Study  of  the  Design  Variables  in  Extractive  Distillation 

A  two-feed  column  with  a  total  condenser  and  a 

19 

partial  reboiler,  has  2C  +  2N  +  12  degrees  of  freedom  ,  where 
C  denotes  the  number  of  components  in  the  system,  and  N  the 
number  of  equilibrium  stages.  The  choice  of  the  2C  +  2N  +  12 
specifications  depends  on  the  purpose  of  the  study,  the  items 


available  for  specification  being  the  following: 

Pressure  in  N  stages  N 

Pressure  in  condenser  1 

Pressure  and  heat  leak  in  reflux  divider  2 

Heat  leak  in  N-l  stages  N  -  1 

Reflux  temperature  or  condenser  load  1 

Total  number  of  stages  1 

Fresh  feed  composition,  temperature  and  rate  C  +2 

Solvent  feed  composition,  temperature  and  rate  C  +2 
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Figure  3.2  Simplified  Flow  Diagram  for  an 
Extractive  Distillation  Process 
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Fresh  feed  location  1 

Solvent  feed  location  1 

Rate  of  overhead  product  1 

One  reflux  ratio  1 

Overhead  concentration  of  the  extractive  agent  1 

Ratio  of  light  key  to  heavy  key  components 

in  bottom  product  1 

Recovery  of  one  of  the  keys  1 

2C  +  2N  +  15 

The  two  operating  variables  that  most  strongly  affect 
the  performance  of  an  extractive  distillation  column  are  the 
reflux  ratio  and  the  ratio  of  the  solvent  feed  rate  to  the 
fresh  feed  rate.  However,  before  these  two  variables  can 
be  studied  and  optimized,  other  variables  relating  to  the 
physical  column  configuration  such  as  the  feed  locations,  the 
total  number  of  stages,  etc.  have  to  be  fixed.  These  can  be 
determined  from  a  simulation  of  the  column  if  it  exists  al¬ 


ready,  or  likely  random  values  for  these  variables  can  be 
specified . 
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IV.  GENERALIZED  MODEL  FOR  DISTILLATION 
4 . 1  The  Pertinent  Equations 

The  purpose  of  any  rigorous  calculation  method  in 
distillation  design  is  the  determination  of  the  unique  set 
of  flow  rates,  temperatures  and  compositions  that  will  exist 
in  the  column.  To  obtain  such  a  unique  solution  a  necessary 
and  sufficient  number  of  variables  has  to  be  specified  in 
advance.  The  basic  concept  used  in  all  distillation  models 
is  the  equilibrium  stage  shown  in  Figure  4.1  where  the  vapor 
and  liquid  streams  leaving  the  stage  are  assumed  to  be  in 
equilibrium . 


The  equations  relating  to  this  model  can  be  expressed 


as  follows: 


1)  The  component  i  material  balance  around  stage  j 


+  F  .Z  .  . 
D  ij 


0 


(4.1) 


2)  The  equilibrium  relationship 


0 


(4.2) 


3)  The  restriction  on  fractional  concentrations 
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L  ID 


(4.3) 
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Figure  4.1  The  Model  Column  and  the  Equilibrium  Stage 
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4)  The  energy  balance  around  atage  j 


L.  ,h.  .  -  (V.  +  W . )  H  .  -  (L.  +  U  . )  h  . 
D-l  D'1  3  3D  3  33 


+  V  .  ,  ,  H  .  , +  F  .H_  -  Q  .  =  0 

j+l  :+l  3  F_.  ^ 


(4.4) 


Equation  (4.2)  eliminates  the  necessity  of  dealing  with  vapor 

3  6 

concentrations  in  the  other  three  equations  .  Also,  the 
overall  material  balance  equation 


V . 
3 


+ 


k=2 


D 


2  <_  j  <_  n-1 


(4.5) 


can  be  used  to  eliminate  one  set  of  the  phase  rates  V.  or  L.. 

3  3 

Thus,  the  above  equations  4.1,  4.3  and  4.4  can  be  respectively 
represented  as  a  material  balance  equation  M,  a  sum  equation 
S,  and  an  energy  equation  E  as  follows. 


ij  ' 

V. 

3 

or 

Lj, 

V 

=  0 

(4.6) 

ij  ' 

V  . 

3 

or 

Lj, 

V 

=  0 

(4.7) 

ij' 

V  . 

3 

or 

Lj, 

V 

=  0 

(4.8) 

This  system  of  equations  is  consistent  since  the  number  of 

independent  variables  is  equal  to  the  number  of  equations. 

Any  rigorous  computational  procedure  aims  at  finding  the  set 

of  x. .,V.  and  T.  that  will  satisfy  these  equations. 
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The  M-equation  provides  a  liquid  composition  vector 
when  the  and  T_.  profiles  are  available.  A  detailed  pre¬ 
sentation  of  the  methods  available  for  the  solution  of  this 
equation  is  given  by  Smith^.  More  recently,  Wang  and  Henke^ 
proposed  an  algorithm  for  the  solution  of  the  M-matrix  that 

does  not  involve  the  subtraction  of  nearly  ecrual  numbers. 

2  6 

Thus,  it  preserves  accuracy  and  is  numericall/  stable 

Equation  (4.7)  is  solved  for  the  temperature  vector 
in  the  column.  The  solution  of  this  equation  involves  a 
trial  and  error  procedure  to  find  the  temperature  at  every 
stage  that  satisfies  one  of  the  conditions  of  equation  (4.3) 
(dew  point  or  bubble  point  temperatures  respectively) .  The 
E-equation  is  the  energy  balance  equation  which  provides  the 
flow  rates  in  the  column.  Except  for  the  case  of  high-boiling 
components,  direct  algebraic  solution  of  equation  (4.8)  is 
found  to  be  quite  adequate. 

The  approach  used  in  this  work  to  solve  equations 
(4.6)  to  (4.8)  is  outlined  in  the  following  sections. 

4 . 2  Solution  of  the  Material  Balance  Equations 

Figure  4.1  also  represents  the  model  column  used  to 

derive  the  working  equations. 

The  column  has  N  equilibrium  stages.  A  partial  or 

total  condenser  is  the  first  stage,  and  the  reboiler  is  the 
t  h 

N  stage.  The  model  is  derived  for  the  most  complex 
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conditions  possible:  each  stage  includes  a  feed  stream,  F ^ , 
one  vapor  and  one  liquid  side  stream,  and  U ^  respectively, 
and  one  heat  exchange  device  with  load  . 

To  derive  the  solution  algorithm,  we  rewrite  the 
material  balance  equation  (4.1)  after  the  substitution  for 
y .  .  from  (4.2)  : 

JT  1-j 

L.  .x.  .  ,  -  (V.  +  W . )  K  •  .x.  .  -  (L.  +  U  . )  x . 


+  V  .  ,  K  .  ,±1x.  .  .  +  F . Z .  .  =  0 

3+1  1,3+1  1,3+1  3  13 


(4.9) 


Collecting  terms 


L  .  n  x  .  ■  -» 

3-1  i/3-l 


(V.  +  W.)K.  .  +  (L.  +  U.) 
LI  3  13  3  3 


x .  . 

13 


+  V .  , , K .  ,±1x,  .  , - 

3+1  1,3+1  1,3+1 


Defining 


B  . 
1 


L  3 


3  13 


C.  =  V  ,  ,K.  ... 

3  3+1  i/D+1 


-F  .  Z  . 

3  13 

(4.10) 

+  (L  .  +  U  .  ) 

3  3  J 

2  <_  j  n-1 

(4.11) 

2  s.  j  <.  n-1 

(4.12) 

A  . 
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D  . 
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=  L 


j-1 


jri 
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-F  .  Z .  . 
3  13 


2  <_  j  <_  n-1 
2  i  j  <  n-1 
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(4.13) 
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and  substituting  into  equation  (4.10) 


A  .  x .  .  ,  +  B . x .  .  +C.x.  . , , 

i  i,:-i  i  id  :  1,1+1 


D  . 

1 

2  <_  j  <_  n-1 


(4.14) 


This  equation  now  lends  itself  to  matrix  notation  as  follows: 
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A  . 
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xu 

C2 

xi2 

B  .  C  . 

X  .  . 

3  3 

11 

A  ,  B  .  C  , 

X  .  , 

n-1  n-1  n-1 

i ,  n-1 

A  B 

x . 

n  n 

in 

D. 


D. 


D  . 
3 


D 


n-1 


D 


n 


(4.15) 


Thus,  the  material  balance  equations  in  the  column 

reduce  to  a  tridiagonal  matrix  notation.  One  algorithm  to 

2  6 

solve  for  {x}  is  offered  by  Wang  and  Henke  ,  derived  from  a 
Gauss  elimination  method.  In  this  method  of  solution,  two 
auxiliary  quantities  p_.  and  q..  are  first  calculated  across 
the  column: 

Pi 


B. 


D1 

ql  =  _ 

B, 


(4.16) 
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With  the  vectors  { p ^  >  and  {q ^ }  available,  the  liquid  composi¬ 
tion  at  the  bottom  stage  is  calculated  first,  then  j  is 

decreased  till  x.,  is  reached.  Thus, 

ll 


x . 


iN 


=  q 


N 


(4.19) 


x  .  .  = 

13 


q.  -  p.x.  1  j  n-1  (4.20) 

3  3  1.3+1 

This  algorithm  provides  a  simple  method  for  the  solution  of 
the  matrix  equation  (4.15). 


4 . 3  Determination  of  the  Temperature  Profile 

When  the  composition  of  a  liquid  and  the  equili¬ 
brium  constants  of  its  components  are  available,  the  prob¬ 
lem  is  to  find  a  value  of  T.  such  that  the  second  condition 

3 

in  equation  (4.3)  is  satisfied.  Several  approaches  are  used 
in  the  determination  of  this  bubble-point  temperature. 

One  approach  is  to  use  Newton's  method.  Rearrang¬ 
ing  equation  (4.3)  in  functional  notation  gives 

C 

f (T)  =  y  K.x.  -  1.0  (4.21) 

i=i  1  1 

Thus,  the  positive  value  of  T  is  to  be  found  such  that 
f (T)  =  0.  Since  each  of  the  ' s  increases  with  temperature, 
f (T)  has  only  one  root. 

If  T  is  the  assumed  value  of  T  for  trial  number 
n 

n,  the  temperature  used  in  the  next  iteration,  T 
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by  Newton's  formula 


n+1 


where 


f'(T)  = 


T  - 
n 


f  (T  ) 
n 

f '  (T  ) 
n 


(4.22) 


C  dK. 

I  *i(— ) 


(4.23) 


i=l 


dT 


=  the  first  derivative  of  f (T) 

Thus,  Newton's  method  consists  in  drawing  successive  tangents 
to  the  f (T)  vs.  T  curve.  Convergence  with  this  method  is 
not  guaranteed.  Depending  on  the  initial  guess  of  T  and  the 
particular  shape  of  f (T) ,  convergence  may  or  may  not  occur. 

If  f  (T)  has  turning  points  or  inflections  in  the  interval  of 

interest,  Tn+]_  maY  further  away  from  the  desired  root  than 

m  37 
Tn  * 

Another  disadvantage  of  this  method  is  the  need 

to  determine  (dK^/dT)  which  is  almost  impossible  if  is  a 

complex  function  of  temperature,  pressure,  and  composition 

3  •  4 

such  as  in  Chao-Seader  or  Pfausnitz  models. 

Another  approach  for  the  determination  of  bubble- 
point  temperatures  is  the  regula-falsi  method.  Here,  two 
arbitrary  values  of  T  are  selected.  Denoting  these  by  and 
Tn+^,  the  next  best  value  Tn+2  can  k>e  obtained  from 


n+2 


T  f(T  ,-.)  -  T  -  f(T  ) 
n  n+1_ n+1  n 

f (T  .  )  -  f (T  ) 
n+1  n 


(4.24) 
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This  technique,  however,  is  inefficient  and  again  convergence 
is  not  guaranteed. 

2  6 

Wang  and  Henke  report  good  success  in  bubble- 

27 

point  temperature  determinations  using  Muller's  method 

This  is  the  method  that  has  been  used  in  this  work  exclusively. 

Muller's  method  consists  in  making  three  arbitrary  estimates 

The  next  estimate  is  based  upon 


—  29 

of  the  root  desired  (T) 


fitting  a  second-degree  Lagrangian  interpolation  formula  to 


the  points 


f  (T  (o)  )  ,  T  (o) 


f (T (1) ) ,T (1) 


,  and 


f (T (2) ) ,T (2) 


and  finding  T^,  where  T^°^,  T^,  and  T^  are  the  estimated 


roots.  The  next  iteration  then  uses  T 


(3)  ,  T (1)  and  T (2)  as 


the  estimates  and  T 


(4) 


becomes  the  new  root,  etc.  The  itera- 


T(n)  _  T (n-1) 


<  e  where  e 


tion  process  is  continued  until 
is  a  predetermined  tolerance  in  the  value  of  the  root. 

The  equation  determining  the  choice  of  the  next 

root  is 


Tk  =  Tk_x  +  (Tk_1  -  Tk_2)*dk  k>4 


(4.25) 


where 


-2*f(Tk_1)*(i  +  a  ) 


t  /b2  -  4*f  (Tk-1)*dk-1*(1  +  dk-l)*C 


(4.26) 


k>  4 


In  equation  (4.26) 


C  =  f<Tk-3>*dk-l  -  f(Tk-2>*(1  +  dk-l>  +  f(Tk-l>  (4'27) 


frs  - 


90n9?i*vnoo  niePB  baa  inalolll^l  el  ,«v8wort  ,90pinrk>9l  mid* 


. bse^nexeup  :*on  si 


rSbortl9m  a'xalluM  pnieu  enollBniinslsb  aiuitroqma*  Jnxoq 


.  ^ (T)  be^iseb  3ooi  9rfi  10 


oi  Biurmoi  aoij»A«4w»«*  - - 

T.(‘£,T)»tb«.  .PI>*,((I,*)»] 

sralle9  ex.  <«*  bn.  ,<X>*  vto'i  -eriw  ,(  T  pnlbnl!  bn. 


T  DflB  T  y 


T  U&&U  aerU  r.ottiisa*!  2xo:i  ®rT  •  — 


(fr) 


T  bns  89^®mx^89  sri* 


,  (X-n)T  _  (n '  »j»  lx:*flu  bsunxinoo  si  eseooxq  noil 

i 


.  joox  sriJ  lo  «oI.v  9ril  nl  sonexs ,oi  bsniraxs isbexq  .  si 
jxsn  sril  lo  aolorto  srfJ  pnlnlmxexeb  nollsup*  e-1? 


el  loox 


(ds.fr) 


eidrfw 

triid 


(dS.fr)  noi^fli/p®  al 


25 


and 

b  =  f(Tk-3>*(dk-l)2  '  f(Tk-2>*(1  +  dk-l>2 

+  f(Tk_i)*(1  +  2*dk_x)  (4.28) 

Equation  (4.26)  will  yield  two  roots  due  to  the  + 
notation  in  the  denominator.  If  one  of  these  roots  is  nega¬ 
tive  it  is  immediately  rejected.  If  both  roots  are  positive, 

N 

the  root  that  yields  a  value  of  £  y.  closer  to  1.0  is 

i=l  1 

accepted . 

Equation  (4.25)  is  iterated  upon  until  the  tempera¬ 
ture  that  satisfies  equation  (4.3)  is  determined.  No  lack 
of  convergence  has  been  experienced  with  this  method.  Of 
importance  in  this  method  is  that  the  starting  values  for 
the  iteration  are  not  special  and  f'(T)  need  not  be  evaluated 
The  method  is  graphically  presented  in  Figure  4.2 


4 . 4  Solution  of  the  Energy  Balance  Equations 


The  conventional  method  of  the  solution  of  the 
energy  balance  equations  has  been  exclusively  followed  in 
this  work.  An  energy  balance  around  the  equilibrium  stage 
shown  in  Figure  4.1  can  be  written  as  follows: 


L.  .h.  ,  -  (V.  +  W.)H.  -  (L.  +  U . )  h . 
j-1  2-1  j  23  3  33 


+  V.  ,,H  .  .  +  F  .  H  -  Q.  =  0 

2+1  3+1  3  Fj  3 


(4.4) 
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S  . 
D 


Newton-Raphson 


Figure  4.2  Graphical  Comparison  of  Muller's  Method 
to  Newton-Raphson 
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rearranging  terms  and  substituting  for  L ^  from  a  material  balance 

(H .  ,  -  h  . )  V  .  ,  -  (H.  -  h.)(V.  +  W.) 

1+1  3  D+l  3  33  3 


-  (h.  -  h.  , ) L .  .  +  F . (H  -  h.)  -  Q.  =  0  (4.29) 

3  D-l  1-1  3  F.  j' 


Equation  (4.29)  is  used  for  the  determination  of 
the  vapor  rate  profile  across  the  column.  A  material  balance 
around  each  stage,  such  as  with  equation  (4.5) ,  provides  the 
liquid  rates. 


4 . 5  Order  of  Solution 

Equations  (4.6)  to  (4.8)  inclusive,  are  solved  in 

t  he  following  order  using  the  computational  procedures  out¬ 
lined  above: 

1)  Assume  an  initial  profile  for  the  vapor  phase  and  a 
temperature  profile  across  the  column  . 

2)  Calculate  by  any  of  the  equilibrium  equations  pre¬ 

sented  in  Chapter  V. 

3)  Solve  the  M-equation,  equation  (4.6),  using  the  Gauss 
elimination  algorithm  (equations  (4.16)  to  (4.20))  to 
obtain  the  liquid  composition  profile  in  the  column. 

4)  Solve  the  S-equation,  equation  (4.7),  to  obtain  the 
temperature  vector  in  the  column. 

5)  Calculate  the  enthalpies  of  the  streams  by  one  of  the 
methods  presented  in  Chapter  V. 
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6) 

7) 


Solve  the  E-equation  for  a  new  set  of  . 

Repeat  steps  2  to  6  until  the  sum  of  the  squares  of 
the  stage  temperatures  between  any  two  successive  itera¬ 
tions  is  less  than  a  certain  value,  usually  0.01  N, 
where  N  is  the  number  of  stages. 


V  lo  198  W9fl  «  lOl  n0il60p9-3  »rfl  9Vi03 
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V.  CALCULATION  OF  EQUILIBRIUM 
AND  ENTHALPY  DATA 

The  solution  of  the  material  balance  equations  in 
a  column  requires  equilibrium  data  for  the  components  being 
distilled.  Enthalpy  data  on  the  phases  in  eauilibrium  are 
necessary  for  the  energy  balance  equations.  The  methods 
that  are  in  use  to  obtain  such  data  will  be  discussed  in 
this  chapter. 

5 . 1  Polynomial  Presentation  of  K  and  H  Data 

In  order  to  render  the  large  quantity  of  experi¬ 
mental  equilibrium  and  enthalpy  data  available  to  a  computer, 
it  is  necessary  to  represent  the  data  in  equation  form. 
Several  curve-fitting  techniques  have  been  tried  for  this 
purpose,  with  polynomial  presentation  being  most  suitable 
due  to  its  ease  of  handling.  Examples  of  this  form  of  data 
presentation  is  the  polynomial  presentation  of  the  equilibri¬ 
um  ratios  of  twelve  hydrocarbon  components  originally  pre- 

5 

sented  graphically  in  the  NGAA  data  book  .  These  charts 

correlated  K-ratios  in  terms  of  temperature  and  pressure, 

the  effect  of  composition  being  taken  into  account  by  the 

additional  parameter  known  as  the  convergence  pressure. 

Another  example  of  this  form  of  data  presentation 

17 

is  the  Yen  and  Alexander  enthalpy  correlation,  whereby 
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enthalpy  departures  of  some  hydrocarbon  components  can  be  cal- 

6 

culated  from  a  polynomial  curve-fit  of  Rossini's  data  using 
Forsythe  polynomials. 

Although  with  increasing  sophistication  in  curve¬ 
fitting  techniques,  data  presentation  by  this  method  might 
gain  wider  applications,  some  disadvantages  are  inherent  in 
them,  the  most  serious  being  the  discreteness  of  the  condi¬ 
tions  for  which  data  can  be  retrieved.  For  example,  K-ratios 
can  only  be  obtained  at  definite  convergence  pressures,  and 
H-data  only  within  defined  limits  of  reduced  properties. 

As  a  result  of  this  inability  to  "predict"  data,  new  informa¬ 
tion  is  needed  every  time  a  problem  is  run  at  new  conditions. 


5 . 2  The  Chao-Seader  Correlation  for  Equilibrium  Ratios 

3 

The  correlation  offered  by  Chao  and  Seader  for 
the  calculation  of  K  ratios  for  hydrocarbons  of  various 
molecular  types  is  amenable  to  machine  computation.  K's  are 
calculated  by  combining  three  factors  as  follows: 


K. 
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y  .  y  . 
li 


(5.1) 


In  this  equation,  is  the  f/p  ratio  in  the  liquid  phase 

is  the  activity  coefficient  in  the  liquid 
phase 


<Jk  is  the  fugacity  coefficient  in  the  vapor 
phase 
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vu  is  calculated  by  combining  a  term  for  the  fugacity  coef¬ 
ficient  of  a  simple  fluid  with  a  correction  factor  for  real 
fluids  as  follows: 


log 


log  yR  +  a)  log  y3"^ 


(5.2) 


where 


log  y? 
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1.22060 


4.23893  -  8.65808T  - 
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(5.3) 


-  3.15224T  -  0 . 025 (P  -  0.6)  (5.4) 

r  r 

The  constants  used  in  equation  (5.3)  are  specific  for  methane, 

nitrogen,  hydrogen,  carbon  dioxide  and  hydrogen  sulfide  and 

common  for  all  other  components  handled  by  this  correlation. 

Oxygen  and  carbon  monoxide  use  the  same  constants  as  methane. 

7 

These  constants  were  tabulated  by  Erbar  .  (Appendix  III) 
y ^  is  calculated  from  the  following  equation: 
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This  is  the  Scatchard-Hildebrand  model  and  assumes  hydrocar¬ 
bon  solutions  to  be  regular  i.e.  to  have  no  excess  entropy  of  mixing. 

The  determination  of  <j>^  uses  the  Redlich-Kwong 
equation  of  state  in  the  form 


In  <j> . 
1 


where 
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and 


B . 
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(5.7) 

(5.8) 

(5.9) 

(5.10) 

(5.11) 

(5.12) 

(5.13) 


Since  (5.8)  is  implicit  in  Z  it  is  solved  by  trial-and- 
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5 . 3  Chao-Seader  Correlation  for  Enthalpy  Calculations 

1 8 

Edmister  et  al  have  extended  the  Chao-Seader 
correlation  to  the  calculation  of  vapor  and  liquid  enthalpies 
in  hydrocarbon  mixtures. 

The  vapor  phase  enthalpies  are  calculated  from 

3  A  2  BP 

- E_  in  ( 1  +  -£L_)  +  1  -  Z  (5.14) 

2  B  Z 


H 


m 


=  I 


y.H.°  + 

1 1  l 


The  first  term  represents  the  enthalpy  of  a  gas  mixture  at 
the  ideal  gas  state  pressure  of  zero.  The  second  term  re¬ 
presents  deviations  from  ideality  due  to  pressure. 

Liquid  phase  enthalpies  are  calculated  from  the 
relationships 

h  =  J  h.  (5.15) 
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and 


3  In  y . 

( - -) 

3T 


P,x. 


3 . ) 2  x  1.8 

l 


(5.18) 


H^°,  the  ideal  gas  state  enthalpy  of  a  component, 
is  calculated  from  a  cubic  equation  in  temperature  as  follows: 


H.°  =  AH.  +  (BH.)T  +  (CH.)T2  +  (DH.)T3  (5.19) 

l  l  i  l  l 

where  the  coefficients  are  tabled  for  all  the  components  that 
can  be  handled  by  the  correlation. 


5 . 4  The  Prausnitz  Equilibrium  Model  Based  on  Wilson's  Equation 

In  1964  Wilson^ proposed  a  model  that  could  calcu¬ 
late  liquid  phase  activity  coefficients  in  multicomponent 
systems  using  only  parameters  obtainable  from  the  binary 

4 

data.  Orye  and  Prausnitz  described  a  vapor-liquid  equili¬ 
brium  model  in  which  they  proposed  to  incorporate  the  Wilson 
equation.  This  equilibrium  model  was  further  developed  by 
O'Connell  et  alxu  and  Lyckman  et  al  ,  and  is  presented  in 
this  section. 

For  equilibrium  between  a  vapor  and  a  liquid  phase, 
the  condition  to  be  satisfied  is  the  equality  of  the  fugaci- 
ties  of  the  two  phases,  i.e. 


(5.20) 
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These  terms  can  be  developed  as  follows^  : 


P  v . L 


f  . 
1 


=  Y 


c  o(po)  ,  i  N 

.x.f  .  exp  ( - ) 

RT 


ill 


(5.21) 


and 


f  . 

l 


V 


=  4> .  y  .  P 

l  J  l 


(5.22) 


Substituting  in  (5. 20) yields  the  equilibrium  equation  pre¬ 


setted  by  Orye  and  Prausnitz  thus: 

f  o(po) 


<f> .  y  .  P  =  y  .  x 
ri  ill 


v.LP 
exp  (  — — ) 
RT 


(5.23) 


From  this  equation 


,  x  v . LP 

<f,°(p0)  exp  f-^-) 


1  1 


Ki  = 


(5.24) 


X 


<t>.p 

1 


The  meaning  and  computation  of  each  term  will  be  developed 
separately . 


The  Fugacity  Coefficient  of  the  Vapor-Phase  (<}>^) 


This  term  can  be  calculated  from 


2  NCP 


In  . 

l 


=  Z  I 

v  j  =  l 


y.B.  . 
Ji  i: 


-  InZ 


(5.25) 


where 

Z  =  compressibility  factor 

v  =  molar  volume  of  the  vapor  mixture 

B  =  second  virial  coefficient  of  the  i-j  inter¬ 


action 
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Z  and  B  are  related  by 

NCP 

(.1 

Z  =  1  +  -1—1- 


NCP 

V 


(5.26) 


The  virial  coefficient  is  calculated  from  Pitzer  and  Curl's^ 

correlation  for  non-polar  gases.  This  correlation  was  re- 

13 


ported  by  Lewis  and  Randall 


For  polar  gases  a  recent  pub¬ 


lication  by  O'Connell  and  Prausnitz  y  offers  a  correlation 
based  on  an  extended  theory  of  corresponding  states. 

The  Reference  Fugacity  in  the  Liquid  Phase  f  °(P°) 

f.°(P°)  is  the  fugacity  of  the  saturated  liquid 
component  i  at  zero  pressure.  It  is  defined"*^  as 


where 


f  o(po) 

l 


-  L_  s 
v .  P 


,  s  .  s  ,  1  1  x 

’  •  4>  •  exp  (-  - ) 

RT 


(5.27) 


P^  is  the  saturation  vapor  pressure  at  T 

s  s 

is  the  vapor  fugacity  coefficient  at  P^  and  T 

v.L  is  the  molar  liquid  volume  at  T 
23 

Lyckman  defines  the  quantity  ip  such  that 


ip  = 


f  o(po) 
i 


(5.28) 


where  P  is  the  critical  pressure  of  i.  A  correlation  is 
ci 

then  presented  where  ij;  is  a  function  of  the  reduced  tempera¬ 
ture  and  the  acentric  factor  of  the  component  w .  Pitzer 
suggests  the  following  form  for  ijj; 


-dxjq  in9OM  6  MMt'-uloq  w*  •  '‘iW"*  6ne  ,iW*J  *d  bS^0q 
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log  ip  =  log  (\p)°  +  co  log  (i^)^  (5.29) 

where  log  (\p)°  and  log  (i|>)^  are  each  functions  of  the  reduced 

g 

temperature.  Lyckman  et  al  present  this  functionality  in 
table  form  and  they  have  been  polynomially  fitted  vs.  T  for 
purposes  of  this  work  as  follows: 


and 


log  (i|/) 


-14.8016  +  49.4529  T  -  71.3050  T  ‘ 

r  r 

+  54.1713  T  3  -  21.0864  T  4 
r  r 


+  3.3213  T 


(5.30) 


variance  =  +0.56048  x  10 


-5 


log  (\p) 


=  -54.2280  +  409.1308  T  -  1538.1941  T 


variance  = 


+ 

3445.8673  T  3 

r 

-  4812.2940 

i-3 

+ 

r 

4221.5419  T  5 

-  2259.4562 

r 

T  6 

+ 

r 

674.3256  T  7  - 

86.0182  T 

r 

8 

+  0 

r 

.16503  x  10"4 

r 

(5.31) 


These  variances  being  the  minimum  in  all  polynomial 

t  h 

fits  up  to  and  including  the  10  degree. 

8  s 

Lyckman  also  offers  a  correlation  for  <f> ^  where 

s 

<j>^  is  a  function  of  and  co ,  and  is  represented  by 


log  <p±s  =  log  U^)0  +  co  log  U^)  1  (5.32) 

The  table  of  values  of  (4>^s)°  and  (<f>^S)'*'  presented  by  Lyck- 
23 

man  have  also  been  polynomially  fitted  against  T^  and  the 
relationships  are  as  follows: 


nl  Ylil.noliwu.1  sxriX  Xneesxq  I.  *•  .e™*”**™* 


T  oeoe.xv  -  ' *  ««.«  ♦  -  °<*>  ^  , 

C  T  £!££.£  + 


e  t  exxe.t  + 
.  ' 


'■<♦)  pol 


f 


d  T  Sd2*.G2£S  - 


*"oi<  x  £oeai.o+ 


iaxraonYloq  lie  <U  apmLatm  arid  pnxsd  asor.sxiav  ee«IT 


' 

' 

...  -  • 


srii  f*u>  T  *snx6i?«s  baidXi  YXX6xmon><Xoq  nosd  oeXs  evert  nsm 


a  :o  lol  fcfi  ***  eqirfi  novels* 


38 


(*±3)° 


53.15801  -  433.7522  T  +  1463.0379  T  ‘ 

r  r 

-  2611.1941  T  3  +  2601.3246  T  4 

r  r 

-  1372.3205  T  5  +  299.57934  T  5 

r  r 


variance  =  +0.6874*10 


-6 


(5.33) 


and 


U+)1 * 


-0.16186  +  2.98895  T  -  14.05742  T  ' 

r  r 

+  27.77032  T  3  -  24.06024  T  4 
r  r 

+  7.48698  T  5 


(5.34) 


variance  =  0.6299*10 


-6 


Combining  equations  (5.23)  and  (5.27),  and  assuming 
the  equality 


-  L  L 

v .  =  v . 

l  l 


(5.35) 


we  obtain 


°  Yixipis*is  exP  ^ 


(P  -  p.°) 


(5.36) 


Liquid  Phase  Activity  Coefficient  y^ 

14 

Wilson  presented  a  model  where  the  excess  free 
energy  of  mixing  is  a  logarithmic  series  in  the  liquid  compo¬ 

sition  and  in  parameters  obtainable  from  the  binary  equili- 

4 

bna  only.  Orye  and  Prausnitz  derive  the  model  through  the 
following  steps: 
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At  constant  temperature  and  pressure,  an  exact 

E 

relationship  exists  between  g  ,  the  excess  Gibbs  energy  of 
mixing  and  y^,  the  activity  coefficient  for  any  component  k 
in  a  liquid  mixture.  This  relationship  is 


3n  *< 


E 


RT  In  yk  =  (• 


•) 


T , P , all  n. (i^k) 
dnk  1 


(5.37) 


where  n^  =  number  of  moles  of  component  i 

nT  =  total  number  of  moles 

From  the  concept  of  athermal  solutions  developed  by  Flory 

,  „  16  E 

and  Huggins 


15 


g  can  be  expressed  as  follows: 
E 


g 

RT 


NCP 

=  l  x  m 

i=l  1 


*i. 


x . 

l 


(5.38) 


where 


x .  v . 

l  l 


NCP 


i=l 


=  volume  fraction  of  i  (5.39) 


x .  v . 

l  l 


and 


v 


=  molar  liquid  volume  of  i 

JL 

From  basic  statistical  concepts,  in  a  binary  solu¬ 
tion  of  components  1  and  2,  the  probability  of  finding  a  type 
2  molecule  around  a  central  molecule  of  type  1  as  compared 
to  that  of  finding  another  type  1  molecule  is 


12 


11 


,  21, 

*2  eXp  ~  ~RT 

,X11. 

X.  exp  -  ( - ) 

RT 


(5.40) 
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where  A^  and  A  are  proportional  to  the  1-2  and  1-1 
interaction  energies 

and  the  exponential  terms  are  Boltzmann  factors. 

With  this  concept  in  mind,  the  volume  fractions  in 

14 

the  Flory-Huggins  equation  are  redefined  into  what  Wilson 
terms  local  volume  fractions 


L  A 11 

x1v1  exp  -  (-^p— ) 


1  L  / 11,  .  L  ,X12, 

Xnv,  exp  -  ( - )  +  X9v?  exp  -  ( - ) 

-1  RT  RT 


(5.41) 


and 


L  .22, 

x^v^  exp  -  (p^r~) 


‘2  2 


L  .22,  ,  L 

x„v„  exp  -  ( - )  +  x.v,  exp 

RT  1  1 


‘2  2 


^21 
(— ) 

RT 


(5.42) 


defining 


,  v2  ,X12  Xll. 

A12  =  —  exp  -  ( - ) 


V- 


RT 


(5.43) 


and 


A 


21 


V1  ,X21  X  22 , 

— f  exp  -  ( - ) 

v0  RT 


(5.44) 


and  substituting  into  (5.41)  and  (5.42)  then 


X1  +  x2A12 


and 


x2  +  x1A21 


(5.45) 


(5.46) 


substituting  these  local  volume  fractions  for  <j>^  in  equation 
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(5.38)  , 


where 


and 


g 

RT 


A 


Di 


NCP 


?  x .  In  (  y  x.A..) 
i-1  1  i=l  3 


NCP 


V  . 

A.  .  =  -1 

^  v.! 

1 


V  . 
1 


v  . 
3 


X  .  .  -  X  . 

exp  -  (-iJ - — ) 

RT 


X  .  .  -  X  .  . 

exp  -  (— - -LI) 

RT 


(5.47) 


(5.48) 


(5.49) 


substituting  equation  (5.47)  into  (5.37),  a  general  expres¬ 
sion  for  the  activity  coefficient  can  be  written  as  follows' 


In  y. 


NCP 

In  (  l  x.A  )  +  1 
j=l  3 


NCP  x.  A., 
v  i  lk 

.  L  NCP 
1=1 


(5.50) 


y  x .  a  . . 

i=l  3  13 


Thus,  the  liquid  activity  coefficient  of  every 
component  in  a  mixture  can  be  obtained  from  the  liquid  com¬ 
position  and  parameters  obtainable  from  the  respective 
binary  equilibria.  The  calculation  of  these  parameters  will 
be  discussed  in  the  following  chapter. 

From  the  above  information  then,  the  following 
steps  are  taken  to  obtain  vaporization  constants  from  equa¬ 
tion  (5.24)  : 
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1)  y^'s  at  any  liquid  composition  are  obtained  from  the 
Wilson  equation  (5.50)  usinq  the  Wilson  parameters. 

The  following  chapter  discusses  means  of  obtaining  opti¬ 
mum  Wilson  parameters. 


2)  f  o(Po)  is  obtained  from  equation  (5.27),  where: 


s 

a)  ,  the  saturated  vapor  pressure  of  component  i,  is 

obtained  by  the  use  of  Antoine's  equation. 


g 

b)  <f>^  is  obtained  from  the  correlations  represented  by 
equations  (5.32)  to  (5.34)  inclusive. 


3)  P  is  the  total  pressure  of  the  system. 

4)  <Jk  is  calculated  by  equation  (5.25).  Simplifications  to 
that  equation  are  presented  below. 

5 . 5  The  Simplified  Prausnitz  Model  for  Vapor-Liquid  Equilibria 

Several  simplifications  and  approximations  can  be 
made  on  the  vapor-liquid  equilibrium  equation  presented 
above  without  affecting  its  accuracy  greatly.  This  is  par¬ 
ticularly  true  if  the  system  pressure  is  low^ . 

The  fully-developed  equilibrium  equation  has  been 
expressed  as  follows: 


RT 


(5.23) 


or 
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The  Fugacity  Coefficient  of  the  Vapor  Phase 

The  rigorous  expression  for  the  calculation  of  <J>^ 
is  equation  (5.25).  However,  as  Eckert‘S  points  out,  the 
determination  of  <}>^  requires  a  trial-and-error  procedure 
since  it  is  expressed  as  a  function  of  y. 

Van  Ness'*''*'  presents  an  excellent  treatment  of 
constant-composition  fluid  behavior  which  will  be  used  here 
to  justify  some  of  the  simplifications  we  can  bring  upon  ^ . 

At  constant  T,  $  is  rigorously  expressed  as 

1  P 

- - /  adP  (5.51) 

RT  o 

RT 

=  —  -  V  (5.52) 

P 

dP 

( Z  —  1 )  —  at  constant  T  (5.53) 

P 

(5.54) 

Z  can  be  written  as  a  power  series  in  P  explicit  in  V  as 
follows11 : 

BP  (C  -  B2)  9 

Z  =  1  +  —  +  - 2 -  P  +  .  .  .  (5.55) 

RT  (RT) 

where  B  and  C  are  second  and  third  virial  coefficients  res¬ 
pectively.  Since  the  virial  ecruation  truncated  after  the 
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second  virial  coefficient  is  a  highly  reliable  equation  of 
state  at  low  pressures'^  ,  then  equation  (5.55)  can  be  writ¬ 
ten  simply  as 

BP 

Z  =  1  +  —  (5.56) 

RT 

Substituting  (5.56)  into  (5.53) 

BP 

In  <p  =  —  (5.57) 

RT 

To  show  the  reliability  of  this  equation  relative  to  other 
more  rigorous  ones,  Van  Ness^  shows  plots  of  the  fugacity 
coefficient  of  SC>2  at  three  different  temperatures  vs.  pres¬ 
sure.  From  these  graphs  it  is  clear  that  below  a  pressure  of 
approximately  6.0  atms .  at  50°C,  equation  (5.57)  is  reliable 
for  this  system  and  at  250°C,  up  to  80  atmospheres  can  be 
tolerated  without  much  loss  in  accuracy.  Considering  that 
the  pressure  of  interest  is  normally  one  atmosphere  in  pro¬ 
cesses  such  as  extractive  distillation  and  that  SC>2  is  a 
highly  polar  gas,  it  was  decided  that  for  our  purposes,  equa¬ 
tion  (5.57)  would  be  used  to  compute  <j> .  The  correlation  de- 

1 2 

veloped  by  Pitzer  and  Curl  and  reported  by  Lewis  and  Ran- 
1  3 

dall  expresses  BP  /RT  as  a  function  of  w  and  T  .  The 
application  of  this  correlation  to  the  system  methyl  acetate- 
benzene-cyclohexane  yielded  values  such  that  1 . 0  2.  <f>  >_  0.950. 

The  general  closeness  of  the  value  of  <jj  to  1.0  has  been  con- 

Q 

firmed  by  O'Connell  et  al  .  Thus,  the  further  simplification 
to  =  1.00  was  made  without  seriously  affecting  the  equilibrium 
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predictions . 


The  Factor  <j>^s  exp  ((v^L/RT)  (P  -  P^s)) 


Several  authors-^,  in  calculating  liquid-phase 


activity  coefficients  from  equilibrium  data  use  the  follow¬ 
ing  equilibrium  equation: 


s 


(5.58) 


The  simplicity  of  this  equation  makes  it  very  attractive. 


exp  (  (viL/RT)  (P-PiS))/ 


•  s 

It  was  thus  decided  to  study  the  ratio  <f> . 


which  is  the  factor  applied  to  the  right-hand  side  of  the 
above  equation  (5.58)  to  give  the  right-hand  side  of  the 
rigorous  equilibrium  equation  (5.36).  For  the  system  methyl 
acetate-benzene-cyclohexane  at  one  atmosphere,  this  ratio 
ranged  from  0.9730  to  1.017  in  all  three  components.  It  was 
thus  concluded  that  for  processes  such  as  extractive  distil¬ 
lation  where  the  system  pressure  is  almost  always  atmospheric, 
assigning  a  value  of  unity  to  both  <|k  ,  the  fugacity  coeffici¬ 
ent  of  the  vapor  phase,  and  <j>^s  exp  ((v^L/RT)  (P-P^s)),  the  cor¬ 
rection  factor  in  the  liquid  phase  would  lead  to  a  negligible 
loss  in  the  accuracy  of  equilibrium  predictions. 

5 . 6  Summary 

Two  methods  of  calculating  equilibrium  K-values 


using  liquid-phase  activity  coefficients  obtained  from  the 
Wilson  equation  are  available. 
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The  simplified  method  can  be  described  in  steps  as 


follows : 

1)  The  equilibrium  equation  is 

y. 


k.  = 

i 


Y  .  P  .  ' 
'  1  1 


(5.59) 


x . 

l 


2)  The  activity  coefficients  y^  are  obtained  from  Wilson's 
parameters  A  ^  and  Ajj_*  To  obtain  these  parameters,  at 
least  one  totally  described  equilibrium  point  in  every 
binary  (x^,  y^,  y 2)  is  needed.  The  y  values  in  this 
equilibrium  point  should  have  either  been  obtained  from 
equilibrium  data  using  equation  (5.58),  or  this  equation 
should  be  used  to  calculate  them  from  the  equilibrium 
compositions . 

3)  P  is  the  total  pressure  of  the  system  and  P^s  is  the 
vapor  pressure  obtained  from  the  Antoine  equation. 

The  rigorous  method  can  be  described  in  steps  as 

follows : 


1) 


The  equilibrium  equation  is 

y. 


K.  =  —  = 

1 

l 


^ipiS*iS  exp  w 


(p  - 


pis> 


(5.60) 


P 


2) 


The  activity  coefficients  y^  are  obtained  from  Wilson's 
parameters  A^  and  A j ^ .  These  parameters  should  be 
obtained  from  equilibria  points  (at  least  one  for  every 
binary)  where  the  activity  coefficients  were  calculated 
by  the  rigorous  method,  i.e.  if  y^  were  obtained  by 
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equation  (5.58),  they  should  be  multiplied  by  the 

L 

i S  exp  (P  -  PiS) 

factor  -  .  These  liquid  activity 


coefficients  should  then  be  used  to  calculate  A. .  and 

ID 

A j ^ .  For  systems  at  atmospheric  pressure,  this  factor 
is  close  to  unity.  However,  in  a  rigorous  calculation, 
this  correction  in  y.  should  be  made. 

'  l 

3)  is  obtained  from  (5.32)  to  (5.34)  inclusively. 

4)  <f>^  is  obtained  from  equation  (5.57). 

Both  the  rigorous  and  the  simplified  methods  were 

used  to  predict  the  vapor  compositions  of  the  system  methyl 

acetate-benzene-cyclohexane.  The  results  of  the  study  in  the 

3  8 

form  of  standard  deviations  from  the  equilibrium  y^ ' s  is 
shown  in  Table  5.1,  and  shows  that  the  two  methods  are  com¬ 
parable  in  accuracy.  This  justifies  the  predominant  use  of 
the  simplified  method  in  the  design  of  extractive  distilla¬ 
tion  columns  in  this  work. 


Table  5.1 

Comparison  of  the  Rigorous  and  Simplified 
Prausnitz  Methods  of  Equilibrium  Predictions 


Pressure  =  1  atm. 


Standard  Deviation 

Rigorous  Simplified 

Model  Model 


methyl  acetate 

benzene 

cyclohexane 


0.03769 

0.00743 

0.00842 


0.04208 

0.00752 

0.00770 
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5 . 7  Calculation  of  Enthalpy  Data  with  Heats  of  Solution 

The  enthalpies  of  the  vapor  and  liquid  phases 
leaving  an  equilibrium  stage  can  be  calculated  from 


H  . 
3 


l  y-  •  h .  . 

i=l  13  13 


(5.61) 


h .  = 

3 


V  x..K.  . 
i=l  13  13 


(5.62) 


where  iL  ^  and  are  vapor  and  liquid  partial  molar  enthal¬ 

pies  of  component  i  in  the  solution.  Negligible  error  is 

involved  in  the  assumption  of  vapor-phase  ideality,  i.e.  that 
—  19 

,  and  the  enthalpy  of  the  vapor  phase  can  be  cal¬ 
culated  from  pure  component  enthalpies.  h^'s,  however,  can 
be  obtained  only  if  partial  molar  heats  of  solution  are 
available  as  a  function  of  composition.  The  component  en¬ 
thalpy  in  solution  can  then  be  expressed  as 


h  .  .  =  h  .  .  +  L  .  . 

13  13  13 


(5.63) 


For  most  hydrocarbons,  h^ .  and  .  can  be  obtained  from  "Data 

22 

Book  on  Hydrocarbons"  .  Enthalpies  for  numerous  other  com- 

20 

pounds  are  reported  in  the  Critical  Tables 

The  partial  molar  heats  of  solution  can  be  ob¬ 
tained  with  sufficient  accuracy  from  activity  coefficient  data 
by  using  the  equation 


d  log  y . 

( - -) 

d  i 

T 


L. 

1 


X' 


2 . 3R 


(5.64) 
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However ,  y^  vs.  T  at  constant  x-data  are  very  difficult  to 
obtain.  Thus,  two  methods  for  the  calculation  of  the  partial 

molar  heats  of  solution  L.  can  be  used. 

1 

1)  The  activity  coefficients  of  i  at  infinite  dilution  in  a 
reference  compound  r  iy^)  can  be  obtained  at  various 
temperatures  by  the  application  of  Pierrotti ' s  correla¬ 
tion^-?  Then,  using  equation  (5.64),  r ' s  can  be  ob¬ 
tained.  The  following  expression,  which  is  derived  from 
the  Redlich-Kister  equation,  is  then  used  to  obtain 
partial  molar  heats  of  solution  at  any  desired  composi¬ 
tion  : 


L  =  H (1-x  )  7  x . (L* .  +  La  . ) 

r  r  y  l  lr  ri 


-*5  T  (L  .  .  +L  .  . )  x  .  x 
.  , .  ii  i i  i 


i^r 


+  ^  l  (L  .  -L  . ) x . 
y  lr  n  l 


2x  (1-x  +x . ) -x . 
r  r  l  l 


I  (l\. 

.  k  .  ii 

^3  J 
i^r 


L  .  . ) x . x  .  (x . -x . ) 
l  j  l  j  l  j 


(5.65) 


where  i,j  and  r  refer  to  components.  This  equation  is 

19 

derived  by  Smith  in  full  detail. 

2)  Equation  (5.65)  can  be  used  without  recourse  to  Pier¬ 
rotti  ' s  correlation.  This  is  done  by  using  "The  Critical 
20 

Tables  " ,  which  are  a  source  of  a  large  number  of  expen 
mentally  determined  infinite  dilution  partial  molar  heats 
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VI.  EVALUATION  OF  PARAMETERS  FOR 


THE  WILSON  MODEL 


6 . 1  The  Model 

One  of  the  methods  used  in  this  work  to  obtain  li¬ 
quid  phase  activity  coefficients  is  the  Wilson  model.  It  has 
been  presented  above  as  follows: 


NCP  NCP  x.A., 

In  y,  =  -  In  (  I  x.A,  .)  +  1  -  T  - - 

k  jii  3  kJ  i=i  NCp 


y  x .  a  . . 

L -  3  i: 


j=l 


For  a  binary  solution,  this  reduces  to 


(5.50) 


A 


In  y 


ln<xl  +  A12x2)  +  x2( 


12 


xl  +  A12X2 


A 


21 


A21X1  +  x2 


(6.1) 


and 


In  y 


•In  (x 2  +  A2ixi^  “  xi  ( 


12 


X1  +  A12x2 


A 


21 


A21X1  +  X2 


(6.2) 


It  can  be  seen  from  these  equations  that  the  liquid  activity 
coefficients  of  the  components  in  a  solution  can  be  obtained 
as  a  function  of  the  composition  in  the  solution  and  the 
parameters  A^ j ,  called  Wilson's  parameters,  which  are  obtain- 
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able  from  the  binary  data.  is  defined  as 

L  A  .  .  -  X  .  . 

exp  -  (-U - ±±) 


A.  .  = 

ID 


JL 


(5.48) 


v . 

l 


RT 


where  A^  is  proportional  to  the  i-j  interaction,  etc.  Equa¬ 
tions  (6.1)  and  (6.2)  show  that  one  experimental  point  per 
binary,  (x^,  y^,  Y2)>rorP'  sufficient  to  yield  the  Wilson 
parameters  A^  an<^  A21  ^or  t^lat  system.  However,  it  is  found 
that  the  magnitude  of  both  the  Wilson  parameters  obtainable 
from  equations  (6.1)  and  (6.2)  and  the  interaction  energies 
obtainable  from  equation  (5.48)  depend  on  the  composition 
point  used  to  determine  them. 

Thus,  two  questions  arise: 

a)  What  are  the  optimum  Wilson  parameters  and/or  interaction 
energies?  i.e.  which  is  the  pair  of  parameters  that  will 
best  predict  the  binary  liquid  phase  activity  coefficients? 

b)  What  is  the  effect  of  using  non-optimum  Wilson  parameters 

or  interaction  energies  on  the  activity  coefficients 

\ 

predicted? 

An  attempt  at  answering  these  questions  is  made  in 


the  next  sections. 
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6 . 2  Determination  of  Optimum  Wilson  Parameters  and 

Interaction  Energies 


Mathematically,  the  equilibrium  equation  using 
Wilson's  model  can  be  presented  as  follows: 


(6.3) 


where  in  equation  (5. 59), for  example,  C  is  the  ratio  of  the 
vapor  pressure  to  the  total  pressure. 


Equation  (6.3)  can  also  be  presented  as 


by  combining  equation  (5.48)  with  (6.3). 


Two  methods  for  the  optimization  of  the  parameters 


in  these  equations  have  been  used  in  this  work: 
a)  The  first  method  is  a  finite  optimum  search  technique. 
Nineteen  points  in  every  binary  (covering  a  liquid  com¬ 
position  range  of  0.05  to  0.95  mole  percent  at  0.05  in¬ 
tervals)  are  each  used  to  yield  a  pair  of  Wilson  para¬ 
meters  and  interaction  energies.  Every  pair  of  these  para¬ 
meters  is  then  used  to  predict  a  range  of  vapor-phase  com¬ 
position.  To  make  this  possible,  the  activity  coefficients 
and  the  temperature  at  every  point,  and  the  total  pressure, 
have  to  be  known.  Each  y-prediction  is  then  compared  to 
the  experimental  values  by  determining  the  standard  devia¬ 
tion  o 
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where  n  is  the  number  of  points  in  the  prediction  range. 

The  optimum  pair  of  Wilson  parameters  (and  interaction 
energies)  is  then  the  one  which  yields  the  minimum  stan¬ 
dard  deviation  |d|.  Table  6.1  illustrates  the  range  in 
the  variations  of  the  parameters  with  composition  of  the 
binary  methylcyclohexane-toluene  as  an  illustration. 

Table  6.2  shows  the  values  of  the  optimum  Wilson  para¬ 
meters  and  optimum  interaction  energies  obtained  by  this 
technique  for  several  binaries.  In  predicting  multicom¬ 
ponent  equilibria  it  was  assumed  that  the  combination  of 
the  optimum  binary  parameters  would  yield  the  best  multi- 
component  prediction. 

b)  The  second  technique  uses  the  concept  of  quasilinearization 

to  predict  the  best  values  of  the  Wilson  parameters 

and  and/or  the  interaction  energies  (A^j  -  A^)  an<3 

(A..  -  A..).  The  data  needed  are  at  least  two  totally 
33 

described  composition  points.  An  initial  guess  of  the 
parameters  being  optimized,  henceforth  referred  to  as 

and  B^ ,  is  made,  a  range  of  y  is  predicted,  and  the 
prediction  is  compared  to  the  available  observed  points. 
Although  only  two  such  points  are  needed,  the  wider  the 
range  of  composition  that  the  available  information 
covers  the:  more  dependable  the  optimum  Wilson  parameters 
and/or  interaction  energies  will  be.  A  least  squares 
criterion  is  applied  to  obtain  a  new  set  of  A^  and  , 
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a  new  y-profile  is  predicted  and  the  process  is  continued 

until  the  difference  between  the  parameters  from  two 

-4 

successive  iterations  is  <JL0 

If  and  denote  the  initial  guess  points,  then 


y(1>  =  f<°>  +  (!!_,  <o>  (Aj<i>  .  A(o), 


3A 


ID 


ID 


ID 


+  (!i_)  t°>  (B<P  -  B<°>) 
3B  D1  jl 

ji 


(6.6) 


where  superscript  (1)  denotes  values  of  and  B_.^  to 

be  determined.  In  more  general  terms,  the  above  equation 
can  be  written  as 
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(n+1)  _  f(n)  +  ( _ j  (n)  ^A(n+1)  _  A (n)  ^ 
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If  a  deviation  is  dfined  as  E,  then 


E  =  y  .  -  <— )  (n)  A.(n+1) 

obs  3A. .  ^ 
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(6.8) 


The  known  terms  are  defined  as  parameters  in  the  equation 


E  =  y  -  PA . .  -  QB . .  -  R 
ID  D1 


(6.9) 
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where 


R  =  fv“'  -  ( 


(n)  _  ,3f  j(n)  A(n)  _  — j  <n)  B  (n)  (610) 

*3  31 


3A.  . 
ID 


3B  .  . 
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3  f 

Q  =  ( - ) 

3B  .  . 
31 


(n) 


(6.11) 


and 


3  f 

P  =  ( - ) 

3A.  . 
13 


(n) 


(6.12) 


squaring  the  error  for  absolute  magnitudes 


E2  =  (y  -  PA. .  -  QB . .  -  R) 2 

13  31 

.  .  .  .  2 

Minimization  of  E  corresponds  to  differentiating  the 

above  function  with  respect  to  A. .  and  B..  and  setting 

*  13  31 

the  derivatives  to  zero 


3E‘ 


=  -  2P* (y  -  PA. .  -  QB . .  -  R) 
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13 


31 


(6.13) 
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(y  -  PA. .  -  QB . .  -  R) 
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(6.14) 
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Since  these  functions  can  be  determined  for  as  many  pre¬ 
diction  points  as  observed  data  are  available,  then 


3E‘ 


3  A .  . 
13 


=  2Py  +  2A±^  J>2  +  26^  J>Q  +  2  ^PR  (6.15) 
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3E‘ 


3  B  .  . 


=  -l2Qy  +  2A.  J>Q  +  2B  £q2  +  2\  QR 


(6.16) 


li 

setting  the  error  functions  to  be  zero,  the  two  equations 

presented  can  be  solved  simultaneously  to  obtain  new 

values  of  A. .  and  B . . . 

ID  D1 

The  iterative  procedure  is  stopped  when  the  absolute  dif¬ 
ference  between  any  two  successive  values  of  A. .  and  B.. 

ID  D1 


-4 

is  less  than  a  certain  tolerance  (usually  10  ) . 

Table  6.3  shows  the  optimum  values  of  A. .  and  B.. 
for  some  binaries  obtained  by  this  method.  A  comparison  of 
the  figures  in  Tables  6.2  and  6.3  shows  that  the  two  tech¬ 
niques  yield  results  quite  different  in  magnitude.  However, 
as  will  be  seen  later,  the  effect  of  this  on  subsequent  pre¬ 
dictions  is  slight. 

Appendix  I  shows  the  expressions  for  the  differen¬ 
tials  and  algebraic  terms  introduced  in  equations  (6.6)  to 
(6.16)  inclusive.  Appendix  II  shows  the  computer  programs 
and  the  computer  outputs  used  in  the  optimization  studies. 


6 . 3  Sensitivity  of  Y-Predictions  to  Non-Optimal  Wilson 

Parameters  and  Interaction  Energies 

To  study  the  effect  of  using  non-optimal  para¬ 
meters  on  the  eventual  vapor-phase  composition  predictions, 
it  was  decided  to  use  equation  (5.59)  to  predict  isobaric 
ternary  acetone-methanol-water  vapor  compositions  at 
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1  atmosphere  using: 

a)  optimum  Wilson  parameters  obtained  by  technique  (a)  using 
isobaric  binary  data  at  1  atmosphere.  The  same  Wilson 
parameters  are  used  regardless  of  the  temperature  of  the 
ternary  point. 

b)  optimum  interaction  energies  obtained  by  technique  (a) 
using  isobaric  binary  data  at  1  atmosphere.  A  new  set 
of  A's  is  calculated  from  equations  (5.48)  and  (5.49)  at 
every  ternary  prediction  point. 

c)  optimum  Wilson  parameters  obtained  by  technique  (a)  using 
isothermal  binary  data  at  100°C.  Effect  of  temperature 
here  is  completely  neglected. 

d)  optimum  interaction  energies  obtained  by  technique  (a) 
using  isothermal  data  at  100°C.  A  new  set  of  A's  is 
calculated  from  equations  (5.48)  and  (5.49)  at  every 
new  ternary  temperature. 

e)  optimum  interaction  energies  obtained  by  technique  (b) 
using  isobaric  data  at  1  atmosphere.  A  new  set  of  A's 
is  calculated  at  every  prediction  point. 

The  values  for  these  parameters  are  reported  in 

Table  6.4.  In  every  one  of  the  above  cases  the  y-predictions 

34  35 

were  compared  to  the  literature  values  '  and  standard  de¬ 
viations  computed. 
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A  conclusion  that  can  be  reached  from  the  results 
in  Table  6.4  is  that  the  interaction  energies  have  no  tempera¬ 
ture  dependency,  i.e.,  once  the  optimum  values  for  a  binary 
have  been  calculated  they  can  be  used  with  confidence  to  pre¬ 
dict  vapor  equilibria  covering  a  wide  range  of  temperatures. 
Moreover,  interaction  energies  obtained  from  binaries  under 
some  conditions  can  be  used  to  predict  equilibria  at  reasonably 
removed  conditions,  e.g.,  the  optimum  interaction  energies 
obtained  from  the  binaries  at  100°C  can  be  used  without  a 
great  loss  in  accuracy  to  predict  vapor  phase  compositions  in 
the  isobaric  (1  atm.)  ternary  where  the  temperature  ranges 
from  58.2°C  to  70.0°C. 

The  Wilson  parameters  on  the  other  hand  are 
strongly  temperature  dependent.  Table  6.4  shows  that  regard¬ 
less  of  the  conditions  in  the  binaries,  their  accuracy  in 
prediction  is  poor  when  the  temperature  effect  is  not  taken 
into  consideration. 

The  two  optimization  techniaues  produce  inter¬ 
action  energies  and  Wilson  parameters  that  are  at  times  quite 
different.  However,  as  Table  6.4  shows,  these  different  op¬ 
timum  values  predict  the  ternary  data  equally  well. 
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VII.  STUDY  OF  MULTICOMPONENT  DISTILLATION  COLUMNS 
WITH  THE  GENERALIZED  COMPUTATIONAL  ROUTINE 

7 . 1  Detailed  Calculational  Procedure 

The  procedure  presented  in  Section  4.5  for  the  study 
of  distillation  columns  can  now  be  outlined  in  detail  as  fol¬ 
lows  : 

1)  Assume  initial  temperature,  vapor,  and  liquid  profiles. 

2)  Use  any  of  the  following  models  for  the  determination  of 
the  vaporization  constants  : 

a)  Polynomial  temperature  function  if  the  behavior  of  the 
components  in  the  system  can  be  considered  independent 
of  concentration 

b)  Chao-Seader  correlation  for  the  range  of  compounds 
for  which  it  is  applicable.  In  this  correlation, 

is  a  function  of  both  temperature  and  composition 

c)  The  Wilson  equation  for  the  determination  of  liquid- 
phase  activity  coefficients  coupled  with  the  use  of 
the  simplified  equilibrium  equation  (5.59). 

For  those  cases  such  as  (b)  and  (c)  where  the  K^'s  are  a 
function  of  both  temperature  and  composition  in  the  column, 
the  composition  profile  is  calculated  for  the  first  itera¬ 
tion  only  by  the  use  of  a  polynomial  or  an  ideal-K^ 
subroutine.  The  latter  subroutine  is  obtained  from  the 
Chao-Seader  model  by  assuming  the  liquid  phase  activity 
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coefficients  to  be  equal  to  the  vapor  phase  fugacity  co¬ 
efficients  . 

3)  Apply  the  tridiagonal  matrix  method  and  the  modified 
Gauss  elimination  algorithm  to  calculate  a  liquid  composi¬ 
tion  profile. 

4)  Calculate  a  new  temperature  profile  applying  Muller's 
quadratic  root  method. 

5)  Calculate  the  new  vapor  and  liquid  profiles  by  solving 
the  energy  balance  equations.  The  component  vapor  and 
liquid  enthalpies  can  be  obtained  by  any  of  the  following 
methods • 

a)  Polynomial  fits  of  the  enthalpies  with  temperature 
such  as  in  hydrocarbons  up  to  c^q* 

b)  Chao-Seader  enthalpy  correlation  for  the  hydrocarbons 
where  the  model  is  applicable. 

22 

c)  Component  enthalpies  from  the  literature  with  the 
necessary  corrections  to  take  the  mixing  effects  into 
consideration . 

6)  Return  to  step  (2)  using  the  newly  obtained  T ^ ' s  and  V ^ ' s 
and  repeat  the  procedure  until  the  sum  over  the  column  of 
the  squares  of  the  temperature  differences  between  two 
successive  iterations  is  less  than  a  certain  value  (usual¬ 
ly  0.01N  where  N  is  the  number  of  stages). 

The  computer  flow  diagram  for  both  composition- 
dependent  routines  is  shown  in  Figure  7.1. 
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Figure  7.1 

Flow  Diagram  :  S.R.  Ideal  K  and  S.R.  Wilson  (w) 
or  Chao-Seader  (C.S.) 
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The  first  part  of  this  chapter  presents  two  test 
problems.  These  problems  were  solved  in  the  literature  using 
the  same  method  outlined  above  and  the  results  were  available 
for  comparison. 

The  second  part  presents  a  stabilizer  problem  where 

.  .  33 

the  material  balance  equations  were  solved  in  the  literature 

by  the  Thiele-Geddes  method  and  in  this  work  by  the  tridiagonal 
matrix  method.  This  provides  an  idea  of  the  comparative  ef¬ 
ficiencies  of  the  two  methods.  The  detailed  statements,  the 
enthalpy  and  equilibrium  data  used,  the  initial  and  final 
profiles,  and  the  product  compositions  for  all  the  problems  in 
Chapter  VII  are  reported  in  Appendix  III. 

7 . 2  Problem  I:  Distillation  of  a  Hydrocarbon  System 

Using  a  One-Feed  Column 

26 

This  is  the  same  problem  reported  by  Wang  and  Henke 

23 

and  Amundson  and  Pontinen  .  The  one-feed  column  is  of  six¬ 
teen  equilibrium  stages,  with  stage  one  a  partial  condenser 
and  stage  sixteen  a  reboiler.  The  computational  routines 
used  with  this  problem  are  the  following: 

-  Tridiagonal  matrix  method 

-  as  a  polynomial  temperature  function 

-  Muller's  quadratic  root  method 

-  Enthalpies  as  a  polynomial  temperature  function 

-  Conventional  energy  balances 
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The  results  from  this  work  duplicated  those  of 
Wang  and  Henke^  in  four  iterations  as  compared  to  the  reported 
seven  This  improvement  is  due  to  the  normalization  of  the 
liquid  composition.  Total  ob'ject  time  for  this  problem  was 
fourteen  seconds  on  an  IBM  7040  computer.  This  is  equivalent 
to  0.04  sec/component/stage/iteration. 

7 . 3  Problem  II:  Distillation  of  a  Hydrocarbon  System 

Using  a  Complex  Column 

This  problem  was  also  studied  by  Wang  and  Henke^0 
2 c 

and  Lyster  ~.  The  system  to  be  separated  contains  eleven 
hydrocarbon  components.  The  column  uses  two  feeds  and  two 
side  streams,  with  a  total  of  21  equilibrium  stages. 

The  routines  used  with  this  problem  are  the  fol¬ 
lowing  r 

-  Tridiagonal  matrix  method 

-  as  a  polynomial  function  of  temperature 

-  Muller's  quadratic  root  method  for  the  determina¬ 
tion  of  the  bubble  point  temperature 

-  Enthalpies  as  a  polynomial  temperature  function 

-  Conventional  energy  balances 

The  program  converged  in  10  iterations,  the  computation  time 
being  0.036  secs/component/stage/iteration  and  the  total  ob¬ 
ject  time  being  1  min.  and  23  seconds. 
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7 . 4  Problem  III:  Use  of  the  Chao-Seader  Correlation 

with  Problem  II 

Since  one  of  the  equilibrium  and  enthalpy  models 
specifically  applicable  to  many  hydrocarbon  systems  is  the 
Chao-Seader  correlation,  it  was  decided  to  use  the  Chao-Seader 
correlation  for  equilibrium  and  enthalpy  data  with  problem  II 
using  the  same  computational  procedure  as  outlined  above. 

The  program  converged  in  7  iterations  with  a  computation  time 
of  0.08  secs/component/stage/iteration  and  a  total  object 
time  of  2  minutes  and  14  seconds. 

Since  the  only  difference  between  problems  II  and 
III  is  the  method  of  obtaining  the  equilibrium  and  enthalpy 
data,  it  can  be  seen  from  the  computation  times  that  the 
Chao-Seader  correlation  takes  almost  3  times  as  long  as  the 
polynomial  model  in  yielding  equilibrium  and  enthalpy  data. 

The  number  of  iterations  in  this  problem  was  7 .  This  pro¬ 
gram  shows,  however,  that  the  Chao-Seader  correlation  can  be 
successfully  used  with  a  complex  column  (more  than  one  feed 
and  two  product  streams) . 

Appendix  V  presents  the  detailed  computer  program 
for  problem  III  as  an  example  of  the  programming  used  with 
the  problems  in  this  chapter. 
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7 . 5  Problem  IV:  Use  of  Tridiagonal  Matrix  and  Muller  Methods 

with  a  Stabilizer  Column 

39 

Tatuch  solved  this  stabilizer  problem  using  the 
Thiele-Geddes  method  to  obtain  composition  profiles.  The 
same  problem  was  solved  in  this  work  using  the  tridiagonal 
matrix  method. 

The  routines  used  with  this  problem  were  as  follows: 

-  Tridiagonal  matrix  method 

-  Chao-Seader  K.  .  '  s 

ID 

-  Muller's  quadratic  method 

-  Chao-Seader  enthalpies 

-  Conventional  enthalpy  balances 

The  program  converged  in  10  iterations  with  a  computation  time 

of  0.12  secs/iteration/component/stage.  The  total  object 

time  was  5  minutes  and  27  seconds  compared  to  3  minutes  and 

3  9 

55  seconds  reported  by  Tatuch  .  However,  these  object  times 

39 

are  not  strictly  comparable  since  Tatuch  does  not  report 
the  initially  assumed  profiles. 
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VIII.  EXTRACTIVE  DISTILLATION:  PAPT  I 

The  work  done  on  the  distillation  of  hydrocarbon 
systems  presented  above  showed  that  the  tridiagonal  matrix 
method  is  adequate  for  the  simultaneous  solution  of  the 
material  and  energy  balance  equations  in  complex  columns  (more 
than  one  feed  and  more  than  two  product  streams) .  This  work 
also  showed  that  the  quadratic  fit  algorithm  for  the  determina¬ 
tion  of  the  bubble  point  vector  across  a  column  (the  Muller 
method)  is  adequate.  The  Wilson  equation  in  Section  6.3  was 
shown  to  be  capable  of  predicting  multicomponent  vapor-liquid 
equilibria  using  data  obtained  from  the  binary  systems  only. 

It  was  thus  decided  to  study  some  extractive  distillation 
columns  using  these  computational  routines.  All  extractive 
distillation  columns  are  two-feed  columns.  The  systems 
handled  in  these  columns  always  show  large  deviations  from 
Raoult's  law  so  that  the  application  of  these  routines  to 
extractive  distillation  would  be  a  good  test  for  their  ver¬ 
satility.  The  flow  diagram  for  the  computational  procedure 
has  been  presented  in  Section  7.1  and  Figure  7.1. 

8 . 1  Problem  V:  Extractive  Distillation  of  the  System 
Methylcyclohexane  (MCH) -Toluene  with  Phenol 

The  problem  is  to  optimally  separate  100  moles/unit 
time  of  a  50-50  MCH-toluene  mixture.  The  binary  separation 
of  MCH  and  toluene  is  extremely  difficult  since  MCH  boils  at 
101°C  and  toluene  boils  at  110. 8°C. 
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Accordingly  phenol  is  used  as  the  extractive 
agent  since  it  satisfies  all  the  properties  specified  for  a 
solvent  in  Chapter  III. 

19 

Problem  V  was  studied  by  Smith  where  he  presented 
the  first  iteration  in  calculating  the  number  of  stages  re¬ 
quired  to  achieve  a  specified  separation.  Since  the  physical 
configuration  of  the  column  is  necessary  prior  to  using  the 
programs  developed  in  this  work,  a  preliminary  study  was  neces¬ 
sary  to  determine  the  stages  at  which  the  solvent  and  fresh 
feed  streams  should  be  introduced. 

The  number  of  specifications  required  for  this  column 
is  2C  +  2N  +  12.  The  choice  made  from  the  2C  +  2N  +  15  specifi¬ 
cations  reported  in  Section  3.6  is  the  following: 


Pressure  in  N-stages  N 

Pressure  in  condenser  1 

Pressure  and  heat  leak  in  reflux  divider  2 

Heat  leak  in  N-l  stages  N  -  1 

Condenser  load  1 

Total  number  of  stages  1 

Fresh  feed  composition  temperature  and  rate  C  +2 

Fresh  feed  location  1 

Solvent  feed  composition  temperature  and  rate  C  +2 

Solvent  feed  location  1 

Rate  of  overhead  product  1 

L/D  ratio  1 


2C  +  2N  +  12 

Figure  8.1  shows  the  result  of  the  study  relating  to 


the  feed  introduction  stages. 
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Solvent  Feed  Plate  Location 

Figure  8.1  Study  to  Locate  Solvent  and  Fresh  Feed  Plates 

Total  Plates  =  21 
Paramete  f  =  l/D 
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This  preliminary  study  shows  that  for  a  21  stage 
column,  the  best  separation  is  obtained  by  introducing  the 
solvent  feed  in  staae  5  and  the  fresh  feed  in  stage  14. 

The  computational  routines  used  with  this  problem 
were  as  follows: 

-  Tridiagonal  matrix  method 

-  from  the  simplified  Prausnitz  model  presented 
in  Section  5.5.  The  liquid-phase  activity  co¬ 
efficients  are  obtained  from  the  Wilson  model 

-  Muller's  quadratic  root  method 

-  Enthalpies  from  pure  component  enthalpies  and 
heats  of  mixing  corrections. 

The  optimized  values  of  the  interaction  energies  in 
the  Wilson  equation  and  the  enthalpy  data  are  reported  in 
Appendix  IV  along  with  the  problem  statement,  the  initial 
and  final  profiles  and  the  product  compositions. 

8 . 2  Enthalpy  Data  in  Extractive  Distillation 

The  enthalpy  data  for  the  systems  separated  by  ex¬ 
tractive  distillation  were  computed  as  follows: 

1)  The  binary  infinite-heats  of  dilution  were  obtained  from 
19 

Smith  for  system  phenol-MCH-toluene  and  from  the 

2  0 

"Critical  Tables"  for  system  acetone-methanol-water. 
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2)  Equation  (5.65)  was  used  to  obtain  partial  molar  heats 
of  solution  at  the  compositions  of  the  stages 

3)  Equation  (5.63)  was  used  to  calculate  enthalpies  of  the 
liquid  phase.  Ideality  of  the  vapor  phase  was  assumed. 

The  equations  representing  the  variation  of  the  enthal¬ 
pies  of  the  ideal  vapor  and  liquid  phases  with  temperature 

22 

were  derived  from  the  Data  Book  on  Hydrocarbons 

8 . 3  Optimization  of  the  Column 

With  this  information  available,  the  optimization 
of  the  operation  variables  of  the  column  such  as  reflux  ratio 
and  solvent  and  fresh  feed  rates  could  now  be  studied.  The 
separation  efficiency  was  characterized  by  one  of  the  follow¬ 
ing  two  parameters : 

-  Recovery  of  toluene 

-  Ratio  of  MCH/toluene  in  the  bottom  product 
Figure  8.2  shows  the  dependence  of  MCH/toluene 

ratio  in  the  bottoms  on  rate  of  solvent  feed  and  L^/D  ratio. 
With  increasing  L^/D  ratio  and  increasing  solvent  rate,  the 
purity  of  the  bottom  product  (the  toluene  product)  improves. 

The  dependence  of  the  recovery  of  toluene  on  the 
same  two  variables  is  shown  in  Figure  8.3.  We  can  note  from 
this  figure  that  for  a  fixed  L^/D  ratio,  a  higher  recovery  is 
sometimes  possible  with  a  lower  solvent  feed  rate.  To  show 
this  effect  clearly.  Figure  8.4  shows  the  variation  of  solvent 
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Figure  8.2  Variation  of  MCH/Tol  Ratio  with  the  Reflux  Ratio 
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feed  rate  with  the  L^/D  ratio  for  toluene  recoveries  of 
76,  77  and  78%.  This  figure  shows  that  there  is  a  particular 
reflux  rate  that  requires  the  least  amount  of  solvent.  For 
a  recovery  of  77%,  an  L^/D  of  5.3  requires  the  least  solvent 
per  fresh  feed  rate  ratio.  For  a  different  length  column,  the 
shape  of  the  curve  will  be  the  same  though  its  location  will 
be  different10. 

Table  8.1  shows  the  numerical  data  obtained  from 
the  optimization  study  of  the  system  phenol-MCH-toluene . 

The  detailed  computer  program  used  with  this  system  is  identi¬ 
cal  to  the  one  used  with  the  system  acetone-methanol-water 
and  will  be  presented  in  the  next  chapter. 
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Table  8 . 1 

Numerical  Results  for  the  Optimization 
of  the  System  Phenol-MCH-Toluene 
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IX.  EXTRACTIVE  DISTILLATION:  PART  II 

9 . 1  Problem  VI:  Extractive  Distillation  of  Acetone-Methanol 

The  possibility  of  using  the  computational  routines 
developed  in  this  work  to  study  and  optimize  extractive  dis¬ 
tillation  was  shown  in  Chapter  VIII  on  the  phenol-methylcyclo- 
hexane-toluene  system.  Since,  however,  it  was  not  possible 
to  rigorously  verify  the  results  obtained,  it  was  decided  to 

a)  simulate 

b)  optimize 

an  existing  industrial  column. 

Chemcell  (1963)  Ltd.  provided  us  with  all  the  re¬ 
quired  information  for  the  study  of  their  acetone-methanol 
extractive  distillation  tower  with  water  as  the  solvent. 

This  information  is  presented  in  Figure  9.1  and  Table  9.1 
presented  below. 

9 . 2  Simulation  of  the  Industrial  Column 

In  the  simulation  of  an  existing  industrial  separa¬ 
tion  column,  the  first  attempt  is  to  measure  the  operating 
conditions  and  performance.  This  will  yield  component  balances 
around  the  column  when  coupled  with  accurate  analyses  on 
the  ingoing  and  outgoing  streams.  The  next  step  is  to  vary 
the  number  of  stages  in  every  section  in  the  column  until  the 
overhead  and  bottoms  compositions  are  duplicated. 
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Figure  9.1 


Diagram  of  Industrial  Acetone-Methanol  Tower 
(Chemcell  Ltd.) 
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However,  obtaining  a  perfect  material  balance  around 
the  column  for  every  component  is  almost  impossible  with  the 
chemical  system  under  study.  This  is  due  to  the  fact  that 
laboratory  analyses  for  acetone-methanol  mixtures  are  rather 
difficult.  As  a  result,  the  two  feed  rates  and  their  composi¬ 
tions  and  the  overhead  rate  and  composition  were  assumed  to 
be  the  correct  ones.  A  material  balance  around  the  column 
with  these  values  showed  that  the  reported  bottoms  flow  rate 
was  3.15%  higher  than  the  calculated  one. 

The  flow  diagram  shown  in  Figure  7.1  is  the  one 
used  with  this  problem  also.  The  specified  degrees  of  free¬ 
dom  are  the  same  as  those  specified  for  the  optimization  of 
the  phenol-MCH-toluene  system. 

The  variables  whose  influence  on  the  simulation  was 
studied  were  the  total  number  of  equilibrium  stages  and  the 
locations  of  the  feeds.  The  other  specifications  were  held 
constant  at  the  same  values  as  reported  in  Figure  9.1  and 
Table  9.1. 

Table  9.2  shows  the  results  from  the  simulation 
and  compares  them  to  the  reported  information  on  the  indus¬ 
trial  column.  The  problem  statement,  initial  and  final  pro¬ 
files  and  product  compositions  are  reported  in  Appendix  IV. 

The  detailed  computer  program  is  reported  in  Appendix  V. 
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9 . 3  Activity  Coefficient  Predictions  with  the  Wilson  Model 

at  Low  Concentrations 


As  was  mentioned  in  Chapter  VII,  the  convergence 
criterion  ordinarily  set  for  distillation  problems  is  that  the 
sum  over  the  column  of  the  squares  of  the  temperature  differ¬ 
ences  between  two  successive  iterations  be  less  than  0.01  times 
N,  where  N  is  the  number  of  stages.  In  this  thesis,  the  con- 


2 

vergence  criterion  set  for  £  ( AT  )  was  1.0  or  less.  This 

j  =  l 

criterion  was  generally  met  with  the  phenol-MCH-toluene  system 
studied  in  the  preceding  chapter.  With  the  acetone-methanol- 
water  system,  however,  it  was  not  possible  at  first  to  bring 


N  9 
l  (AT2) 
j=l  N 
When  £ 

j=l 


below  a  certain  minimum  value,  often  close  to  10.0. 
2 

(AT  )  reached  that  point,  it  would  increase  again 


for  a  few  iterations,  and  oscillations  resulted.  This  trend 
to  oscillate  occurred  in  all  cases  except  when  the  solvent 
feed  rate  and  the  reflux  ratio  were  very  low.  In  examining 
this  difficulty  several  parameters  that  were  thought  relevant , 


r  2 

including  l  (AT  )  were  plotted  vs.  the  iteration  m  the  man- 

j=l 

ner  of  Figure  9.2.  From  this  graph  it  was  noticed  that  the 

N  2 

trend  in  the  variation  of  7  (AT  )  was  closely  followed  by  a 

3=1 

similar  trend  in  the  bottom  compositioh  of  acetone,  component 
1,  and  was  not  followed  by  any  of  the  other  plotted  variables. 
It  was  thus  hypothesized  that  as  the  concentration  of  acetone 
in  the  bottoms  decreased,  which  would  be  the  expected  trend 
toward  convergence  acetone  being  the  top  product,  the  ability 
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of  the  Wilson  model  to  predict  its  activity  coefficient  was 
severely  curtailed.  This  would  consequently  lead  to  poor 
predictions  of  the  vapor  phase  compositions,  and  the  accuracy 
of  the  bubble-point  temperatures  would  thus  be  affected. 

This  trend  would  probably  continue  until  the  increasing  ace¬ 
tone  concentrations  reached  a  level  where  better  y  predictions 
by  the  Wilson  equation  were  possible. 

To  test  this  hypothesis,  a  program  was  written  to 
test  the  thermodynamic  consistency  between  each  two  succes¬ 
sive  points,  where  the  composition  of  both  acetone  and  methanol 
-n  the  system  varied  from  0.0005  to  0.5120  The  Wilson  model, 

with  the  optimum  binary  Wilson  parameters,  was  used  to  predict 
the  activity  coefficients  of  every  component. 

The  consistency  test  was  one  proposed  by  McDermott 
40 

and  Ellis  and  can  be  written  as 


D  = 


(xij  +  xi,j+i)(log  Yi,j+i 


log  y  .  .  ) 

i/D 


where  C  is  the  number  of  components,  j  and  j+1  two  successive 
points  and  D  is  therefore  directly  related  to  the  errors  in 
activity  coefficients  of  the  points  concerned.  Tables  9.3  and 
9,4  show  some  of  the  trends  revealed  by  the  study. 

Table  9.3  shows  (at  constant  concentrations  of  ace¬ 
tone  in  the  ternary)  the  concentration  of  methanol  and  the  in¬ 
consistency  parameter  D.  A  direct  relationship  between  methanol 
concentration  and  the  inconsistency  is  clear.  Table  9.4  shows 
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(at  constant  methanol  concentration)  the  concentration  of 
acetone  and  the  inconsistency  parameter,  and  the  trend  here 
is  a  decrease  in  the  inconsistency  with  increasing  acetone 
concentration.  These  results  seemed  to  prove  that  our  hypo¬ 
thesis  about  the  mechanism  that  triggered  the  oscillations  in 
the  convergence  criterion  was  correct.  We  could  add  to  that, 
however,  that  is  is  not  the  low  concentrations  of  acetone  in 
the  system  'per  se '  that  is  the  cause  of  erroneous  predictions, 
but  rather  the  low  concentrations  of  acetone  when  accompanied 
by  relatively  large  concentrations  of  methanol.  To  avoid  fur¬ 
ther  lack  of  convergence,  the  concentration  of  acetone  in 
the  last  two  stages  (reboiler  and  one  above)  was  set  equal  to 
zero.  In  no  case  previously  had  any  of  these  concentrations 
been  larger  than  5%,  so  that  the  errors  involved  in  this  ap¬ 
proximation  should  be  acceptable. 

9 . 4  Optimization  of  the  Acetone-Methanol-Water  Tower 

The  effect  of  variations  in  the  solvent  feed  rate 
and  the  reflux  ratio  on  the  performance  of  the  column  was 
studied  with  this  system  also.  The  column  used  contained 
25  ideal  stages,  with  the  solvent  feed  being  introduced  in 
stage  8  and  the  process  feed  in  stage  20.  This  configuration 
was  chosen  for  the  closeness  of  the  acetone  concentration  in 
D  to  that  of  the  actual  column.  Figure  9.3  shows  the  depen¬ 
dence  of  percent  methanol  recovery  on  the  reflux  ratio  at 
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different  solvent  feed  rates.  The  recovery  of  methanol 
figures  are  only  relative  due  to  the  setting  of  the  acetone 
concentration  in  B  to  0.000. 

It  is  clear  from  this  curve  that  a  reflux  ratio  of 
5.0  is  optimum  for  all  solvent  feed  rates  studied.  Chemcell 
uses  an  L^/D  ratio  of  4.13,  though  it  is  understood  that  a 
certain  amount  of  internal  reflux  takes  place. 

Figure  9.4  shows  the  dependence  of  the  solvent  feed 
rate  on  L/D  ratio  at  the  particular  methanol  recovery  of  97%. 
Again  it  is  shown  that  there  is  a  particular  reflux  rate  that 
requires  the  least  amount  of  solvent. 

The  detailed  results  of  the  optimization  work  is 
presented  in  Table  9.5.  Comparing  these  results  to  Table  9.1 
shows  that  the  use  of  50.0  moles/hr.  of  water  (instead  of 
115.56)  at  an  L^/D  ratio  of  5.0  should  yield  an  acetone  con¬ 
centration  in  D  of  0.9522,  compared  to  the  industrial  values 
of  0.957. 
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Table  9 . 5 


Optimization  of  System  Acetone-Methanol-Water 


Total 

no.  of  plates 

=  25 

Solvent  feed  location 

=  8 

Fresh 

feed  location 

=  20 

Fresh 

feed  rate 

=  492.87 

Water 

concentration  in  D 

0.100 

moles/U . T . 


Solvent 

Feed 

moles/U . T . 


Relative 
,  %  Recovery 

r  of  Methanol 


Water 

Concentration 
in  D 


Acetone 
Concentration 
in  D 


25.0 

1 

87.99 

0.0901 

0.7040 

3 

95.64 

0.0206 

0.9021 

5 

97.27 

0.0210 

0.9264 

7 

96.68 

0.0174 

0.9208 

9 

95.88 

0.0143 

0.9115 

50.0 

3 

96.32 

0.0235 

0.9100 

5 

99.43 

0.0285 

0.9522 

7 

99.12 

0.0254 

0.9506 

9 

98.59 

0.0226 

0.9453 

100.0 

3 

96.09 

0.0243 

0.9071 

5 

98.98 

0.0282 

0.9479 

7 

99.49 

0.0279 

0.9559 

9 

100.0 

0.0286 

0.9645 
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X.  DISCUSSION  OF  RESULTS 

10.1  Equilibrium  Data  Sources 

Several  correlations  providing  vapor-liquid  equili¬ 
brium  data  have  been  used  during  this  work.  Problems  I  and 
II  used  a  polynomial  correlation.  The  main  goal  with  those 
problems  was  to  check  the  soundness  of  the  developed  programs 

by  comparing  the  solutions  to  published  information.  Both 

26 

problems  duplicated  the  results  of  Wang  and  Henke  ,  problem 
I  converging  in  4  rather  than  the  reported  7  iterations  due  to 
our  normalization  of  the  liquid  composition  at  the  end  of 
every  iteration. 

The  only  difference  between  problems  II  and  III  is 
that  while  problem  II  uses  polynomial  data,  problem  III  studies 

3 

the  same  problem  using  Chao-Seader  equilibrium  data.  It  can 
be  seen  from  the  results  in  Appendix  III  that  temperature  dif¬ 
ferences  higher  than  6°F  on  some  stages  can  be  caused  by  the 
use  of  the  different  data  sources. 

A  relatively  new  equilibrium  correlation  is  that 

4 

proposed  by  Orye  and  Prausnitz  .  In  this  model,  the  liquid 

phase  activity  coefficients  are  calculated  by  the  Wilson 
14 

equation  .  This  model  was  used  in  problems  V  and  VI.  The 
systems  treated  in  those  problems  were  common  enough  that  no 
difficulty  was  encountered  in  finding  binary  equilibrium  data 
for  them.  At  least  one  point  in  every  binary  is  necessary  in 
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the  determination  of  the  Wilson  or  the  interaction  energy 
parameters  which  are  then  used  in  the  prediction  of  the  liquid 
phase  activity  coefficients. 

10.2  Discussion  of  the  Wilson  Model 

Although  only  one  equilibrium  point  in  every  binary 
is  necessary  for  the  determination  of  the  interaction  energy 
parameters,  it  can  be  seen  from  Table  6.1  that  the  values  of 
the  interaction  energy  parameters  and  consequently  their 
ability  to  predict  activity  coefficients  vary  with  the  compo¬ 
sition  point  at  which  they  are  determined.  This  necessitated 
their  optimization  to  determine  the  pair  of  interaction  ener- 

4 

gies  which,  when  used  with  the  Prausnitz  model  ,  predicted 

the  equilibrium  data  with  the  minimum  standard  deviation. 

14 

The  Wilson  model  for  the  calculation  of  liquid 
phase  activity  coefficients  offers  two  outstanding  advantages 
over  other  methods : 

1)  The  liquid  phase  activity  coefficients  in  a  multicomponent 
system  can  be  calculated  from  parameters  obtained  from  the 
binary  equilibria  in  the  system.  This  is  an  important 
advantage  since  binary  equilibrium  data  are  more  readily 
available  than  multicomponent  data. 

2)  As  has  been  shown  in  Chapter  VI,  the  interaction  energy 
parameters  are  independent  of  temperature,  i.e.,  once 
the  optimum  interaction  energies  have  been  calculated  for 
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a  system  under  particular  conditions,  they  can  be  used 
to  calculate  liquid  activity  coefficients  for  the  same 
system  under  removed  temperature  conditions.  This  is  a 
strong  point  in  favor  of  the  Wilson  model.  It  makes  the 
model  particularly  useful  in  the  design  of  isobaric  dis¬ 
tillation  columns  where  large  temperature  variations  occur 
in  the  column. 

The  relative  independence  of  the  interaction  ener¬ 
gies  from  temperature  has  been  well  demonstrated  in  Table  6.4 
where  optimum  interaction  energies  obtained  from  isobaric 
binaries  at  100°C  predicted,  using  the  Prausnitz  model,  the 
vapor-liquid  equilibria  of  the  same  binaries  at  temperatures 
ranging  from  58.2  to  70.0°C. 

10.3  Relative  Timings  of  the  Equilibrium  Models 

Table  10.1  shows  the  relative  timing  for  the 
equilibrium  ratio  correlations  used  in  this  work.  From  this 
table,  it  is  evident  that  the  time  required  increases  with 
the  complexity  of  the  correlation  used.  The  table  shows 
that  it  takes  roughly  three  times  as  long  to  use  the  Chao- 
Seader  correlation  as  it  does  to  use  polynomial  equilibrium 
ratios.  The  Prausnitz  correlation  using  the  Wilson  model 
takes  roughly  four  times  as  long  as  the  polynomial. 
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10.4  Comparison  of  Enthalpy  Data  Sources 

Other  than  the  polynomial  fit,  the  two  main  sources 

3 

of  enthalpy  data  in  this  work  have  been  the  Chao-Seader  model 
and  ideal  enthalpies  with  adjustments  made  for  heats  of  solu¬ 
tion.  The  advantage  of  the  Chao-Seader  model  over  other  semi- 
empirical  data  sources  is  the  fact  that  the  parameters  calculated 
in  the  process  of  determining  equilibrium  ratios  can  also  be 
used  to  calculate  enthalpy  departures  for  the  same  components. 

The  second  method  of  enthalpy  data  calculations 
has  been  used  with  the  two  extractive  distillation  problems 
V  and  VI.  The  two  methods  for  the  calculation  of  heats  of 
solution  were  presented  in  Section  5.7.  The  second  method, 
which  avoids  the  use  of  the  Pierotti  correlation,  was  used 
exclusively  in  this  work.  As  will  be  shown  later,  the  success 
of  the  simulation  work  done  with  the  acetone-methanol-water 
system  in  problem  VI  shows  that  the  enthalpy  data  used  were 
adequate . 

10.5  Method  of  Solution  and  Convergence 

The  method  exclusively  used  in  this  work  to  solve 

the  material  balance  equations  and  determine  the  component 

distribution  in  the  column  was  the  tridiagonal  matrix  method 

2  6 

coupled  with  a  Gauss  elimination  algorithm  .  The  method  is 
numerically  stable  and  the  subtraction  of  nearly  equal  quanti- 
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ties  is  avoided.  Also,  no  matrix  inversion  operations  are 
necessary  which  was  the  disadvantage  of  the  matrix  method 
presented  by  Amundson  .  Muller's  method  for  the  determina¬ 
tion  of  bubble  point  temperatures  was  also  exclusive  to  all 
other  methods.  In  most  cases  convergence  occurred  in  the 

second  iteration.  The  use  of  the  Muller  method  has  been  com- 

28 

puted  by  others  to  cut  solution  times  by  at  least  a  factor 

of  three  over  the  Erbar  technique. 

Problem  IV  used  both  of  these  algorithms  to  study 

a  stabilizer  problem  which  has  already  been  studied  using 

39 

the  Thiele-Geddes  method  for  the  solution  of  the  material 
balance  equations.  The  computation  time  for  the  tridiagonal 
and  Muller  algorithms  was  5  minutes  and  29  seconds,  as  com¬ 
pared  to  3  minutes  and  55  seconds  for  the  Thiele-Geddes 

method.  These  timings,  however,  are  not  strictly  comparable 

•  3  9 

since  the  initial  assumptions  used  by  Tatuch  were  not  re¬ 
ported.  Our  initial  assumptions  and  results  are  reported  in 
Appendix  III. 

The  distinct  advantage  of  the  tridiagonal  matrix 
solution  is  the  ease  with  which  it  can  handle  additional 
feeds  and  sidestreams.  The  solution  model  makes  provisions 
for  a  feed  stream,  a  liquid  and  vapor  sidestreams  and  a  heat 
exchange  device  in  every  stage. 
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10.6  Possibility  of  Column  Simulation  and  Optimization 

The  results  of  the  simulation  of  an  industrial  column 
were  presented  in  Chapter  IX.  The  system  was  an  acetone- 
methanol  mixture  being  extractively  distilled  with  water  and 
the  top  product  was  acetone.  The  information  supplied  by 
Chemcell  Ltd.  was  reported  in  Figure  9.1  and  Table  9.1.  This 
information  can  be  summarized  as  follows: 

-  Composition  of  feed  and  product  streams. 

-  Flow  rates  and  temperatures  of  the  feed  and 
product  streams. 

-  Temperatures  at  a  number  of  points  in  the  column. 
This  information  allowed  the  assumption  of  initial  vapor, 
liquid,  and  temperature  profiles.  The  following  step  in  the 
simulation  was  to  duplicate  the  product  compositions  by  vary¬ 
ing  the  number  of  stages  and/or  the  feed  plate  locations. 

The  same  reflux  ratio  and  feed  rates  as  in  the  industrial 
column  were  used. 

Once  the  product  compositions  were  duplicated,  the 
optimization  of  the  reflux  ratio  and  the  ratio  of  the  feeds 
required  the  holding  of  the  total  number  of  plates  and  the 
feed  locations  constant.  This  is  tantamount  to  assuming  con¬ 
stant  plate  efficiencies  regardless  of  the  magnitudes  of  the 
operating  variables,  which  limits  the  value  of  the  optimiza¬ 
tion.  Accordingly,  if  the  optimum  values  of  the  operating 
variables  are  widely  different  from  their  values  in  the 


laiidaubni  ns  do  noxdsXumia  arid  do  adXi/ga*  9fi? 


bns  39iaw  rfdiw  belllSkib  ylavidosxdxa  prided  aiudxim  Xonsrfdam 

.1.8  aldsT  bns  1.8  sxupi'S  ni:  banoaBi  asw  .  bdJ  IXaomariD 


jawollod  as  bssilisminue  ad  nso  noxdsnraodni 


amsaida  douboiq  bns  baad  do  noidxaoqmoD 


bns  bead  arid  do  eainds-xaqmad  bns  adds*  woll  - 

. bijibbiSb  douboxq 

n  5  dS  8'9St- 


»^oqsv  IsidinX  do  noiiqflinaEs  arid  bawoXIs  noidsimodni?  siriT 

■ 


-y isv  yd  snoidXsoqmoo  doubosq  arid  adsoxXqi/b  od  asw  noidslr/nua 


>X  sdslq  bead  arid  ?o\bns  aepsda  do  isdmun  9rfd  pni 
Isx-xdenbni  arid  ni  as  sads?  bead  bns  oxdsx  xnlde?  amse  ariT 

" 


do  oidsi  arid  bn&  oil  a?  xnXdex  arid  do  n  c  * 
arid  bns  eadslq  do  xednu/n  Xsdod  arid  do  pnibXori  arid  basinpai 

arid  do  89budinpsm  arid  do  aaaXbxspad  aaxonaioidda  adsXq  dnsda 
±jSiraxdqo  9rid  do  anisv  9rid  adxmiX  rfoxriw  ^eaXdsixsv  pnidsxeqo 


pnidsxaqo  9fid  do  aauXsv  mufflidqo  arid  di  , yXpnibioooA  .noid 


arid  ni  aenlsv  ixarid  mold  dnsiaddib  yXabiw  bib  BBldBiTLBV 


102 


simulated  column,  a  new  simulation  should  be  attempted. 

A  comparison  of  some  of  the  simulated  and 
industrial  temperatures  is  shown  below: 

Simulation  Industrial 

Condenser  input  (V2)  133.24°F  135. 5°F 

Bottoms  product  194. 3°F  approx.  194°F 

The  closeness  of  the  simulated  and  the  industrial 
temperatures,  shows  that  the  method  of  calculating  component 
enthalpies  by  adjusting  ideal  enthalpies  with  heats  of  mixing 
is  quite  adequate. 

A  comparison  of  the  optimization  work  reported  in 
Figures  8.4  and  9.4  shows  that  the  optimum  operating  condi¬ 
tions  in  the  phenol-MCH-toluene  system  occur  in  a  much  smaller 
range  of  reflux  ratio  than  in  the  acetone-methanol-water 
system.  Thus,  the  optima  in  the  MCH-toluene  system  might  be 
considered  only  local  optima  valid  for  that  short  range  of 
reflux  ratios  only,  and  that  some  higher  combination  of  reflux 
and  feed  ratios  might  be  a  universal  optimum.  Although  this 
is  possible,  the  reflux  ratios  in  this  system  were  investi¬ 
gated  up  to  L^/D  values  of  30.0  without  such  a  combination 
occurring,  and  it  is  doubtful  whether  higher  values  of  L  /D 
would  be  of  any  practical  interest. 
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XI  -  SUGGESTIONS  FOR  FUTURE  WORK 

Future  work  in  this  field  should  proceed  the  following 
directions : 


1.  The  development  of  an  equilibrium  model  where  the 

parameters  used  can  be  derived  from  intrinsic  proper¬ 
ties  of  the  components  in  the  system.  In  the  Wilson 
model  and  X^  are  only  'proportional'  to  the  i-j 

and  i-i  interaction  energies  respectively,  and  the 

interaction  energy  parameter  (X^j  -  X^)  cannot  be 

4  10 

rigorously  defined  ’  ,  so  that  the  main  concept  is 

not  based  on  any  fundamental  properties  of  the  compo¬ 
nents  i  and  j.  The  development  of  a  model  where  this 
is  the  case  would  have  wide  applications  in  chemical 
engineering  in  general  and  extractive  distillation  in 
particular . 


2.  Investigate  the  possibility  of  using  the  parameter 

(X.  .  -  X.  .)  to  determine  the  best  solvent  for  the 
lj  ii 

extractive  distillation  of  a  system.  Since  we  know 
that : 


a)  The  action  of  a  solvent  is  to  reduce  the  inter- 
molecular  bonds  in  a  system  thus  leading  to  more 
mobility  and  a  higher  volatility,  and  that 
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b) 


the  1-2  interaction  energy  increases  as  stronger 

4  10 

molecular  and  hydrogen  bonds  exist  in  a  binary  ' 


It  should  then  be  possible  to  predict  solvent 
action  a  priori  by  deciding  on  the  relative  magnitudes  of 
the  1-2  interaction  energies  between  the  solvent  and  the 
light-key  component  and  between  the  solvent  and  the  heavy- 
key  component.  The  wider  the  difference  between  these 
magnitudes,  the  more  successful  will  be  the  solvent. 

If  this  idea  is  feasible,  and  if  the  concept 
( X i j  can  be  based  on  a  more  fundamental  basis,  choosing 

a  solvent  for  extractive  distillation  would  lose  the  element 
of  art  and  experience  in  it  and  tend  to  be  based  on  basic 
principles . 


3.  Develop  a  deeper  insight  into  the  predictive  properties 
of  the  interaction  energy  parameters  (A^“  A^)  .  Our 
work  shows  that  widely  different  magnitudes  of  this 
parameter  are  able  to  predict  vapor-liquid  equilibria 
equally  well.  This  could  be  due  to  any  of  two  reasons 
or  their  combination: 


a)  Minimum  standard  deviation,  which  is  the  parameter 
we  have  used  to  decide  on  the  predictive  ability  of 
the  interaction  energies,  is  an  insensitive  indication 
of  how  well  two  vectors  compare. 
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b)  Different  magnitudes  of  the  interaction  energies, 

which  cannot  be  a  priori  calculated,  are  inherently 
able  to  predict  adequate  equilibrium  data. 

The  question  then  arises:  Is  it  possible  in  a 
ternary  system  for  example,  to  use  two  pairs  of  interaction 
energies  at  their  optimum  values  and  one  grossly  inaccurate 
pair  and  still  predict  good  ternary  vapor-liquid  equilibria? 

If  this  is  the  case,  then  equilibrium  data  on  each  one  of 
the  binaries  involved  in  multicomponent  systems  might  not  be 
necessary.  This  would  be  a  milestone  in  equilibrium  prediction 
with  very  limited  data. 

4.  At  the  present  stage,  a  tedious  literature  search  is 
often  necessary  to  locate  data  on  the  ideal  enthalpies 
of  the  components  being  studied  and  the  binary  heats  of 
solution  at  infinite  dilution.  The  development  of  a 
generalized  liquid  enthalpy  correlation  that  will  take 
enthalpy  departures  due  to  heats  of  solution  into  account 
would  be  very  helpful  in  extractive  distillation  design. 
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DEFINITION  OF  SYMBOLS 


constant  in  various  expressions 

constant  in  enthalpy  correlation  of  Chao-Seader 
constant  in  various  expressions 
second  virial  coefficient 

constant  in  enthalpy  correlation  of  Chao-Seader 
number  of  components 

constant  in  enthalpy  correlation  of  Chao-Seader 

constant  in  enthalpy  correlation  of  Chao-Seader 

reference  fugacity:  fugacity  of  the  pure  saturated 
component  as  a  liquid  at  zero  pressure 

excess  Gibbs  energy  of  mixing 

feed  rate,  moles/U.T. 

the  fugacity 

enthalpy  in  the  vapor  phase,  BTU/lb.mole 

enthalpy  in  the  liquid  phase,  BTU/lb.mole 

vaporization  constant 

liquid  flow  rate,  moles/U.T. 

partial  molar  heat  of  solution 

number  of  stages 

number  of  components 

saturation  vapor  pressure 

heat  load,  BTU/U.T. 

liquid  sidestream  flow  rate,  moles/U.T. 

vapor  flow  rate,  moles/U.T. 

3 

molar  liquid  volume,  cm  /mole 
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vapor  sidestream  flow  rate,  moles/U.T. 

liquid  composition,  mole  % 

vapor  composition,  mole  % 

feed  composition,  mole  % 

compressibility  factor 

activity  coefficient  in  liquid  phase 

c 

fugacity  coefficient  at  P 
fugacity  coefficient  in  vapor  phase 

proportional  to  interaction  energy  between  i  and  j , 
cal/gm.mole 

Wilson  parameter 

f/p  ratio  in  liquid  phase 

acentric  factor 

solubility  parameter 
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APPENDIX  I 

EQUATIONS  FOR  THE  OPTIMIZATION  OF  THE 


INTERACTION  ENERGIES  IN  THE  WILSON  MODEL 
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APPENDIX  II 


COMPUTER  PROGRAM  AND  OUTPUT  FOR 
OPTIMIZED  INTERACTION  ENERGIES 
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$  I P  FTC  WILSON  NOLTST *NODEC< 
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C  CALCULATION  OF  OPTIMUM  binary  INTFRACTION  energies 
C  USING  THE  TECHNIQUE  OF  QUAZ ILINEARI7ATION 

c 

IPO  FORMAT  H  X  ,  I  ?  ) _ 

10?  FORMAT ( 1 X*?l HINITIAL  GUESS  A  =,Fl0.4*RH  ,  R  =,F10.4) 

10s  FORMAT (1X»1bHITERATI0N  NO.  ,  I  b  ,  ?X  ,  ^>HA  =*F10.4*rH  ,  P  =,Fi0.4) 

_ 104  FORMAT ( 1 HL ) 

1 01  FORMAT ( 1 X ,4 FT  ? .5 ) 

105  FORMAT ( IX *64H  X  ORS  Y  OBS  Y  CALC 

_ 1  ARS  FRRQR) _ 
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C  NSFTS  TS  NO.  OF  SYSTFMS _ 

C  FMT  IS  NAMF  OF  SYSTEM 

C  FMW  IS  MOL  WT.  OF  COMPONFNT  1 

c  SMW  IS  MOL  WT*  OF  COMPONENT  ? 

C  N  IF  NO.  OF  data  POINTS  IN  THF  SYSTEM 

C  RHO,  ALPHA,  and  BFTA  ARF  CONSTANTS  IN  the  DENSITY 

C  EQUATIONS  OF  THE  Rf.sPFCTIVF  COMPONENTS _ 

C  A  amo  B  ARE  THF  INITIAL  GUFSSFS 

C  Y »  X  »  C  *  AND  T  ARF  THF  DATA  POINTS  WHERE  C  1 9  GAMMA 
C  M  IS  a  COUNTER 

C  AS  FAR  A S  POSSIBLE*  THF  SYMBOLS  USED  ARF  THF  SAME  AS 

C  IN  THE  EQUATIONS  THEY  RfPRESENt  IN  THE  PREVIOUS  APPENDIX 
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PDAFH (  T  )  =  (  1  .0-X  1  (  T  )  )*  ( 1  .O/FTA  (  I  )  -1  .0  /  (X  1  (  I  )  +RFTA  (  I  )  )  ) 

F0A<^H( T ) =ALPHA < I ) /FT A ( T ) 

Y0  ALC  (  I  )  =C  (  I  )  *X  1  (  T  )  *  FDA  FH  (I  >  *FXP  (  PDAFH  (  I  )  ) 

PRROR (T)=APF(YCALr(I)-Vl(T)) 

WRlTF(A*10O)Xl (  I  )  , Y1  (  I  )  , YCALr ( I  )  •  F  RROR (  I  ) 

P6  CONTINUE _ _ _ 

WRITE (6 ♦ 1 ) 

11  F I GQFQ=0 .0 

- IGPR»0.0 _ 

c  !  GPY1  =n  .n 
FIGPF0  =  0  .o 

g IGOR=o  .  r> _ _ _ 

g IGQYi =n  ,o 
C I GPO=o • o 

OQ  12  L*J  _ 

ALPHA ( T ) = ( V 1 ( I  )  /  V? (I) ) *FXP ( A/ ( RR*r ( t ) ) ) 

P  F  T  A ( I )  =  ( 1 .0-X 1(1)  ) * ( V? ( I )/Vl  (T)  )*FXP ( R / ( R  R  *  T ( T )  )  ) 

FT  A  (  I  )=X1  (  I  )»,ALPHA(  I  )+M  .0-X1  (I)  ) _ _ 

P  D  AS  H ( I  )a(l«0-Xl(I  M:»C  1«0/ETA(  I  >-l«n/.(Xl  (I  )+PFTA(  I  n  > 

SOASH( T  > =  A L P H A (  I  )  /  FT A ( T  > 

YCALC (  I  ) =C (  I  ) *X]  (  T  ) *  c  n  A  FH (  I  I *^XP (PDAFH (  I  )  ) 

7  f  t  a  (  y  )  =  ^XP  (  PnA  SH  (  T  )  )  *  aLPh A  (  !  }*(  i  .  o-x  1  (  I  )  )  /  (RR*T  (  I  )*FTA  (  I  )**2  ) 

THFTA ( I  )  =  ( 1 • 0-X 1  ( T )  >*RFTA  < I >*£XP ( PDASH< I  )  ) / ( RR*T ( T )*( XI  ( I  )+RFTA { I 
I*»?  )  ) _ _ _  _ 

P (  I  )=C(  I  )  *X  1  (  I  )  *  Z  E  T  A (  I  )  * ( 1 • 0-X 1  (  I  )  *  S  D  A  S  H  ( I )  ) 

Q(  I  )=C(  I  )  *  X 1  (  I )  *  g  D  A  F  H  ( I ) *THF  T  A (  I ) 

p(  i  )=vrALr(  t  )-a*P(  t  )-p*g(  i  )  _ _____  _ _ 

S T  GQSG*S T  GO SQ+Q ( T )**P 
FTGPR=FIGPR+P( I )*R( I ) 

F  T  GPY1  =  F  I  GPY1  -f  P  (  I  ) *Y 1  ( T  ) _ _ _ _ _ 

SIGPSQ=s TGPSQ+P (  I  ) **? 

F  I  GQR  =  F I GQR+Q ( I ) *R ( I ) 

SI GOY 1«SIGQY1+QU  )*Y1  ( I  ) _ 

c TGPO=F IGPQ+P (I ) *Q( I) 

IP  FONT  I  NUF 

Pi  =  (  SIGPSQ*SIGQR-SIGPSQ*SIGQYl  +  STGPQ*5  TGPY1  -  F  TGPO*.FTGPR  )  /  f  FTRPO**: 
1-F IGPFO*FIGQFQ) 

A  1  =  ( SIGPY1-R1*SIGPQ-F IGPR ) /figpfg 
WR  I TE  (  6*10  3L )  M  VaI 

IF(ABF(A1-A).LE.1.0F-0?)G0  to  16 

A  =  A  1 _ _ _ _  _ _ 

q  =  R  1 

GO  TO  11 

1£  IF(ABS(B1^B)»LE»1>OE-0^1GO  to  if 
A  -  A  i 
p  =  P1 

GO  TO  11 _ _ _ _  _ 

15  '*'R  I  tf  ( 5  *  1  04  ) 

WR  I  TF  (  6  ♦  1 05  ) 

/,n  nr>  20  T  =  1.‘' _ _ _ _ 

ERROR  (  I  )  e ABF  ( YCALCU  )  -Y1  (11) 

URITF(6»10Q)X1  (  I  )  *Y1 ( I  )  ,YCALC( T )  .FRROR (  I  ) 

20  CONTINUF  _  _ _  _ 

WR  I  TE ( 6  » 106 ) 

Fn  CONTINUF 
FND 
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ACETONE-WATER  SYSTEM 


INITIAL  GUESS  A  =  466.0000 


3  =  1653.0000 


I TERA  T I08 

NO. 

1 

A 

= 

545. 1388 

t 

3 

as 

1595.2224 

HERAT  1 0\ 

NO. 

2 

A 

= 

566.4876 

y 

8 

s 

1582.7030 

ITERATION 

NO. 

3 

A 

= 

569.6239 

y 

B 

~ 

15/9.8546 

ITERATION 

NO. 

4 

A 

= 

5  70.4456 

t 

B 

= 

1579.1616 

IT  ERA T ION 

NO. 

5 

A 

5  /0.6365 

t 

6 

= 

1578.993  / 

ITERATION 

NO. 

6 

A 

= 

570.6831 

y 

B 

- 

1578.9529 

! TERA  T I  ON 

NO. 

7 

A 

= 

5  70.6947 

y 

B 

= 

1578.9429 

ITERATION 

NO. 

8 

A 

- 

5  70 . 69  /3 

t 

8 

= 

1578.9406 

ITERATION 

NO. 

9 

A 

s 

5  /0.69  8  1 

r 

B 

s 

1578.9400 

X  03S 

0. 10CC0G00 
0.20000000 
0. 300  0  00  0  0 
0 . 4 GO 0 00  GO 
0. 50C000  00 
0. 6000 CO  00 
0. 7000C0  00 
0 . 80 COCO  00 
0 . 9CC0C0  00 


Y  OBS 

0. 76000000 
0.79000000 
C. 81500000 
0 . 83000000 
0. 84C00C00 
0.86000000 
0. 87000000 
0.90000000 
0. 94000C00 


Y  CALC 
0. 78132080 
0 . 74836183 
0.76713257 
a  0.82070781 
0.83688446 
0.91405127 
0.91921479 
0.91942493 
0.93908454 


ABS  ERROR 

0.021 32080 

0.0416381/ 

0.04786743 

0.00929219 

0.00311554 

0.05405127 

0.049214/9 

0.01942493 

0.00091545 


ACETONE-METHANOL  SYSTEM 


INITIAL  SUESS 


ITERATION  NO. 
ITERATION  NO. 
ITERATION  NO. 
ITERATION  NO. 
ITERATION  NO. 
ITERATION  NO. 
ITERATION  NO. 
ITERATION  NO. 
ITERATION  NO. 


X  OBS 

0. LCOCGOOO 
0.20000000 
0 . 3C0CCU  00 
0.4000  OU 0  0 
0 .50000000 
0.60000000 
0. 70000000 
0.80000000 
0.90000000 


A  =  -220.9000  ,  B  =  596.0000 


A  =  -250.7115 
A  =  -249.7374 
A  =  -2 50.2654 
A  =  -250.3610 
A  =  -250. 39  o5 
A  •=  -250.4070 
A  =  -250.4116 
A  =  -250.4127 
A  =  -250.4L30 


6  =  608.2608 
B  =  610.7667 

8  =  611.6919 

b  =  611.9687 

8  =  6L2.056L 

6  =  612.0833 

0  =  612.0926 

b  =  612.0954 

B  =  612.0963 


Y  OBS 

0.  13  000  000 
0. 30000000 
0.42000000 
0.51000000 
0.60000000 
0.67500000 
0. 73COOOOO 
0.80000000 
0. 89000000 


Y  CALC 
0.  1974  3172 
0 . 29693533 
0.41453840 
0.49012692 
0. 58743420 
0.69760516 
0 . 74952030 
0.80860055 
0.91286954 


ABS  ERROR 
0.01743  L  72 
0.00306467 
0.00546  IbO 
0.0198/308 
0.01256580 
0.022608  1  6 
0.0  1952031 
0.00860056 
0.02286953 
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APPENDIX  III 


PROBLEM  STATEMENTS,  INITIAL  AND  FINAL  PROFILES, 
AND  PRODUCT  COMPOSITIONS 


FOR  PROBLEMS  I  TO  IV  INCLUSIVE 


hi  xiawaqqA 


^sajiaoaq  jAma  q*ia  jaitiwi  ,  2TnaM3TAT2  Majaoaq 
awoiTiaoqjMOD  Touaoaq  awA 


3VI8UJ0MI  VI  OT  I  2MajQOflq  H03 


PRORLFM  NO.  1 
ONF-F^FH  FRACTIONATOR 
PRORLFM  STATEMENT  AND  DATA 


A.  COLUMN  DATA 

PRESSURF  =  ? 50,0  PSIG 
NO.  OF  PLATE<;=  16 

conqfnsfr  typf=  partial _ 

OVERHEAD  PROOUCT=  22,6  MOL  VAPOR/ll.T. 
ROTTOM  PRODUCTS  77.4  MOL/U.T. 

RrCLUx  RATF=  154.3  MOL/U.T. 

FFED  PLATE  NO.®  P 


R.  FFFD  DATA 

TFMPFRATIlRF=  ??5  F 

RATF=  1^0  MOL/U.T. 

VAPOR  C0NTFNT=  16.4  MOL/lOO  MOL 

composition  in  mole  fr act  ton 

COMPONENT  COMPOSITION 

C2 

O' 

c 

. 

c 

CR 

0,20 

C  4 

0  •  P  7 

CP  0  •  3  F 

C  6  0,^5 


C.  EQUILIBRIUM  Ai'MD  ENTHALPY  DATA  ARE  POLYNOMIAL  TEMPERATURE 

FUNCTIONS  RFPORTFD  pv - 

AMt INDSON  ,  N  . R  .  ,  AND  PON  T  I NEN  ,  A  .  J.  , 

7^n,  ( 1  o*P 8  ) 


IND.  FNG.  CHFM.,50, 


r  .  '  ’  j 

'  POT  *  rtO! 

T  •  /’HA  TH  ■  T  a  T?  •>  JMOl-  9 


I  ’  ’  i 

■"> t  .  ;5<;  *qfuj-v.*qq 

=?^TAJ°  1 
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°  irajq  nr.1 


*  T  AH  IT 
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PROBLEM  NC.  1 


INITIAL  PROFILE 


STAGE _ TEMP  VAPOUR  RATE  JL I QU  I .0  RATE 


F 

MOLES/U* 

T. 

MOLES/U. 

T. 

1 

1 1C.  CCO 

C. 226CC00E 

02 

0. 154300GE 

03 

2 

12  3  .CCO 

C.  1765CCCE 

03 

0. 1 543000E 

03 

■3 

136  .COO 

C.  1 769C00E 

03 

0. 1543000E 

03 

4 

14  5  .  C  C  C 

0.  1  769CGCE 

03 

0. 1543000E 

03 

c 

162 .CCO 

C. 1769C0GE 

03 

0. 1543000E 

03 

6 

175. CCC 

C.  1765C00E 

03 

0. 1 5  4  3  0  0  0  E 

03 

7 

188  .CCO 

•  C.  1769C00E 

0  3 

0. 1543000E 

03 

8 

201.  CCC 

C.  1  769C0CE 

03 

0. 2379000E 

03 

S 

214. CCO 

C. 1605CC0E 

03 

0.237900CE 

03 

1C 

227. CCO 

0. 1 6O5C00E 

03 

0.2379000E 

03 

1  1 

24C  .CCO 

C.  1605C00E 

03 

0.2379000E 

03 

12 

253. CCO 

0. 1605C00E 

03 

0.2379CC0E 

03 

12 

266 .COO 

C. 16C5CC0E 

03 

0 . 2379000E 

03 

14 

275. CCO 

C. 1605000E 

03 

0.2379000E 

03 

15 

292 .CCC 

C.  1605C00E 

03 

0.2379000E 

03 

16 

305. CCO 

0. 1605CCCE 

03 

0. 7740CG0E 

02 

FINAL  PROFILE 

4 


STAGE 

TEMP 

VAPCUR  RATE 

LIQUIC  RATE 

F 

MOLES/U. 

T. 

MCLES/U. 

T . 

1 

118.921 

0. 226C000E 

02 

0. 1543000E 

03 

2 

125.553 

0.  1769CC0E 

03 

0. 1572289E 

03 

•2 

137. 8C2 

0. 1798239E 

03 

0. 1548853E 

03 

4 

146 .565 

C. 1774853E 

0  3 

0. 1495772E 

03 

c 

157.638 

0. 1721772E 

03 

0. 1424170E 

03 

6 

171.470 

0. 165C17GE 

03 

0 . 1336754E 

03 

7 

188.842 

0.  1562754E 

03 

0.  123255  IE 

03 

8 

211.049 

0.  1458552E 

03 

0.207779 7 E 

03 

c 

> 

218.629 

0.  1  303797E 

03 

0.211 3498  E 

03 

1C 

225.597 

0.  1  339498E 

03 

0.213912 9 £ 

03 

1  1 

232.279 

C.  1365129E 

03 

0.2163277E 

03 

12 

238.656 

C.  1389277E 

03 

0.2185323E 

03 

12 

244.926 

C.  141 1323E 

.  03 

0.2  199648  E 

03 

14 

251. 88C 

C. 1425643E 

03 

0.2197723E 

03 

15 

261 .242 

0. 1423723E 

03 

0.217 1495E 

03 

16 

275.303 

C.  1397495E 

03 

0. 77400C0E 

02 

SUM  (C6LTA  T**2)  AT  ITERATION  A 


0.068 


■ 


- 


PRODUCT  STREAM  CGMPGSI TIG MS 


COMPONENT  CONDENSER  STREAM 

C0  C.5048424E-01 

C3  0.9136172E  CO 

C.  0.35851  19E-C1 

4 

C5  _ 0 ,47406 5 3E-Q4 


BOTTOM  STREAM 
0.  1961257C-06 
0. 10G1785E— 01 
0,473 59 12E  00 
0, 451841 7E  00 


C. 1407975E-07 


0. 645491  IE-01 


' 


FINAL  PROFILES  AS  REPORTED  IN 

LITERATURE  ( 

REF.  26  ) 

_ STAGE 

.  TEMP  _ 

vapour  Rate _ 

LIQUID  RATE  _ _  _  _ 

F 

MOLES/U.T. 

MOLES/U.T 

• 

1 

119.020 

0.2260000E  02 

0.1 543000E 

03 

2 

130.140 

0.1769000E  03 

0.15  5  5400E 

/ 

03 

3 

138.090 

0.17814 OOE  03 

0. 1525400E 

03 

4 

146.950 

0.1751400E  03 

0 . 1470500E 

*  V 

03 

c; 

158.100 

0  •  16965 OOE  ^3 

0.1399200E 

•  > 

03 

<  * 

6 

171.940 

0 . 1625200E  03 

0 .1313000E 

n  3 

7 

189.260 

0.1539000E  03. 

0.1209000E 

03  ' 

8 

211.520 

0.143  5000E  03 

0 • 2064100E 

03 

9 

218.820 

0 . 1290200E  03 

0.2090600E 

03 

10 

225.660 

0.13 1 6600E  03 

0.2112900E 

03 

11 

232.290 

0.133  8900E  03 

0.21361 OOE 

03 

12 

238.660 

0 •  136  2 1 OOE  03 

0.2 158300E 

03 

13 

244.930 

0 . 1 3843 OOE  03  .. . 

0 .2173400E 

03  .  ...  „ 

14 

251.870 

0.1399400E  03 

0.2172500E 

03 

15 

261.190 

0. 1398500E  03 

0.2147600E 

03 

16 


275.710 


0.1 373600E  03 


0.7740000E,  n2 


nS' 


PRORLFM  NO.  2 
TWO-FFFD  FR ATT  I ONATOR 
PRQPL^M  RTATFMFNT  AND  DATA 


A.  column  data 

PRrcQitRF=  2  c  A  .  o  PSIG 
NO.  OF  PL ATES=  21 

C ONDFNSFR  TYP^=  PARTIAL _ 

FIRST  FFFD  RATF=  24.5  MOL  of  saturated  VAPOR/II. t. 
FIRct  fffd  PL  a  tf=  7 

SFCOND  FEFD  R ATE=  6^.^  MOL  OF  SATURATED  L  IQlJ  I  Q/i  I .  J  . 
■COND  fffd  PLATE*  1 3 

FIRct  F I DF  STREAM=  15.0  MOL  OF  LIQUTD/M.T. 

F  I R  ft  cjqfftREAM  Platf=  4 _ _ 

SFFOND  S  I  DFSTRFAM=  25.0  MOL  OF  VAPOR  /l !.  T  . 

SFCOND  fidESTRFAM  PLAtf=  16 

OV^RHFAD  PRODUCT**  91*0  MOL  \/APOR/m.t.  _ _ _ 


REFLUX  RATE 

=  7] .806  M0L/U.T. 

R.  F^Fn  DATA 
FIRST  FFFD 
SFCOND  FFFD 

TFMPERATURF=  117.7 
TFMPFRATURF=  268. 

F 

5  F 

COMPONENT 

FIRST  FFFD 

COMPOS  I T ION 
MOLFS/II.  T. 

FFCONO  F  F  F  n 

COMPOS  I T ION 

MOLFS/U.T. 

METHANE 

?  .000 

o  .  ooo 

ETHANE 

1 o . ooo 

o  .  noo 

propylene 

6  .  ooo 

,  .l.ooo 

propanf 

1 o . non 

7.ooo 

I  cO-pi  itanf 

1  .ooo 

4.^00 

n-rutanf 

5  .  OOO 

1  7  .  Ann 

N-PFNTANF 

o.5oo 

15.200 

n-hfyanf 

0.000 

9  .  ono 

HFPt ANF 

0.000 

4.5oo 

Oct  a  n  f 

o  .  ooo 

4.300 

OtHFRf (  26o 

F )  o.onn 

2.500 

34  « 500  65.5  OO 


C.  EQuTlT RR 1 1  IM  AND  ENTHALPY  data  ARE  Polynomial  TEMPERATURE 
FOOT  I ON S  REPORTED  RY - 

LYSTER.  U.N.  ,  ET  AL.  ,  PFT.  REF.  08%  No.  6.  _?  ?  1  ,  NO.  7% 

151,  AND  NO."  I  o’,  130  (1Q5S)' 
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•  . 
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r  r  *  TAjq  ni^i 

•  T«i»\ot  tj  ii  j*  .a  *  -1^-?  imi?  T->,fq 

J>  e  IT/  J  '  'AM  T?  ‘Jn!  5  T  *>  -1  1  3 
•T."\-  '  '  ',/  -J(t  jr  ‘  ",  r  -  A  ^  ‘  "j  n  j  3  ’  '  M5 

Ar  =  TAJG  A  iVr?  -'lp  AM03‘3? 

.t.io  :  G/y  i  M  -  .  r  r  =  T3U'~0Vjq  AlMq?\/0 

.7  .  1  *\  •  M  a  Y-  nqjj 


UAfi'  A  -  M 

1  X 

.  v  r  r  =  R  q  r  '  v  m  •■  m  t 

1  q  3  3  t  ’ 1  •  3 

1  A. 

t  -i it a* qqyqr 

o  ^03  is 

oq^q  hmooip 

Tl43|/<0  l  •  3 

Mom^oq  03 

MO T  T  T  POqMOO 

.  T  .  \  '  ij 

.  .  \ 

AAA#ft 

OAO.  r 

AAft.r 

AAA  f 

^MAHt3 

.  r 

AAA  .  > 

■RM^Jyqo^q 

AAA.T 

oaa.o  r 

3/uq  >qq 

A  AA  #  ^ 

oo a  #  r  . 

T  '  -•  - 

AAA  #  -  (• 

A  A  A  .  c 

qn  a  y  'o-i/ 

.  f 

OOP  #  A 

o  T  ^ r  *— 1/1 

A  A  A  #  p 

ooo  .  A 

3M  A  v3H— M 

aa?  .  A 

AAA  .  O 

AT93H 

OAT  .  4\ 

OAO  ,  o 

VM  T30 

.  r 

oo A  *  A  f  q 

a  ^  r >pqqHT0 

or  .  P* 

A  A  f*  ,  *’ 

if  JAIK.OJ--.YJ  3HA  ATAO  Y^JAhTH  J  Ci  a  MU  fflfl  I  J I  UU  =| 

- V  'Ml  MTV  ■  'OlT'iiq 

•TA-C  .  ,  r  V\  ,*r  ,  .  :yy.  ,t:  .  .  T  ...'•*  ,fl  -JMYJ 

(  -  'or  )  •  r  ,0  r  .0  :  a  .  r  >  f 


INITIAL  PROF  I  L  f 


STAGE 

TFMP 

VAPOUR  1 

RATF 

’  LIQUID  RATF 

F 

MOL F  F/i  i 

x 

•  • 

MOLFS/U. 

T 

I  • 

1 

8°.°o 

^ . 77000F 

07 

0 .7 1 8  1  AF 

0  7 

2 

98.50 

°  .  °4  8 1  OF 

02 

0.71 81  OF 

02 

3 

117.00 

0.040 ] nF 

07 

o.71 8 1 0  F 

07 

4 

135.50 

0. P48 10F 

07 

0.56 81 OF 

02 

3 

154.00 

0.94  8 1  Of 

07 

0.56810F 

07 

6 

172.50 

0,048 1  OF 

07 

0.66 81  OF 

02 

7 

101.00 

0.0401  OF 

02 

0 . 568 1  of 

07 

8 

7  °9 . 5  0 

n,6oo ] nF 

02 

0.5681  OF 

02 

Q 

228.00 

0.607 1 oF 

07 

0 .568 1  OF 

07 

r° 

246.50 

0.607 i or 

07 

0.668 1  of 

07 

1 1 

765.00 

0.6O31OF 

0? 

0.568] OF 

07 

l? 

283.50 

0.603 1  OF 

02 

0 .56810F 

07 

13 

283.50 

0.603 1  OF 

07 

0 . 1  ?  7  3  ]  F 

03 

14 

32°  •  50 

0.603 1 OE 

07 

o.l 7231F 

07 

15 

039.00 

o.603 1  OF 

02 

o.i  223  IE 

03 

16 

057.50 

0.607  1  OF 

0  7 

o.l  2  7  7 1 F 

07 

17 

776. °0 

0.853  1  OF 

02 

o.l  7  7  3 1 F 

03 

18 

094.50 

0,857  1  OF 

07 

0.17231F 

03 

IQ 

410.00 

0.857100 

0  7 

0.1  7  7  3 1 f 

07 

20 

431 .50 

O.85710F 

02 

0 . 1  2  2  3 1  F 

03 

21 

450.00 

0.857 1 oF 

02 

o.l2?31F 
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FINAL  PROFILE 


STAGE 

TFMP 

F 

VAPOUR  RATF 
MOLFS/U. T. 

LIQUID  RATF 

MOLFS/U.T. 

1 

7  7.38 

0 . 73O00F 

0? 

0 .71 81  OF 

07 

2 

94.87 

0.94810F 

02 

0 .7776oF 

02 

3 

103.69 

0.95760F 

02 

0.71510F 

o2 

4 

111.11 

0.945  1  OF 

02 

0.54620F 

02 

5 

120.75 

0 . 976  20F 

07 

0.57490^ 

07 

6 

131.14 

0 . 904 9 of 

07 

0 . 5  06  3  OF 

0  7 

'  T 

141.21 

0.886 7 op 

0  7 

0.48300F 

0  7 

8 

1 69.09 

0.51  ROOF 

0  2 

0 .49740F 

02 

Q 

185.9° 

0.53740F 

02 

°.c0] 1  OF 

02 

IO 

198.05 

0.536] OF 

07 

0 .4967oF 

07 

11 

208.71 

0.53] 7  OF 

02 

O.47R30F 

0? 

12 

277.68 

0.51770F 

07 

0.475005 

07 

TT~  ~ 

753.79 

0 .460O0F 

07. 

o.l.066°F 

07 

14 

757.65 

0 . 4  4  6  9  0  F 

02 

o.l 0872F 

03 

15 

267.79 

0.46 7 7 OF 

02 

0.1 0925F 

03 

16 

769.56 

0.47750F 

07 

0 . 1 I006F 

0  7 

17 

779.70 

0 . 73O60F 

07 

O  .  1  1  0  0  4  F 

07 

18 

283.97 

0 . 73040F 

07 

0.10918F 

03 

T9' 

794.58 

0.77180F 

02 

0.1O715F 

03 

20 

311.34 

0.701  5  0  F 

07 

0.1o]94F 

03 

21 

744.73 

0.649 4 OF 

07 

0.77O0OF 
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PRODUCT  CTRFAM  COMPOSITIONS 


CO'  •  v  -  CONDFNSFR  RECYCLE  ROT  TOM  STREAM 


\’:r  THANS 

o . OOF  P 

n, aaoo 

l-  *J*  H  A  iV 

o . noon 

?  R  C  P  y  L  r  N  F 

n. 

0.0000 

DROP AN r 

o  •  5  3  6  4 

o.ooo7 

I  SQ- Pi  IT  M  F 

o  .  roof 

0.017? 

V  _  p  u  T  a  \>  F 

0 . 0  0  5  * 

0.1471 

N  —  P  c  N  T  A  N  - 

r.oon^ 

" . 0^63 

n-mfx anf 

n .  on  Do 

0.7131 

HF PT  A\ir 

n  .  0000 

0.1141 

OCTANF. 

o  .  000O 

0.1127 

OTHPRS(S6A  F) 
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FINAL  PROFILES  and  COMPOSITION  as  REPORTED  in  LITERATURE  (REF.  26) 


STAGE  TEMP  _  VAPOUR  RATE 


F 

MOLES/U.T. 

1 

77.290 

0.2300000E  02 

2 

94.770 

0.O48O600E  02 

3 

103.610 

0.9534200E  02 

4 

111.070 

0.9461700E  02 

5 

...120*780. 

....  Q. 9.2.7.13 00E  .02 

6 

131.270 

0.9053700E  02 

7 

141.410 

0.8869100E  02 

*■  \ 

8 

169.450 

0.52^6900E  02 

Q 

186.400 

0.5386300E  02 

V  > 

10 

198.520 

0.5424200E  02 

11  . 

209.140 

0. 5382700E  02 

12 

222.910 

0.5 199700E  02 

13 

253.420 

0.46 5 9000E  02 

< 

14 

257.560 

0 . 4447900E  02 

< 

15 

262.51 0 

0.4580O00E  02 

. 

16 

269.050 

0 . 46 8 1 400E  02 

.17. _ 

275.440 

0.7265200E  02. 

18 

282.950 

072658000E  02 

19 

293.300 

0771 80500E  02 

20 

309.880 

076977500E  02 

21 

343.330 

064589000E  02 

COMPONENT 

LIQUID 

VAPOUR 

SIDESTREAM 

SIDESTREAM 

methane 

0.26010E-01 

0. 19480E-10 

ETHANE 

0.88760E  00 

0.2081 0  E— 04 

PROPYLENE 

0.29773E  01 

0.21480E  00 

propane 

0.88691E  01 

0.19773E  01 

I -BUT  ANE 

0 * 994 1 OE  00 

•  0.32555E  01 

N-BUTANE 

0.12404E  01 

0.13705E  02 

n-pfntane 

0.54800E-02 

0.42934E  01 

HEXANE 

0.41940E-06  . 

0.11099E  01 

HEPTANE 

0.34570E-09 

0.28010E  00 

OCTANE 

0.24280E-12 

0.13390E  00 

OTHERS ( 360 ) 

0.16310E-18 

0. 29660E-01 

MA3*T23Qie 


*n-3onns.o 


10  3X9*8F.O 


<•  -  /-  .  .  «  *-«.  _ 


PRORLFM  NO.  7 
TWO-FFED  FRACTIONATOR 
proplfm  FTATFMFNT  and  data 


a.  COLUMN  DATA _ 

PRESSURE*  750.0  psig 
NO*  OF  PLATFs=  21 

condfnsfr  typf=  partial _ _ _ 

FIRST  FEED  RATE=  74.5  MOL  OF  SATURATED  VAPOR/!!. T. 
FTRft  FFn  PL ATF=  7 

SECOND  FEED  RATE*  65.*  MQL  OF  SATURATED  L  IQi  i  I")/’ : .  T  . 
SFCOND  Fffo  PLATE*  17 

FIRST  FIDF  STRFAM=  15. n  MOL  OF  LIOUID/U.T. 

FIRST  CIDFSTRFAM  PLATP=  4 _ _  _ _ 

SFCOND  SIDESTREAM*  76.0  MOL  OF  VAPOR/U.T. 

SFCOND  SIDESTREAM  P  L  A  TF  =  16 

OVERHEAD  PRODUCT*  22*3  MQL  VAPOR/U.T. 

REFLUX  Ratf=  71 • 806  MOL/U.T. 


n.  fffd  data 

FIRST  FrFn  TFMPFRAT!  iRF  =  117.7  F 


SFCOND  FFFD 

TFMPERATURF=  268.5  F 

COMPONENT 

FIRST  FEFD  SFCOND  FFFD 

COMPOSITION  COMPOSITION 

MOLFS/ii.T. MQLFS/U.T. 


mfthanf 

ethane 

PROPyLFNF 

?  .  ooo 
l^.onn 
6.000 

o  .  a  no 

o  .  O OO 

]  .  PRO 

PROPaNF 

I  SO  — RUT  AN  F 

N-R|  ITANF 

1 o . ooo 

1  .odd 

7.000 

7  .  ooo 

4  •  n  o  o 

17.  ^o 

N-PfNjANf 

0.500 

16.200 

n-hfyanf 

0.000 

O.noo 

HFPtaNF 

0.000 

4  •  5  RO 

OCTANE 

o .  ooo 

4.300 

OTHFR  S  (  360  F) 

o  .  ooo 

3. 5oo 

34.500 

65.500 

C.  EQUILIBRIUM  AND  ENTHALPY  DATA  ARE  CHAO-SEADER  DATA 

THE  CONSTANTS  FOR  WHICH  ARE  REPORTED  IN  THE  NEXT  PACFS 
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component  Physical  propfRttfc 
FOR  CHAO  SFAOPR  CORRF L A T  T ON 


COMPONFNT  c 

TC 

PC 

W 

delta 

VI 

H2C 

67?. 48 

1 806. 6 

o. 0868 

8.8oo 

43  .  ] 

nitrogfn 

_ 227. 16 

.492.3  . 

.0.0  266. 

.  4.440 

_  5.3  .n  . 

CO? 

647  56 

1 071.8 

0.1768 

7.120 

44.o 

METHANE 

343 . on 

67^  .  1 

n • 

5.680 

5  ?  •  o 

ethylene 

509.60 

742.1 

n . 0940 

6.030 

6  1.0 

fthanf 

560 . 30 

709.8 

o . ] 064 

6 .050 

68.0 

PROPyLFNF 

667.20 

667.0 

0.146], 

6.430 

79.0 

DROPANF 

666.90 

..617.4 

0 . 16  3.8 

6  •  4  o  0 

A4.0 

1  R 1 1 T  F  N  F 

7^5.3 

683.0 

O.2o86 

6.760 

96.3 

I-R!'TANF 

784.70 

629.1 

o. 1826 

6.730 

lo5.6 

N-R(  ITANF 

766 .30 

660.7 

0,19^ 

6 .7oo 

101.4 

I -PFNT  ANF 

8  2  0  .  p 

483.0 

0.21 °4 

7.021 

117.4 

n-pfntame 

846.60 

489.5 

0.2387 

7.021 

116.1 

N-HFXANE 

014.20 

440.0 

0.297? 

7.?66 

131 . 6.. 

N-HFPTANF 

972.30 

396.8 

0.3403 

7.430 

147.6 

N-OCTANE 

1024.90 

362.1 

0.399? 

7.551 

163.6 

N-NONANF 

1073 . oo 

332.0 

0.4439 

7  .640 

1  79.6 

N-DFCANF 

1114.00 

304. 0 

0.4869 

7.721 

1  os  .  o 

CHAO  , 

SEADFR  fNtHALPy 

constants 

COMPONENT  6 

AH 

RH 

CH  DH 

H  2  c 

o. 365921 OE 

04 

0 .8033766P 

oi 

o.65o89o4F“o3  o.3?5o9o?F-06 

N ITROGFN 

o.  3192124F 

04 

0 .69564386 

Ol- 

O.3103818E-O4  0. 235O02PF-06 

C02 

o.  3357740F 

04 

0.84401] OF 

01 

0  •  3 0 4 0 5 4 1 6  —  0  ?  -0.78976896-06 

METHANE 

0. 3683894E 

04 

0 . 7880800E 

0 1 

0 .35O6719E-02  ,  ,  0 . 2.4 4.6.4 72E- 0.6. .. 

FTHYLFNF 

o . 379199?E 

04 

0.91 44806F 

01 

0  .78 1  93r4E-o2  -0 . l 6 ?2679E-05 

ETHANE 

0.4220] ?6E 

04 

0  .  ]  120242E 

0? 

0*89381 ?3E-o2  -O.8O66164E-06 

PROPYLFNF 

o.4714762E 

04 

0 . 1350? 1 4F 

01 

0.1 1 46875F-ol  -0. 1 797251 F-06 

propane 

0. 5068633E 

04 

0 . 1510799F 

op 

0 .1540639F-O1  -0. 2O38696F-06 

1  bcjtfnf 

0.5934714E 

04 

0. 1768871F 

01 

0 .1  763632F-ol  -0.3901157E-05 

I -BUTANF 

0.6  op  3.06,7  E 

04 

0. 1999786E 

. 02- 

.0*2.0 441  34F-01  -0.4-1.461  l?F-06  . 

N-Ri  ITANF 

0 . 6654675F 

04 

0 . 206446^^ 

02 

O  .  1  R66?P0F-ol  -O.3O16703F-O6 

i-pentanf 

0. 7475643F 

04 

0 . 24654?  8F 

0? 

0  .?A46?48F-ol  -n . 443734FF-06 

n-pfntane 

o. 8039759E 

04 

0.2518989F 

0? 

o.2348587E-ol  -0 . 4 1 74458E-05 

N-HEXANF 

0 . 94  7 5  o  ?8  E 

04 

0 . 2952267F 

o? 

0 .2914070F-O1  -0.591  1 897F-06 

n-hfptanf 

o. 1056589F 

05 

0 . 3331 526F 

02 

o.3?787 4 6E_ol  -0 .6696 126F-05 

n— oct  a  nf 

0. 1226219E 

05 

0 . 3899065F 

02 

0.3841 118E*~0l  -0.794 

n-nonanf 

o.l 366069F 

05 

0 .437361  8F 

02 

0 ,43ol7?6F-ol  -O.8954763F-06 

n-dfcanf 

o.  1505435F 

05 

0 .4848591 F 

02 

0.475615  8  F~o 1  -0 .9868597F-06 
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CHAO-SFADER  constants  for  thf 

CALCULATION  OF  LIQUID  FUGAClTY  COEFFICIENTS 


GFNFRAL 

mfthanf 

hydROgfn 

A  (  1  ) 

n  .  ^75748nOF 

01 

n.p4^R4^nr>F 

r>  i 

0 . 1  OR71  roof  01 

A  (  2  ) 

761  OIF 

01 

-0. ?24F8nonp 

01 

0 . 1 0207200F  01 

A  (  3  ) 

-o.  40R«snoooF 

01 

-n. ?40R4^n^F 

00 

-'i.SAn'-'pQFQC-n] 

A  (  4  ) 

n .  ?02?0QnnF 

oi 

o.?12nno^^F 

-op 

n.6?ROnnnoF-00 

A  (  5  ) 

o . 

“0.223  00  nr»np 

-07 

0  . 

A  (  6  ) 

0 .84?7n000F- 

oi 

0. 1  048600^8 

n  n 

n  .  p  SR49QQQF-07 

A  (  7  ) 

n  •  26667 OOOF 

on 

-0. 8601 ononp 

n  1 

n  . 

oc 

< 

-n  .  0] 13R000F 

on 

0. 

0 . 

A  (  9  ) 

-0 .76550000F- 

n  1 

0. 

0  . 

Ann) 

0 . pRRononoF- 

01 

0 . 

n  . 

A  (  1 1  ) 

-n  #425807nnF 

oi 

0. 

0 . 

A  (  1  ?  ) 

0 .86580800F 

oi 

0. 

0  . 

A  (13  ) 

-o .  1??06000F 

01 

0. 

0  . 

A  (  14  ) 

-0.01 57040OF 

01 

0. 

0 . 

A  (  15  ) 

-n  .  pFononooP"- 

01 

n. 

n  . 

A  ( 1  6  ) 

n  . 

n. 

0 . 

nitrogen 

C02 

H  2  S 

A  (  1  ) 

0.27365534E  01 

“0. 3no6n874E 

0? 

0.30581 71 nF 

01 

A  (  2  ) 

“0.1 98 1831 OF  01 

0.6 1 409863F 

01 

“0 .7649J006E 

01 

A  (  3  ) 

— 0.51487289F  00 

0.45763220^ 

07 

0 . 87457045F 

on 

A  (  4) 

0.42470 9 88F-01 

“0.27303012F 

07 

-0 . 1464 70 96 F 

n  1 

A  (  5  ) 

-n  .  78 1 4380OF-O2 

0.6O1R7548F 

n  1 

0 .4^7^>47S6F 

no 

A  (  6  ) 

-O.2O474700F-01 

0.36838431 E 

on 

-0  .96721  F7op 

or 

A  (  7  ) 

0.714Q5P48F-01 

“0.6791681 IE 

00 

0 . 147776 4 8E 

01 

A  (  8  ) 

0. 

0. 1 5546365F 

on. 

“0 « 5  Q24773PF 

on 

> 

0 . 

0 . 

O.ORRROnpOF 

on 

A  (  10  ) 

0 . 

0.8956284OF- 

-0 1 

“0.26678483F 

00 
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0.  
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0. 

0. 
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0. 
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PROBLEM  NC.  '3 


INITIAL  PROFILE 


STAGE 

TEMP 

VAPOUR  RATE 

LIQUIC  RATE 

F 

MOLES/U. 

T. 

MULES/U. 

T  . 

1 

8C.CCC 

C.23CCCC0E 

02 

0. 718C600E 

02 

2 

101. SCO 

C.948C6C0E 

02 

0. 7 18C600E 

02 

3 

lll.CCO 

C.9430600E 

02 

0.718C6G0E 

02 

A 

118. SCO 

C.9480600E 

02 

0. 568CoG0E 

02 

126 .CCC 

C . 94806C0E 

02 

0 • 5680600E 

02 

6 

135. 5CC 

C.948C600E 

02 

0  •  5680600E 

02 

7 

143. OCO 

C.948060QE 

02 

0. 56806G0E 

02 

8 

15C.5C0 

C. 603060 CE 

02 

C. 568C6GGE 

02 

9 

174. CCO 

C.603C6C0E 

02 

0. 56806G0E 

02 

IC 

189. 5C0 

C.603C6C0E 

02 

0.568060GE 

02 

1  1 

201. CCO 

0. 6030600E 

02 

0 . 568060GE 

02 

12 

215. 5C0 

C.6030600E 

02 

0 . 56806C0E 

02 

13 

247. CCO 

C.60306C0E 

02 

0.  I223060E 

03 

14 

251. 5CQ 

0.603C600E 

02 

0. 1223060E 

03 

15 

257. CCO 

C.6030600E 

02 

0. 1223060E 

03 

1  6 

264.500 

C  «  60306C0E 

02 

0. 1223060E 

03 

17 

272. CCO 

0. 853C6C0E 

02 

0.  1223060E 

03 

18 

28C.5C0 

0. 853060GE 

02 

0. 1223060E 

03 

19 

292. CCO 

0. 853C600E 

02 

0.12230606 

03 

2C  , 

311. 5C0 

0. 8530600E 

02 

C.  1223060E 

03 

21 

349. CCO 

0. 853C6C0E 

02 

0. 370000GE 

02 

FINAL  PROFILE 


STAGE  TEMP  VAPOUR  RATE  LIQUID  RATE 


F 

MCLES/U. 

I. 

MOLES/U. 

T. 

1 

8C.569 

0 . 23C0C00E 

02 

0.7 180600E 

02 

2 

101.647 

C.948C600E 

02 

0.72 1244  IE 

02 

•2 

UC. 772 

C. 9512 44  IE 

02 

0.7119993E 

02 

4 

117.774 

C.9419994E 

02 

0 . 5430282E 

02 

C 

126.175 

0 . 923C283E 

02 

0.521835 86 

02 

6 

134.862 

C.9018358E 

02 

0 • 5024468E 

02 

7 

143. 24C 

C. 8824468E 

02 

0.4871 02 7E 

02 

£ 

166.605 

0*522102 7E 

02 

0.4880472E 

02 

9 

181. 51C 

C. 523047 3 £ 

02 

0.4853166E 

02 

1C 

192. 7C3 

C. 5203 1*6  7E 

02 

0 .47920886 

02 

1  1 

202.819 

C.5142089E 

02 

0.4623972c 

02 

12 

216.396 

C.4973973E 

02 

0.4099080E 

02 

12 

248.359 

C  .44490806 

02 

0. 10339596 

03 

14 

252.590 

C.4139586E 

02 

0. 1C46293E 

03 

15 

257.814 

C.4262927E 

02 

C.  105S403E 

03 

16 

■264.856 

C.4384831E 

02 

0. 1070336E 

03 

17 

271.837 

C.  7C0336CE 

02 

0. 1075631E 

03 

18 

28C.165 

C.7056309E 

02 

C.  I073619E 

03 

19 

291.853 

C. 7C36186E 

02 

0. 1059080c 

03 

2  C 

31C.813 

C . 689C8C0E 

02 

0  •  1008616E 

03 

21 

348.489 

C.6386163E 

02 

C . 3700000E 

02 

SUM 

(DELTA  T**2) 

AT  ITERATION  7 

=  C.840 

. 


oooooes.o 

at  dl 

3  3C8SOS.?,0 _ .c _ 
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PRODUCT  STREAM  COMPOSITIONS 


COMPONENT 

CONDENSER 

BOTTOM 

LIQUID 

VAPOUR 

RECYCLE 

PRODUCT 

SIDESTREAM 

SIDESTREAM 

METHANE 

0. 1041E-01 

0.2305E-15 

0.2529E-  2 

0. 73  70£- 10 

ETHANE 

0.2I07E  CO 

0.2; r  6E-09 

0.634  7  E-0  1 

0. 1866E-05 

PROPYLENE 

0.2262E  CO 

0.61 ' Ot-04 

0.1934E  00 

0.9834E-02 

PROPANE 

0.5349E  00 

C.9-.  74E— 03 

0.5915E  00 

0.85 9 IE- 01 

I -BUTANE 

0. 1 104E-01 

0. 1849E-01 

0.654 1 E— 0 i 

0 .  1 302  E  00 

N -BUTANE 

0 . 6747E-02 

0.1402E  00 

0.831 6E-C I 

0.5407E  00 

N- PENTANE 

0 . 3 147E-05 

0.30650  00 

0.4771E-03 

0.  1 732 E  00 

N- HEXANE 

0. 1726 E- 10 

0.2142E  3 

0. 3353E-07 

0.422  8  E-0 1 

HEPTANE 

0.2139E-14 

0. 1140E 

0. 4105E-10 

0. 1095E-01 

OCTANE 

0. 2753 E- IB 

0.1122E  C 

0.6351E-13 

0. 5538E-02 

OTHERS  (360) 

0*28 16E-25 

0. 9340E-01 

0.5109E-1B 

0. 1462E-G2 

■ 


■V 
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FINAL  PROFILES  AND  COMPOSITION  AS  REPORTED  IN  LITERATURE  (REF.  26) 


STAGE  TFMP 


1 

2 

4 

5 

6 

7 

8 
Q 

10 

UL 

12 

13 

14 

15 

16 


18 

19 

20 
21 


r 

77.290 

94.770 
103.610 
111.070 
120.780 
131.270 
141 .410 
169.450 
186.400 
198.520 
_2j09 .140 
222.910 
253.420 
257.560 
262.51 0 
269.050 
275 .440 
282.950 
293.300 
309.880 
343.330 


VAPOUR  RATE 
MOLES/U.T. 

0 • 2300000E  02 
0.9480600E  02 
0.95342 OOE  02 
0.9461700E  02 
0.9271300E  02. 
0.9055700E  02 
0.88691 OOE  02 
0.52^6900E  02 
0.5386300E  02 
0.5424200E  02 
0 .53827 OOE  02_ 
0.5 1 99700E  02 
0.4659000E  02 
0 . 4447900E  02 
0.4580Q00E  02 
0.4681400E  02 
0.726 5 2 OOE  02 
07265  8  OOOE  02 
077180500E  02 
076977500E  02 
064589000E  02 


i 


COMPONENT 


LIQUID 

SIDESTREAM 


vapour 

SIDESTREAM 


methane 

FTHANE 

PROPYLENE 


PROPANE 
I -BUTANE 


N-BUTANE 

n— pentane 

HEXANE 


HEPTANE 


OCTANE 
OTHERS ( 360 ) 


0.26010E-01 
0.88760E  00 
0.2977 3 E  01 
0.88691E  01 
0.99410E  00 
0.12404E  01 
0  •  548  OOE- 02 
0.41940E-06 
0 . 3  4  5 7 OF— 09 
0.24280E-12 
0. 16310E-18 


0. 19480E-10 
0.20810E-04 
0.21  48 0 E  00 
0.19773E  01 
0.32555E  01 
0.13705E  02 
0.42934E  01 
0.11099E  01 
0.28010E  00 
0.13390E  00 
0.29660E-01 


MA3ATS3012 


. 


lo  3990X1,0 

no  aoro3£.r 


. 
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PROBLEM  NO.  4 

starilizer  problem 

PROBLEM  STATEMENT  AN Q  DATA 


A.  COLUMN  DATA 

PR PS. BUR r=  3  4  5.0  PS  I  A 

NO.  OF  PLATES*  17 

rfflux  rati 

0=  ”3.  1.40Q 

distillate 

RATE=  0.44423 

MOL/IJ.T. 

FEED  PLATE 

NO  •  =  9 

*  T 

R.  FPFD  DATA 

.  } 

FLASH  TEMPERATURE*  321.83 

F 

<  i 

FrFn  LIQUID 

FRACTIONS  0.5550 

1 

FEED  R AT E  = 

1.00  MOL/IJ.T. 

» 

COMPOS  I T ION 

IN  MOLE  FRACTION 

) 

component 

TOTAL  FEED 

FLASH  LIQUID 

FLASH  VAPOUR 

N I trogfn 

0. 5426F-03 

0. 8  lORoF-^4 

3.111  8  2  F  —  0  2 

HAS 

o • 1 1 5 oE~oo 

0  •  4  9  6  9  9  E  —  O  1. 

3. 1 9651F-00 

CO  2 

0. 1 194F-01 

0.331  1 9  F  —  0  2 

0.2269 4 F  —  01 

methane 

0. 8410E-01 

3. 16673F-oi 

3 . 1 6823F-30 

ETHAN F 

o .4069F-0 I 

3. 1 3322E-3 1 

0.74834F-O1 

propane 

0. 7379E-M 

0.42773E-A1 

0.1 1248 E -00 

I  —  R i  ITANE 

0. 2  9  3  OF -0 1 

0.227 9 2F  —  31 

0.3741  6  F  —  0 1 

M-Bi  ITANF 

o. 9333P-01 

3.77259E-31 

o.l 1 3  36F  — no 

i -pentane 

0. 5209F-01 

0.53857E-31 

O.49883E-01 

n -pentane 

s.6131 F - o 1 

0.655S3R-O1 

0. 56023F-O1 

HEX ANF 

0. 1 145F-00 

0 . 1467  5E-30 

0.74246E-O1 

HEPTANE 

0.1 264^—0 1 

0. 1.  8  2  74E-3  0 

0.561 77F-01 

OCTANE 

0.6782F-01 

0.13663E-^3 

3.1 9449F-01 

NONANE 

0.3961E-01 

0.65587E-0 I 

0  •  720  8  2  E— 3  2 

DEC AN E 

0. 8953E-01 

0*1 5  298E-00 

3.1 03P3E-31 

C.  EQUILIBRIUM  AND  ENTHALPY 

DATA  ARE  GIVEN 

BY  THF  CHAO- 

S F A o E R  CORRELATION  AND  A.Rf 

reported  for  prorlfm  3 

•*  M  f 


f  -  ,  ,  o 


^-BSA<vgT.f> 


.  fl  I 


PROBLEM  NO*  4 
INITIAL  PROFILE 


STAGE 

TEMP 

VAPOUR  RATE 

LIQUID  RATE 

F 

MOLES/U. 

T. 

MOLES/U. 

T. 

1 

140.000 

0. 4440000E 

00 

0. 1390000E 

01 

2 

170. COO 

0. 1835000E 

01 

0. 1390000E 

01 

3 

190.000 

0.  1835000E 

01 

0. 1390000E 

01 

4 

210.000 

0. 1835000E 

01 

0. 1390000E 

01 

5 

230.000 

0. 1835000E 

01 

0. 1390000E 

01 

6 

250.000 

0. 1835000E 

01 

0. 1390000E 

01 

7 

27C.000 

0. 1835000E 

01 

0. 1390000E 

01 

a 

270.000 

0. 1835000E 

01 

0.  1390000E 

01 

9 

290.000 

0. 1835000E 

01 

0. 1945000E 

01 

10 

33C.0C0 

0. 139000CE 

01 

0. 1945000E 

01 

u 

350. COO 

0. 1 390000E 

01 

0. 1945000E 

01 

12 

370.000 

0. 1 390000E 

01 

0.  1945000E 

01 

13 

390.000 

0. 1 390000E 

01 

0. 1945000E 

01 

14 

410.000 

0.  1 390000E 

01 

0. 1945000E 

01 

15 

430.000 

0. 1390000E 

01 

0. 1945000E 

01 

16 

450.000 

0. 1 39C000E 

01 

0. 1945000E 

01 

17 

470.000 

0. 1390000E 

01 

0.5550000E 

00 

3T4*  GIUOU 
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FINAL  PROFILE 


STAGE 

TEMP 

VAPOUR  RATE 

LIQUID  RAIE 

F 

MULE  S/-J. 

T. 

XQLES/U.T 

• 

1 

144.998 

0.444 2 300E 

00 

0  .  1  3 9  4  800E 

01 

2 

20C  .86  1 

0.  1839030E 

01 

0. 153 54  6  0 E 

01 

3 

221.263 

0. 1 979690E 

01 

0. 15  7  7217E 

0  1 

A 

231.855 

0.202  144/E 

01 

0. 1  5  6  8  7  3  6  E 

01 

5 

24C .  143 

0. 201  296t>E 

01 

0. 1533695E 

01 

6 

248.973 

0.197  / 9 2 5 E 

01 

0 . 1  4  /  3  3  1  9  E 

01 

7 

260.350 

0.191 7549E 

01 

0 . 1366  1 4 2 E 

0  1 

8 

277.991 

0. 18 10 372 E 

0  l 

0. 1  145  7 9 8 E 

01 

9 

313.363 

0.  1  590028E 

01 

0. 1698264E 

01 

10 

34 C. 893 

0.  1  142494E 

01 

0. 1923782E 

01 

1  1 

355.696 

0. 1368012E 

01 

0.202 59 38E 

0  l 

12 

367.  104 

0.  1470168E 

01 

0.20781 9 9E 

01 

13 

378 . 298 

0.  1522429E 

01 

0.2103 3 06E 

01 

14 

391 .228 

0.  1 5 4  7 5 3 6 E 

01 

0.2 10491 OE 

01 

15 

408.094 

0.  1 549140E 

01 

0 .20 7  44  3 3E 

0  1 

16 

432.408 

0. 1518663E 

01 

0. 200664 3E 

0  1 

17 

469 . 146 

0. 1450873E 

01 

0. 555770GE 

GO 

SUM 

(CELTA  T**2) 

AT  ITERATION  17 

__ 

0.938 

. 

PRODUCT  STREAM  COMPOSITIONS 


COMPONENT 

NITROGEN 

CONDENSER  STREAM 

0.1 226082E-02 

BOTTOM  STREAM 

0.  1  1  59296E- 10 

H2S 

0.260119 A E  00 

0. 51 53156E-05 

C02 

0.2699206E-C1 

0. 4390543E-08 

METHANE 

0.19010  34E  00 

0. 682  37  84E-08 

ETHANE 

0.92023  15E-01 

0. 5990067E-07 

PROPANE 

0.1668011E  00 

0. 7923181t-04 

I  - BUTANE 

N  -  BUTANE 

U.6434291E-01 

0.  153474  IE-02 

0.1 9243  7 2E  00 

0.  1481454E-01 

I  PENTANE 

0.4459652E-0 2 

0. 09850 36E-01 

N  -PENTANE 

0.1 526564E-02 

0.1087113E  00 

HEXANE 

0.6  0316  76E-05 

0.2052455E  00 

HEPTANE 

OCTANE 

0.35  7  73  92E-07 

0. 2266508E  00 

0.84262 1 5E-10 

0.1215938E  00 

NONANE 

0.30244 4 4E-12 

0. 7101077E-0 l 

DECANE 


0.5967150E-14 


0.1605038E  00 


:1 
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LITERATURE  VALUES  (REF.  39) 


THIELE  GEDDES  CALCULATION  PROCEEDURE 
BASIS  ONE  LB-MOLE  FEED 
COLUMN  PRESSURE=  345.000  PSIA 


DIS 

TILLATE  TEMP=  146.123 

NO  OF  THEORETICAL  PLATES=  16 

REFLUX  ra 

T 10=  3.14091  DISTILLATE 

RAT  E=  0.44423 

FEED  PLATE  N0=  8 

stage 

VAPOUR 

LIQUID 

TEMPERATURE 

1 

i. 

1.83951 

1.53474 

202.1 1 1 

z 

1.97897 

1.57328 

222.727 

-> 

o 

2.01751 

1.56150 

233. 726 

2.00573 

1.52339 

242.574 

5 

1.96762 

1.45883 

252.103 

6 

1.90306 

1.34198 

264.407 

7 

1.78621 

1.65781 

283.787 

8 

1.65704 

1.90883 

336.363 

9 

1.35306 

2.02336 

351.323 

10 

1.46759 

2.08364 

361.942 

i  I 

1.52787 

2.11882 

371.555 

12 

1.56305 

2.  13529 

381.768 

13 

1.57952 

2. 13252 

394.058 

14 

1.57675 

2. 09780 

410.391 

15 

1.54203 

2.02750 

434.226 

16 

1.47173 

0.55577 

470.361 

% 

CONDENSER 

HEAT  L 0AD=  0.10755995E  05  BTU/LB  MOLE 

* 

FEED 

REBOILER 

HEAT  LQAD=  0. 14896920E  05  BTU/LB  MOLE 

FEED 

COMPOSITIONS 

IN  MOLE  FRACTIONS 

V 

BMP ONE NT 

FEED 

DISTILLATE 

BOTTOMS 

NilTROGEN 

0 . 54259359 E— 03 

0.  122142  22E-02 

0.20757247E-L1 

H2S 

0.11502984c  00 

0. 25893895E  OC 

0.20417795E-05 

C02 

0. 1 1937059E— 0  1 

0. 26871287E-0 1 

0. 11  195078E-08 

'ETHANE 

0. 84102007 E—  01 

0.  1 8932044E  00 

0. 14462437E-08 

ETHANE 

0. 40694519 E—  01 

0. 91606643E-01 

0. 17935442E-07 

’ROPANE 

0 . 7  3  7  9  2  7  2  3  E—  0  1 

0.  16605930E  00 

0.43258454E-04 

-BUTANE 

0 . 29300054E-0 i 

0. 64459787E-01 

0. 1 1965750E-02 

-BUTANE 

0 . 93 3 2 6097 E— 01 

0.  1  94492  lie  00 

0. 12463225E-01 

-PENTANE 

0 . 520'8  8985E-01 

0.  5  18774 4 8 E- 02 

0.8957755 7 E-01 

-PENTANE 

0. 61313076 E- 01 

0-  1834 1 58  8 E— 02 

0.  1 0885509E  00 

EXANE 

0 . 1  1 4 4 8 7 2 5 E  00 

0. 81200844E-05 

0. 20599 1 4 1 E  00 

EPTANE 

0.1 2 84243 IE  GO 

0.  5 3 3 4 6 4 5 7 E - 0 7 

0.22747631E  00 

w  .  /  \  •  .  .3 

•  :  */  •  j  /  1  C»  .  .1  __  f 

W  v  O  «  •» .J  i  .  >  •  7  u  ..  U  . 

0.  1  24 7 964  6 E— 0  9 

0. 1220366 7 E  00 

*  J .  A  i\  u 

0  .  ^9609^-32 c—  0  1 

0.53066324E-12 

0.7  126941 7E-0 1 

-CANE 

0  c  3  9  5  2  7  9  4  2  E—  0  1 

0. 1 1019513E-13 

0. 161088 4 IE  00  * 

3TA'jrjX3lT3«03HT  30  0^ 
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PROBLEM 


APPENDIX  IV 

STATEMENTS ,  INITIAL  AND  FINAL  PROFILES, 


AND  PRODUCT  COMPOSITIONS  FOR  PROBLEMS  V  AND  VI 


VI  xiana«m 


PROBLEM  MO.  5 

fxtRactivf  DI ft t  lla  t  t ON 
6Y6TEM  PHONOL-MCH-TOLt  IFNF 


a.  COLUMN  0  A  T  a 

PRESSURE-ATMOSPHER I c 
NO.  OF  PLATFS=21 

CONDFNSFR  TYPF=TOTAL _ 

OVFRHFAD  PR0DUCT=62 .oo  MOL/ll.T. 
BOTTOM  PRODUCT5  247  .PI  MOL/ll.T. 
RFFLl'X  RAT  10  =  5  •  n 
FFED  PLATES*  5*14 


B.  FFFD  DATA 

TFMPFR AT!  IRF  of  SOLVENT  FFFD=244.49  F 
RATF  OF  SOLVFNT  FFFD  =  200. 0  MOL/JJ.t. 
TEMPERATURE  OF  FRFSH  FFFD*  241.71  F 
RATF  OF  FRFSH  FFFP=  100.0  MOL/U.T. 


COMPONFNT  COMPOSITION  COMPOSITION 
SOLVENT  FEED  FREfh  FFFD 
PHENOL  0*9910  0*0000 

MCH  O*0f\nn  o#50OO 

TOLUFNE  n.^oQo  o . raoo 


C.  EQUILIBRIUM  DATA  ARE  OBTAINFD  FROM  THE  SIMPLIFIED 
FRAUSNITZ  MODEL  USING  WILSONS  equation* 

THE  OPTIMIZED  INTERACTION  ENERGIES  ARE  IN  CAL /GM. MOLE 
WHERE  N  =  PHENOL  M  =  MCH  T  =  TOLlJFNF 


I NTPR ACT  I  ON  FNFRGY(N,L)=  155*2845 
INTERACTION  FNFRGY(L,N)=  967.799 
INTERACTION  ENERGY ( N  * M ) =  60.278 

INTERACTION  ENERGY (M *N ) *  1875.423 
INTFRACTI ON  FNFRGY(M,L)=  278.607 
INTERACTION  FNFRGY ( L ,M ) =-  72.162 


D.  ENTHALPY  DATA  ARE  OBTAINED  FROM  IDEAL  ENTHALPY 
FUNCTIONS  WITH  HEAT  OF  SOLUTION  CORRECTIONS 
THE  INFINITE  DILUTION  HEATS  OF  SOLUTION  aRf  IN  BTU/Lr-MOL 


HF AT  OF  DIL.(N*M)  =  HEAT  OF  DIL.(M*N>=  4627.7 
HEAT  OF  DIL.(N*L)  =  HF  AT  OF  DIL.(L,N)=  27Q7.7 
hf A t  OF  PIL(M*L  )  =  HEAT  OF  n  I  L .  (  L  , M )  =  7.  nor 
THESE  VALUES  ARF  OBTAINED  FROM  SMITH (6) 
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initial  profile 


STAGE 

TFMP 

VAPOUR  1 

RATF 

LIQUID  RATF 

F 

MOLF  c-/\ « 

•  T. 

MOLF 6 /U 

•  T. 

1 

216.18 

0. 

0.26046F 

03 

2 

216.64 

0.31264F 

03 

0 .26761F 

03 

3 ' 

216,4? 

°.7O080F 

03 

0 .2  6407 F 

O  3 

4 

220.^0 

0. 30616F 

03 

0 .249100 

03 

5 

240.14 

0.3011  OF 

03 

O ,4?o96F 

03 

6 

2  4  r .61 

o.32^0KF 

03 

0.41 869F 

0  3 

7 

240.89 

O.3207PF 

03 

0.41 654F 

0  3 

8 

241 .26 

0.31863F 

03 

0 .41461E 

03 

9- 

241.6? 

•  3  1  6  6  0  F 

03 

0.41 2  6 1 F 

03 

in 

241.96 

0.31470 E 

03 

0.41084F 

03 

11 

242.26 

0.31793F 

03 

0 .40922F 

03 

12 

242.4? 

0.311 3  1  F 

03 

0 .40784F 

03 

13 

241.89 

0.309036 

03 

0 .40732F 

03 

14 

237.3 2 

O.30O41F 

03 

0 . 6o3ooF 

o3 

1  5 

237.66 

0. 3o6o9F 

03 

0 . 6  O094F 

03 

16 

238.08 

0.303O3F 

03 

0 .49832F 

03 

17 

238.64 

O.30O41F 

03 

0 .49606F 

03 

18 

2^9.43 

0. ?Q7 14F 

03 

O.49101F 

03 

IP 

240.80 

n.  793  1  OF 

03 

0 .48601F 

03 

20 

246.37 

0.288  1  OF 

03 

0 .48077F 

03 

1  TRT*  3  3  0*28?86E  05  07T9T9IF  03 


final  profile 


STAGE 

tfmp 

VAPOUR  1 

RATF 

LIQUID  RATF 

F 

MOLF  9 /!J 

.T. 

MOLES /i.  I. 

T. 

1 

216.12 

0. 

0.26  o46F 

03 

2 

216.48 

0.31264F 

03 

0 .25759F 

03 

~T~ 

216.49 

o.3no68Tr" 

03 

0.2641  IF 

03 

4 

221 .4? 

0 . 3  06  2  OF 

03 

O.24906F 

03 

6 

244.49 

0.301 1 5F 

03 

0 .47348F 

03 

6 

244.86 

0.32667F 

0  3 

0.471 1  IF 

o5~ 

7 

245.2? 

0.32320F 

03 

0 . 46  8  84F 

03 

8 

245.60 

0 . 32093F 

03 

0 . 46669" 

03 

0 

245.97 

o.3l 6  78  F 

o3 

0 .46  466F 

03 

1  0 

246.33 

0.31 676F 

03 

0.46278F 

03 

11 

246.66 

0.31487F 

03 

0.461 04F 

03 

12 

246.85 

0.513 13 F 

03 

o.45956F 

0  3 

13 

246.30 

0 . 3 1 1 6  5  F 

03 

0 .45898F 

03 

14 

241.31 

0.311 07F 

03 

0 .55405F 

03 

16 

241.63 

0.30614F 

03 

0 . 6  5  2  0  3  E 

03 

16 

242.06 

O.30412F 

03 

0 . 54947F 

03 

17 

242.62 

0.301 56F 

03 

0 .646266 

09 

18 

243.40 

0.29836E 

03 

O  .642316 

03 

19 

244.77 

0. ?9440F 

03 

0 .537426 

03 

20 

249.45 

0.28961F 

03 

0.53282^ 

03 

2  1 

2  /  6 . 6  1 

“TT.  284  9  1  F 

03 

0 .247916 

03 

1  1 1  JA  1-  Jtf  f 
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product  strfam  composition*; 


COMPONFNT  CONDFNSFR  STRFAM  B0TT04  S T R F 1 M 


PHFNOL 

n  .onm 

n . 7097 

MCH 

0.785? 

n.  076 ft 

Tf/turNF 

'  C • 2 1 45 

n  #  7  f>  0 

rfcovfRv 

■or  rotr  fmf1 

"ACM  T  FVFTj  =  n  .  7  /  6  6 

p-O  V 
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PRORLFM  NO.  6 . 

EXTRACT  I VF  DISTILLATION 
system  ACF  TONF-METHANOL-watFK 


A. 

COLUMN  DATA 

P  RF  s SI  iR  F  -  A  T MO  S PHER  I 

NO.  OF  dLAtfs=25 
condfnsfr  typf=total 

OVFRHFAD  PRODUCT  =os • 09 

BOTTOM  PRODUCT  =575.34 

RFFLUX  RAT  10  =  4 • 1 3 

FFFD  PLATES*  0*20 

B. 

FFFD  DATA 

TFMPFRATURE  OF  SOLVFNT  FFFO=l4i .307  F 

RATE  OF  SOLVENT  FFFD=  115.56  MOL/U.T. 

TFMPFRATURE  OF  FRFSH  FFFD=  167.71  F 

R ATF  OF  FRFSH  FEED*  402.87  MOL/U.T. 

COMPOSITION  IN  MOLf  FRACTION 

COMPONENT  COMPOSITION  COMPOSITION 

SOLVENT  EFFD  FRESH  FFFD 

ACETONE  0.0000  0.0718 

MFTHANOL  o.oooo  o.n972 

WATFR_ l.oooo  0  .  8310 

C. 

EGlI  I  LIBRIUM  DATA  ARE  OBTAINED  FROM  THE  SIMPLIFIED 

PRAUSNITZ  MODEL  USING  WILSONS  EQUATION. 

THE  OPTIMIZED  INTERACTION  ENERGIES  ARF  IN  CAL/GM. 
WHERF  N=ACETONF  M=MFTHANOL  L=WATFR 

MOLE 

INTERACTION  ENERGY (N  ,L ) =-164.870 

INTERACTION  ENFRGY(L,N)=  1679.652 

I NT FR ACT  ION  ENERGY ( N  ,M ) =-31 5 .893 

INTERACTION  ENERGY ( M  *N  )  =  1710.198 

INTERACTION  ENERGY ( M  *  L ) =  241.057 

INTERACTION  FNFRGY(L,M)=  -*54.830 

D. 

ENTHALPY  DATA  ARE  OBTAINED  FROM  IDEAL  ENTHALPY 

FUNCTIONS  WITH  HEAT  OF  SOLUTION  CORRECTIONS 

THE  INFINITE  DILUTION  HEATS  OF  SOLUTION  ARE  IN  BT 
HF  A  T  OF  DIL.(N,M)  =  HEAT  OF  DIL.(M,N>=  o.nono 

U/Lr-MOL 

HF AT  OF  DIL.(N*L)  =  HEAT  OF  DTL.(L,N>=  4520.0 

HEAT  OF  DIL(M*L)  =  HF  A  T  OF  DIL.(L,M)=  0600. n 

THESF  VALUES  ARE  OBTAINED  FROM  INTERNATIONAL  CRITICAL 

TARLFS ( ^2 ) 


TJ 


.  <  R' 

*  1 1  r  •a  f  -  ■>  < '  0  T  {0  A  v’  ^  T 


^  _  •  ;T  "  K,  '  . 

- 

. 

* 

nr  ,f  s?-*TA.J'  rn  ^  1 


l'  T  A  V  * 

,?!f  * 

3  TV -Yd  r  sOTR  ; '? 

.  .  ’ 

•  tV'  ’fS^v  ■  J-  :  ^ i  1  •  ttI r'oq-.oo 

T 

- 

. 


nf  r» 


rv  A  A  r 


a aT  *>.< 


1T  •,»  q.ip  3HT  '  ,)'-  !  C  WIAT  10  AT*0  I .  I  V  >\ -i  I  < 

n^,oi  --  MO?  J  T  fW  r  iM  tTI/PUAW 


\JAO 


m 

IT  '* 


=  J 


IOVRIR  kiO!  T >1  1T^  I  R5T  IT^11 


Jj  i  Ml  ,k  = 

.  r  q  .  [  -  =  (  J 

. 

. 

r .  r  w  =  ( *  i 

r  •  •  f  A  r 


'  M  0  A  Ml  T  /•  I 
)  v  H  ':  •  OTTO*'  OT -U 

>  yom  v  m  utoambt^i 


•  —  \  ' 
-  ■  ■  . 


v  .1  4  J  -  u  I  KMH  WUT80  m-  ta  1  *;^T"T  ■  ., 
,T  •  .-I i  1  MUItllJU?  ^0  TA  ■  I  MT!  •  ?T>0IT1l  1 

...  '  • 

■  •  '  ' 

.  ,  ~ 

' 


T  1  MO 


PROBLEM  NC.  6 


INITIAL  PROFILE 


STAGE  TEMP 

F 


VAPOUR  RATE 
MCLES/U. T. 


1 

122. OCO 

C. 

2 

122. CCO 

0.  185CCGUE 

03 

2 

1 2  2  .  C  C  0 

C.  1  85CGCCE 

0  3 

4 

122. CCO 

0. 185C00UE 

03 

1 2  2 . 0  C  0 

C. 185CCCGE 

03 

6 

122. CCO 

C. 185CC00E 

03 

7 

14C.CC0 

C. 1 8CCCCCE 

03 

E 

14C.CC0 

C. 18CCC00E 

03 

5 

14C.0CC 

0. 1 8CCCC0E 

03 

IC 

14C.CC0 

0. 1 8CCC0GE 

03 

1  1 

14C . C  C  C 

0. 180CCC0E 

03 

12 

14C.CC0 

0. 180CCCCE 

03 

13 

14C.0C0 

C. 180CCC0E 

03 

14 

14C .CCO 

0. 18CCCCGE 

03 

15 

158. CCO 

0. 18CCC00E 

03 

16 

158. CCO 

0. 18C0CCOE 

03 

17 

158. OCO 

0. 18CCCCCE 

03 

ie 

158. CCO 

C. 18CCC00E 

G3 

15 

158. CCO 

C. 18CCC00E 

03 

2C  % 

176. CCO 

0. 180CC00E 

03 

2  1 

176. OCO 

0. 18CCOCOE 

03 

2  2 

176. COO 

C. iSCCOOOE 

03 

22 

176. CCO 

C. 1 800C0GE 

03 

2  4, 

176. CCC 

o. lacooooE 

03 

25 

176. OCO 

0. 1 8CCC00E 

03 

LIQUID  RATE 
MGLES/U. T . 

0  . 1 5200GGE  03 
0.140000GE  03 
0  •  140000GE  03 
0. 14000000  03 
C.  14C0CGGE  03 
0.140000CE  03 
0.1400000E  03 
0 . 25CCG00E  03 
0. 250C0G0E  03 
0.2500000c  03 
0 • 250C000E  03 
0 . 250C0C0E  03 
0.25000 COE  03 
0.25G00G0E  03 
0  •  25000G0E  03 
0 •  25000G0E  03 
0.2500000E  03 
0. 2500000E  03 
0 . 2500000 E  03 
0.7500000b  03 
C . 7500000E  03 
0. 7500000E  03 
0.7 500000E  03 
0 . 7500CGGE  03 
0.7500000E  03 
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FINAL  PkCFILE 


STAGE 

TEMP 

VAPOUR  RATE 

LIQUID  KATE 

F 

MOLES/U. 

T. 

MCLES/U. 

T. 

1 

133.185 

C. 

0. 1366617E 

03 

2 

133.237 

C.  169751 7E 

03 

0. 1364129E 

03 

■2 

133.298 

0.  1  695029E 

03 

0. 1361092E 

03 

4 

133.379 

0.16919 9 2E 

03 

0. 1357152t 

03 

C 

133.499 

0.  1 688052E 

03 

0. 1351557E 

03 

t 

133. 7C0 

C.  16824  5  7E 

03 

0. 1 342344E 

03 

7 

134. 1C9 

0. 1673244E 

03 

0. 1321399E 

03 

e 

135.247 

C.  1652299E 

03 

0.1 19288  1 E 

03 

9 

141.305 

C.  152378  IE 

03 

0 . 2484  692  £ 

03 

IC 

141.441 

0. 1659992E 

03 

0.2481724E 

03 

i  1 

141.612 

C. 1657024E 

03 

0.247  8304E 

03 

12 

141.84  7 

0. 16536C4E 

03 

0.2473930E 

03 

13 

142.195 

0. 1649230E 

03 

0 -  246  7 1 2 1 E 

03 

14 

142.765 

C.  1642421E- 

03 

0.2454726E 

0  3 

15 

143.787 

C. I63C026E 

03 

0  •  2432269E 

03 

16 

145.686 

0. 1 607569E 

03 

0.2397327E 

03 

17 

148.990 

0. 1 57 2 62 7 E 

03 

0.2353199E 

03 

16 

153.879 

C. 1528499E 

03 

0 . 2304975E 

03 

19  . 

16C.C29 

0. 1 48 02 75 E 

03 

0.225348 3 E 

03 

2  C 

167.717 

C. 1428783E 

03 

0. 7354395E 

03 

2  1 

173.355 

G. 1600996E 

03 

0 . 7353034E 

03 

2  2 

178.462 

0. 1599634E 

03 

0. 7355438E 

03 

23 

182.617 

0. 1602038E 

03 

0.7350413E 

03 

24 

187. 2C5 

0. 1597013E 

03 

0.7330283E 

03 

25 

194.310 

0. 1576884E 

03 

0.575  3400E 

03 

SUM  (DELTA  1**2)  AT  ITERATION  34 


0.  539 


. 
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PRCCOCT  STREAK  COMPOSITIONS 


CCK  PCX  ENT 
ACETONE 
METHANOL 
WATER 


CONDENSER  STREAM  BOTTOM  STREAM 

0.9637215E  00  0. 

0.5497G82E-02  0. B34608CC-0  1 

C.30781A2E-01  0.9165392E  00 


RECOVERY  GF  METHANCL  ACHIEVED 


1.C0232 
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APPENDIX  V 

EXAMPLE  COMPUTER  PROGRAMS 
-  PROBLEM  III 


PROBLEM  VI 


SIBFTC  THESIS 

DIMENSION  T ( 25 ) >V( ?5 ) »PROFL (25>*G(25)*U(25)»R(?5> 
DIMENSION  VK(15>?5)»F(25)>C(?5),D(?5*15)*P(25)*Q(?5) 
DIMENSION  X ( 15*25) *SUMX(25>  *AT(25)*Y(]5*25)  ,BBMIX(25> 
DIMENSION  ZC( 25 > *AAOR (25) *A ( 17*5 ) *TC( 15 ) *PC ( 15 ) »W( 15 ) 

DIMENSION  DELTA( 15  )  *V1  (  15  )  >DDELMX (25)  *SUMY (25 > _ 

oa  FORMAT ( 1 HI  *  44 X  *  1 3HPR0BL EM  NO.  3) 

99  FORMAT (1X*I2*F5«1*F6*3*F7*3) 

100  FORMAT ( 1HJ»43X»15HINITIAL  PROFILE) 

101  FORMAT ( 1HK*20X*5HSTAGE*6X*4HTEMP*10X*11HVAPOUR  RATE*RX 
1 ♦11HLIQUID  RATE) 

_ L02  FORMAT  (  1X»32Xj>  1HF  *  1  2X  *  10HMOLES/U  .  T  .  ,  l  OX  ,  10HMOLE  S/U  .  T  . L ) _ 

103  FORMAT ( IX *20X *  1 2 *6X ,F8 .3 *2 ( 6X *E14. 7 > ) 

106  FORMAT ( 1 HK  *44X  *  1 3HF IN AL  PROFILE) 

104  F0RMAT(1HK*21X*29HSUM  (DELTA  T**2)  AT  I TERAT I  ON , IX  *  I  2  * 

11X * 1H=,F6.3 ) 

105  FORMAT ( 1X*2F5« 1  ) 

.  110  FORMAT  (  1  X  >_EX4-«uB_l _ 

111  FORMAT ( 1H1 *34X »27HPRODUCT  STREAM  COMPOSITIONS) 

112  FORMAT ( 1HK , 1 9X  *  9HC0MP0NENT *  2X  * 9 HCON DENSER  *6X  *  6H BOTTOM ♦ 

16X  » 6HVAP0UR  *  6X  * 6HL  IQ(J  I  D  ) 

115  FORMAT ( IX *32X *7HPR0DUCT  » 6X * 7HPR0DUCT  *4X *  1 OHS  I  DESTREAM ♦ 3X * 
llOHSIDFSTREAM  ) 

_ 113  FO RMAT( lHJ»29X»4(E11.4t2X> ) _ 

NP  =  2  1 
NCP=1 1 
NPPaNPrl 

Read (5  *99)  (j*T(j)*v(j) *profl ( j ) , j=i  ,np  > 

sw=o.o 

_ W R  I.TE  (  6  »  98  ) _ 

WRITE (6  *100) 

WRITE (6 *101 ) 

WRITE  (6.*  102)  _ 

DO  901  J  =  1 *  N  P 

901  WRITE (6  *  103)  J * T ( J ) ♦ V ( J )  * PROFL ( J ) 

_ DO  13  I  =_1_*  NCP _ 

DO  13  J=1*NP 
13  D ( J  *  I ) =0  •  0 

C  _ 

C  JI  IS  FEED  LOCATION.  D  IS  FEED  COMPOSITION  VECTOR  IN  MOLE 
C  PER  UNIT  TIME 
£ 

J I  =7 
N  J  =  1 3 

READ(5>105)  (  (D( J  I  *  I  )  *D(NJ*I )  )  *I=1*NCP> _ 

CALL  FPTSW 

92  IF (SW.EQ.O.O)  GO  TO  77 

_ GO  TO  123 _ ‘ _ 

123  DO  1 24  J  =  1 *  NP 

124  CALL  CHAOS(T*X*Y*VK*SW»J*BBMIX*ZC*AAOB*A*TC*PC*W* 

1DFLTA*V1 »DDELMX )  _ _ _ _ 

GO  TO  78 

77  DO  7  J=1 *NP 

DO  7  1  =  1  *NCP _ 

CALL  KVAL(T,VK* I *J*SW) 

7  CONTINUE 

78  DO  8  J=1 >NP 
G ( J ) =0  •  0 

8  U(J)=0.0 

C  
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U  REPRESENTS  LIQUID  SIDE  STREAMS.  U(l)  IS  THE  DISTILLATE 
C  RATE  IN  MOLES  PER  UNIT  TIME. 

C 

G{ 16  )=?5 .0 
U  (  4  )  =  1  5 . 0 

_ DO  17  1  =  1, NCP _ _ _ 

DO  10  J  = 1 » N  P 

10  R ( J ) = — (  ( V ( J )  -  i  ( J ) ) *VK ( I  * J )  +  * PROFL ( J )+U ( J ) ) ) 

F ( 1 )=0.00 
DO  11  J  =  2  *NP 
JJ=J-1 

_ II _ E  (  J  J  =PROFL  (  Js. _ 

DO  12  J=1.*NPP 
12  r( J)=V( J+l }*VK( I»J+1) 

C ( NP ) =0.00 

c 

C  THE  ABOVE  MATRIX  IS  SOLVED  FOR  THE  COMPOSITION  OFILE  BY 

C  A  GAUSS  ELIMINATION  METHOD._ _ _ 

C 

P( 1 )=C(1  )/B(l  ) 

Q(  1  )=D(  1  *1  )  / R  C  1  )  . . . 

DO  14  J  =  2  * N P P 

14  P(J)=c(J)/(B(J)-(E(J)*P(J-l)>) 

_ P  (  NP_1=_Q . 00 _ _ _ 

DO  15  J  =  2  » N  P 

15  Q( J) = ( D ( J » I )-( E ( J) *Q( J-l ) ) ) / ( B ( J )-( E( J >*F  ;-l ) ) ) 

X (  I  ,NP ) =0 ( NP ) 

DO  16  J= 1 *NPP 
L  =NP  — J 

_ 16 _ X  (  I  ,  L  )  =0  ( L ) - ( P ( L ) *X (  I  »  L  +  l  )  ) _ _ 

17  CONTINUE 
C 

C  COMPOSITION  PROFILE  IS  BEING  NORMALIZED 
C 


DO  19  J=1»NP 

SUMX { J) =0.00 

DO  19  I  =  1 » NCP 

SUMX (J)=SUMX(J)+X( I ,J) 

19 

CONT  I NUF 

DO  97  J  = 1  , N P 

DO  97  1=1 ,NCP 

X (  I  , J ) =X ( I , J) /SUMX ( J) 

97 

CONTINUE 

DO  75  J  = 1 » N P 

DO  75  1=1 *NCP 

# 

75 

Y  (  I , J)=VK( I ,J)*X ( I  ,  J ) 

IF ( SW.EQ.O. 0 )  GO  TO  125 

GO  TO  127 

125 

DO  126  J=1,NP 

126 

CALL  CHA0S(T»X,Y,VK,SW,J,8BMIx,ZC,AA0B,A,TC,PC,W, 
1DELTA , VI , DDELMX ) 

SW=SW+1 .0 

JJJ=SW 

GO  TO  78 

127 

SW=SW+1 .0 

JJJ=SW 

CALL  TNEW ( X  » J  » AT  *  SW  *  Y  *  T  >  A »  TC  »PC  » W  *  DELT  A  * VI  * AAOB  »  BBM I X  » 
1 ZC  » DDELMX ) 

DO  131  J  =  1  ,NP 
DO  131  1=1, NCP 


. 

or  oo  (o.n.oa.w?.)ii 
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131  Y(  I  ,  JJ=X ( I  *  J ) *VK ( I ,  J) 

DO  150  J=1»NP 

SI )MY  (  J  )  =  0  o  00 

DO  150  1=1 *NCP 

SUMY ( J ) =  SUMY ( J )+Y (  I  * J ) 

150  CONTINUE  _ 

DO  151  J  =  1  *  NP 

DO  151  I  =  1  *  N  C  P 

Y(  I , J)=Y ( I  *  J ) /SUMY ( J ) 

151  CONTINUE 
SUMDT  =  0 . 0 

_ DO  .47  J=1  » NP _ 

47  SUMDT=SUMDT+( AT ( J)-T( J)  )**2 
IF(SUMDT.LE.l.O)  GO  TO  37 
GO  TO  43 

48  DO  49  J=1»NP 
T ( J ) = AT ( J ) 

49  CALL  CHAOS  f  T  »  X  »  Y»VK»SW»J*BRMI  X  * Z  C  *  A  A  0  R  ,  A  » T  C  , P  C  *  W.., _ 

1DELTA*V1 *DDELMX ) 

CALL  PROFYL(V*PROFL, J  * X  * T ♦ SW * Y *  A , TC *PC * W ,DELTA * VI *AAOB* 
1BRMIX  *ZC*DDELMX »NP ) 

CALL  FPTSW 
GO  TO  92 

37 _ L'R  IT  E.{  6  ♦  106) _ : _ 

WR ITE(6*101  ) 

WRITE! 6 ♦ 102 ) 

DO  904  J=1 *NP 

904  WRITE (6»103)  J ♦ T ( J ) , V ( J ) > PROFL ( J ) 

WR I TE ( 6 ♦ 1 04 )  JJJ, SUMDT 

_ W  R_IJ.EJL6_i.lJJJ, _ 

WRITE (6 *112) 

WRITE (6 *115) 

DO  903  1=1*11  _ 

905  WRITE <6 *113)  X(I,1),X(I*21)»X(I,4)*Y(I*16) 

STOP 

END 
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SIBFTC  RAY  NODECK *NOL I  ST 

SUBROUTINE  KVAL ( T * VK , I , J ,SW ) 

DIMENSION  T ( 25 ) >VK ( IS »?5 ) »AA ( 15 >  *AR ( IS ) *  AC ( 15 ) 

DIMENSION  ALPHA (  1 5  >  2  5  ) 

C 

C  THIS  SUBROUTINE  CALCULATES  IDEAL  COMPOSITION-INDEPENDENT 

C  VAPORIZATION  CONSTANTS. 

C 

101  FORMAT (1X*F8.5*F8.5»F7.5) 

IF(SW.EQ.O.O)  GO  TO  70 
GO  TO  11 

70  I  F_(  I  •  EQ  •  1  .AND •  J  •  EQ  •  1_) _ .GO_LO_JLO _ 

GO  TO  11 

10  DO  12  11=1*11 

12  READ ( 5  *  1 01  I  AA( II )  >AR( II > *AC< II> 

A A ( 9  )  =  — A A ( 9 ) 

AA ( 1 0 ) =-AA ( 1 0 ) 

_ AA(.11)=-AA(  11.1 _ 

AB ( 1 ) =-AR ( 1 ) 

AB ( 2 ) =-AB ( 2 ) 

ABO)  =— AB  (  3  ) 

A  R  (  4  )  =  —  A  B  (  4  ) 

AB ( 5 ) =-AB ( 5 ) 

_ AB  ( 1 1  )=-AB  ( 11  ) _ _ _ _ _ _ 

AC ( 7 ) =— AC ( 7 ) 

AC ( 8 ) =— AC ( 8 ) 

11  ALPHA ( I  9  J ) = A A (  I ) +AB ( I  > * ( T ( J  > / 1 00 . 0 ) +AC ( I ) * (  ( T ( J ) / 1 00 • 0  > **2  > 
VK (6  » J ) =“0. 0419  5+0. 34  08  7* ( T (J)/100.)-0.01815*((T(J)/100.) 

1**2 )+0. 02872* ( (T(J)/100.)**3) 

_ VKJIvJ)=  VK(  6»  J)*ALPHA  (  I » J) _ _ _ 

RETURN 

END 
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SIBFTC  HELEN  NODECK *NOL I  ST 

SUBROUTINE  CHAOS ( T  >X  > Y  ,  VK , SW , J , BBM I  X , 2C  ,  AAOB ♦ A > TC , PC * 

1W, DELTA, VI  * DDELMX ) 

C 

C  THIS  SUBROUTINE  CALCULATES  COMPOSITION  DEPENDENT  K-VALUES 

c  by  ihe  chao-seader  method. _ 

c 

DIMENSION  T(25)>X(15>25)*Y(15>25)*VK(15*25)>A(17*5) 

DIMENSION  TC( 15 >  *PC( 15) >W( 15 )  >DELTA( 15 ) *V1 (15 )  >FNU(15>25) 
DIMENSION  ARK (15) *  RRK (  1 5  )  *  PH  I  (  1 5  *  2  5  ) *GAM( 15  *25 )  >VAPL (15*25) 
DIMENSION  AAQR ( 25 ) *BRMIX ( 25 ) *2C ( 25 ) » DDELMX (25 ) *XXF (15*25) 
DIMENSION  YYF( 15 *Z5 ) *TE(25 )  .  _ 

200  FORMAT ( IX »40HMCRE  THAN  THIRTY  TRIALS  ARE  REQUIRED  FOR* 

11  AH  2  CONVERGENCE) 

201  FORMAT ( 1 X *F7 . 2  * F6 . 1  * F6 .4  * F5 . 3 ,F5 . 1 L 
19R  FORMAT ( 1 X  * 2E 1 5 . 8  ) 

NCP=11 

I F ( SW . E Q . 0 .0. AND. J.EQ. 1 )  GO  TO  299 _ _ 

GO  TO  300 

299  DO  75  1=1*16 

75  RFAD(5»199)  (A( I»IJ) *IJ=1»2)  \ 

READ (5  >2  01 )  ( TC( I )  »PC ( I  ) * W (  I ) , DELT A ( I ) , VI ( I ) , I  =  1 *NCP) 

300  P=250.0 

_ T  ( J_)_=_Ti.J  )  +459 .6  _ : _ : _ 

DO  101  I = 1 >  NCP 
IF (  I  .EQ. 1  )  GO  TO  11 

IF(I.GT.l)  GO  .T0..1J2 _ _ _ _ _ _ 

11  IJ=1 
GO  TO  14 

12  _ I  J  =  2 _ _ _ 

GO  TO  14 

14  TR=T( J)/TC( I ) 

PR=P/PC (  I  ) 

FPO  =  A ( 1 > I J)+A( 2  * IJ ) /TR+A (3  » I J)*TR+A (4* I J)*TR**2  +  A( 5 , I J)*TR**3 
l  +  ( A(6*IJ )+A(7* IJ)*TR+A(8>IJ)*TR**2 )*PR+< A ( 9 *  I J)+A( 10* IJJ*TR) 

_ 1  *P  R**2  -  A  LOG  (PR)  /  2_.  302.5.8 _ 

FPI=(A< 11*IJ)+A( 12* I J ) *TR+A (13*1 J) /TR+A ( 14 ♦ I J )*TR**3+A ( 15  *  I J > 
1* ( PR-0 . 6  )  ) 

101  FNU ( I  *  J )  =  ( FPO+W ( I ) *F P I )*2. 30258 _ 

AM  IX  =  0. 0 
BM  I X  =  0 . 0 

_ DO  102  1= 1 >NCP _ 

TR=T( J) /TC( I ) 

ARK ( I )= (0.4278/ (PC ( I )*TR**2.5) )**0.5 

BRK  (  I  )  =  0 .0867/  {  PC  (  I  )*TR  ) _ _ 

RM I X  =  BM I X+BRK ( I )*Y( I  * J) 

102  AM  IX=A.M  I  X  +  ARK  (  I  )  *Y  (  I  »  J  ) 

_ A  0  B  =  AMJLX  *  AMIX/BMIX _ 

CON  1=2.0 

C0N2=BMIX*P +0.001 
DO  106  K K  =  1  ,30 
2 AS= ( CON 1+CON 2  )  /2.0 
H=BM I X*P /2 AS 

_  Z  =  l  ,0/(  1  .0-H) -AOB*  (  H/  (  1.0+HM _ 

DEL=2-2 AS 

IF ( ABS(DEL) .LE.1.0F-04)  GO  TO  107 
IF(DEL.LT.O.O)  GO  TO  104 
IF (DEL. GT. 0.0)  GO  TO  105  ' 

104  CON  1 =2 AS 
GO  TO  106 


-  - 
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105  CON?=ZAS 

106  CONTINUE 

107  AAOB ( J ) =  AOB 
BRM  I X ( J ) =BM I X 
ZC ( J)=Z 

DO  108  1=1,  NCP 

108  PHM  I,J)  =  (Z-]  .0)*(BRK( I  )/BMlX)-ALOG(Z-RMlX*P)-AOR*(  (2.0* 
1 ARK ( I) ) / AM  I X- ( BRK ( I )/BMIX)  ) *ALOG ( 1 • 0  +  H ) 

C1=0.0 
C2=0 • 0 

DO  10Q  1  =  1  » NCP 

D  E  L  M  X  =  X  (  I  ,  J)*V1  (  LI _ 

C1=C1+DELMX 

109  C2=C2+DFLMX*DFLTA (I ) 

DFLMX=C2/C1 

DDFLMX ( J ) =DFLMX 
C1=1.987*T( J)/1.8 

_ DO  110  1=1, NCP  _ •_ _ _ 

GAM (  I  ,  J ) =Vl (  I ) *( DELTA ( I ) -DPLMX  >**?/Cl 
VAPL( I ,J )  =  ( F N U ( I , J ) +G AM ( I ,J)-PHI ( I  ,J)  ) 

110  VK( I *J)=EXP(VAPL(  I  ,J)  ) 

T ( J ) =T ( J ) -459 .6 
RETURN 

_ END _ 


$  I BFTC  FDWARD  N0DFCK  *NOL I  ST 

subroutine  TNEW(X,J,AT*SW,Y,T,A*TC*PC*W*DELTA,V1 ♦AAOr, 
1RRMIX.ZC.DDELMX ) 

DIMFNSION  TT ( 2  5  )  *T(25)*AT(25)*X('15*25>  *VK (15*25)  *D( 25 ) 
DIMENSION  S ( 20  *  20 ) *  SUMY (25)  *Y(15»25).XY(25),A(17»5) 


r 

Clio  1  WIN  \  ID  >  »  ^  V-  V  I  T  '  >  WJ  ID  >  »  UlL  1  rt  1  1  “  1  >  V  ±  l 

DIMFNSION  BRM I X ( 25 ) »ZC(?5) *DDELMX ( 25 ) 

ID  >  *  MMV  D  \  <LD  •  __  _  . 

c 

c 

r~ 

THIS  SUBROUTINE  USES  MULLER* S  QUADRATIC  ROOTS 
CALCULATE  THE  NEW  TEMPERATURE  PROFILE. 

method  to 

.  N  C  P  =  1 1 

57 

J=1 

JJ  =  1 

YX  X  =  T  ( J  )  .  ... 

14 

YX  =YXX 

F J J= JJ— 1 

TT(JJ)=(T(J)-10.0)+FJJ*10.0 

X Y ( J ) =TT ( J J ) 

CALL  CHAOS(XY,X,Y,VK,SW*J»BBMlX*ZC*AAOB*A*TC*PC*W* 

1DELTA.V1 .DDELMX )  \ 

12. 

SUMY ( J ) =0.00 

DO  12  1=1 *NCP 

SUMY( J)=SUMY( J)+VK(  I  * J)*X( I *J) 

4 

S( J*JJ)=SUMY( J)-1.00 

JJ=JJ+1 

IF (JJ.EQ.4)  GO  TO  19 

19 

GO  TO  14 

K  =  4 

K  I  K  =  1  _  _  __  __  .....  .  ..  _______ 

99 

D(K-1 )=(TT(K-l)-TT(K-2) ) /(TT(K-2)-TT(K-2) ) 

R  =  S( J>K-3)*(D( K-l ) **2 )-S(J*K-?)*( 1.0+D(K-1 )  ) **2  +  S ( J  ♦ 

1  K—  T  )*M  •  0+  (  2  .  *D  (  K- 1  )  )  ) 

C=S( J»K-3)*D( K-l )-S( J,K-2)*( 1. 0+D ( K-l ) )+S( J 
E=-2.0*S( J.K-1 ) * ( 1 • 0+D ( K-l ) ) 

AA=B**2-4.0*S(  J>K-1  )  *D  (  K-l )  * (  1 . 0  +  D  (  K-l  )  )  *C 

*  K-l ) 

ABP=B+ ( ARS ( AA ) **0.  5 ) 

ABM=B-< ABS(AA)**0.5) 

RC1=E/ARP 

RC2=E/ABM 

TA  =  TT(K-1)  +  (TT(K-1 ) — T  T ( K— 2 )  )  *BC  1 

TR  =  TT ( K-l )+(TT(K-l ) -TT ( K  — 2 ) )*BC2 

IF(TA.GT.1000.0. AND. TB.GT. 1000.0)  GO  TO  120 

IF (TA. LT. 0.00. AND. TB.GT. 1000.0)  GO  TO  120 

IF ( TA.GT. 1000.0. AND. TB.LT. 0.00)  GO  TO  120 

IFITA.LT. 0.00. AND.TR.LT. 0.00)  GO  TO  120 

IF (TA.LT. 0.0. OR. TA.GT. 1000.0 )  GO  TO  70 

GO  TO  87 

120 

87 

TA  =  YX 

GO  TO  98 

T  (  J  )  =T A 

CALL  CHAOS (T  *X*Y*VK,SW* J *BBMIX  * ZC * AAOR » A *TC 
1DELTA»V1 *DDELMX ) 

SUMY { J ) =0.00 

*PC.W* 

81 

DO  81  1=1 *NCP 

SUMY( J)=SUMY( J)+VK( I  *  J ) *X ( I  * J ) 

SA=SUMY( JJ-1.00 

70 

IF (TB.LT .O.O.OR. TB.GT  .1000.0)  GO  TO  92 

GO  TO  70 

T ( J  )  =TB 

' 
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CALL  CHAOS ( T  *X  » Y  * VK  *SW  * J ,RRM I X  > ZC ♦ A AOR , A ♦ TC ♦ PC *W ♦ 
1  DELTA  * VI *DDELMX ) 

SUMY( JJ =0.0 
DO  83  1=1 *NCP 

S3  SUMY ( J ) =  SUMY ( J ) +VK ( I  * J)*X( I *J) 

_ S  B  =  S  U  M  Y  (  J  )  - 1 . 0  0 _ 

IF(ABS(SR) •LT.ABS(SA) )  GO  TO  91 
GO  TO  92 

91  TF( ARS(SB) .LF. 0.001 )  GO  TO  97 
K  =  K+1 

TT (K-l ) =TB 

_ S( J*Krl ) =SB _ 

KIK=KIK+1 

IF (KIK.EQ.20)  GO  TO  120 
GO  TO  9° 

92  IF( ABS(SA) .LE. 0.001 )  GO  TO  98 
K  =  K  +  1 

_  TTJ.K-1  )=TA _ . _ ; _ _ 

S ( J  »  K-l ) =S A 
K I K  =  K I K+ 1 

IF (KIK.EQ.20)  GO  TO  120 
GO  TO  99 

97  AT ( J ) =TB 

_ T ( J ) =YXX _ 

GO  TO  15 

98  AT  (  J  )  =  T A 

T  (  J  )  =YXX  . .  . 

GO  TO  15 

15  IF (J.EQ.21)  GO  TO  56 

_ G 0  TO  58 _ 

58  J=J+1 

GO  TO  57 
56  RETURN 
END 
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SIBFTC  COt  ITU  N0DECK*N0LIST 

SUBROUT  I NE  FN  THAL (Y»X»T *BBM I  X ♦ ZC  *  A AOB ♦  A  *  TC  * PC ♦ EV ♦ EL * SW ♦ 

1V1 *DDFLMX»DELTA*J*W) 

C 

C  THIS  SUBROUTINE  CALCULATES  THE  LIQUID  AND  VAPOR  STREAM 
C  ENTHALPIES  BY  THE  CHAO-SEADER  METHOD. 

C 

DIMENSION  TC ( 15  > »PC ( 15  J  * A( 17  *5 ) *HFC ( 15  >  *HAC ( 1 5 )  *Vl ( 15  > 
DIMENSION  DELTA (IS) *DDFLMX(25) *HL(15)*T(25)*W(15>*EV(25) 
DIMENSION  EL ( 2  5 ) *CH( 1 5  >  *XXF ( IS ,25 )  *  YYF (15*25)* TF(25)*TT (25) 
DIMENSION  HIDL(15)»AH(15)*BH(15)*DH(15)»Y(15*25) 

DIMENSION  B  (  1 5  »  2  5  L»  BBM I  X  L25±_*  ZC  (  25  )  *  A  AOB  (  25l*  X  (  15  *  2  5  J _ 

TT ( J  )  =T ( J5+459.6 
1C]  FORMAT ( 1X*4E14.7) 

ICS  FORMAT ( IX* 12*2 (2X»F10.3 ) ) 

P  =  2  5  0  •  0 
NCP= 1 1 

_ I.FJ  S.W. .  EQ .  2  .O.AND.  J.EQ  ..l.L_.GO . TO— SLS _ 

GO  TO  100 

9Q  READ ( 5  *  1 01 )  ( AH ( I ) * BH ( I > *CH ( I ) , DH ( I > ♦ I =1 *NCP ) 

100  DO  1  1=1 *NCP 

1  H I DL ( I ) = AH ( T )  +BH  ( I )*T( J)+CH( I ) *T ( J ) **2+DH ( I )*T( J)**3 

SUMHV  =  0  o  0 

_ _ _ DO _ 3 _ 1  =  1  »NCP_ _ . _ _ _ _ _ ___________ 

2  SUMHV=SUMHV+Y( I *J)*HIDL( I ) 

H=BBMIX(J)*P/ZC(J) 

HVC=(1.5*AAOR( J)*ALOG( 1 . 0+H ) +1 .0-ZC ( J ) ) *1 . 987*TT ( J ) 

EV( J)=SUMHV-HVC 
SUMHL  =  0  o 0 

_ D0_4__  1  =  1  ♦  N  CP _ 

I F  (  I  .EO.-l  )  GO  TO  11 
I F ( I «GT  « 1 )  GO  TO  12 

11  IJ  =  1  _ _ _ _ _ _ _ 

GO  TO  16 
12  I J=2 

_ GO _ TQ  .16 _ 

16  TR=TT(J)/TC(I) 

PR  =P /PC (  I  ) 

CON  1 =  — A ( 2  *  I J ) /TR**2  +  A (2* I J)+2.0*A(4  *  I J )*TR+3.0*A ( 5  *  I  J)  _ 

l*TR**2+( A(7* I J)+2.0*A(R*I J)*TR)*PR+A( 1 0 ♦ I J ) *PR** 2+ ( A ( 1 2 
2*IJ)-A(13*IJ) /TR**2+2.0*A( 14*1 J)*TR**2 >*W( I ) 

_ CON  1=  CON  1*2.302 58/TC(  I  ) _ _ 

HFC ( I ) =CON 1 

HAC ( I )=-( (VI ( I )*( DELTA ( I ) “DDELMX ( J ) ) **2 ) / ( 1 .987* ( TT ( J > 

i**2.n  )*i  .8  ..  . . . . . . 

HL ( I ) =— ( HFC ( I )+HAC(  I )  )  *1  .98  7*TT ( J ) **2  +  H I DL ( I ) 

4  SUMHL  =  SUMHL+X ( I  * J)*HL ( I) 

_ EL ( J ) =SUMHL _ 

T( J)=TT( J)-459.6 
RETURN 
END 
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SIBFTC  KATHY  NODECK ♦ NOL I  ST 

SUBROUTINE  PROFYL (V,PROFL,J,X,T,SW,Y,A  *tc*pc,w,delta,vi * 

1 AAOB  *BBMIX*ZC*DDFLMX»NP ) 

C 

C  THIS  SUBROUTINE  CALCULATES  NEW  VAPOR  AND  LIQUID  PROFILES 
C  BY  MAK  ING  CONVENT  I  ONA L  ENTHALP Y  BALANC E3_._ 

C 

DIMENSION  FFV( 25  )  ,  FEL (25)*Y(15*25)  , SUMY ( 25 )  ,T (25 ) 
DIMENSION  SH( 15*25 ) *CH( 15*25 ) *EV(25  >  *EL (25  ) *X (  15  *25  ) 
DIMENSION  V ( 25 ) *PROFL(25) *F ( 25 ) *0 ( 25 ) *P ( 25 ) *U ( 25 ) *HF ( 25 ) 
DIMENSION  EF1(15)*EF2(15)*VK(15*25)*A(17*5)*TC(15)*PC(15) 
DIMFNS ION  W ( 1 5 ) »DFLTA( 15) »V1 (15) ,  APR ( 2 5 ) * BBM I X ( 25 ) 
DIMENSION  ZC( 25  )  *DDELMX( 25)  *TF(25)  ,XXF(  15*25) *YYF( 15,25) 

105  FORMAT(1X*I2*2(2X*F10,3) ) 

106  FORMAT ( 1 X  »F7. 5  ) 

108  FORMAT ( 1 X , 1 HJ ♦ 5X , 1 HT * 7X * 1HV * 7X * 1HL ) 

109  FORMAT ( IX* 1 2 ♦ 1 X , F7 . 2  *  1 X  * F7 . 2  *1X,F7.2 ) 

_ NLP  *2 1 _ 

NCP=1 1 
NPP=NP— 1 

IF (SW.FQ.2.D)  GO  TO  71 
GO  TO  70 

71  DO  85  I = 1 *NCP 
_ DO  85  J= 1 *NP 

XX F (  I  ,J)=0#0 
85  YY  F ( I  *  J ) =0 • 0 
DO  69  1  =  1, NCP 
J=  1 3 

69  RE AD ( 5 , 1 06 )  XXF(I,J) 

75  DO  76  1=1, NCP 

J  =  7 

76  READ ( 5 , 1 06 )  YYF(I,J) 

70  DO  19  J=1,NP 
SUMY( J) =0.00 
DO  19  1=1 *NCP 

19  SUMY( J)=SUMY( J)+Y( I ,J) 

DO  67  J=1 *NP 
DO  67  1=1, NCP 
67  Y (  I,J)=Y(I,J) /SUMY ( J ) 

DO  74  J= 1  ,NP 

74  CALL  ENTHAL( Y ,X , T ,BBM I  X ,ZC , A AOB , A  * TC ,PC , EV , EL ,SW ♦ V 1 , 
1DDELMX,DELTA*J,W) 

V(  1  )  =23.0 
V( 2 )=94.806 
PROFL ( 1 ) =71 *806 
PROFL (21 ) =37 • 0 
DO  51  J= 1 *NP 
F ( J ) =0 • 00 

TF ( J ) =0 • 0 
Q( J) =0.00 
P( J)=0.0 
U( J)=0.00 
51  HF ( J ) =0 • 0 

_  F( 7) =34.50 _ 

F( 13)=65.50 
U(4)=15.0 
P( 1 6 ) =25  .0 
J  =  7 

TF(7)=117.7 

_ CALL  CHAOS ( TF , XXF , YYF , VK , SW , J, BBM IX ,ZC * AAOB , A *TC ,PC*W, 
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(<;.v^*xr*s#v^,xi#^.v^.x[fsifxr)  TAvqoq  p°  r 

K  =  qn 
rr=q:>n 
f-qM=qqn 
rv  ot  no  (o.s.oq.w? t qi 

oy  ot  oo 

qoHff=i  ?r  oa  rv 
qw* r=L  p«  no 
o.o* ( t. i ) oxx 

0 .Os ( L* I ) 1 YY  PP 

qon. r  =  i  Pd  oa 
p  r=u 

(L.nixx  (<^or*?)oAqq  pd 
qoHff*i  dY  oo  ?y 


Y  =  L 

( l *  t  ) qy y  (^or^)OAqq  dT 

qnt  r=L  p  r  oo  oy 

on. os  (L)YMi ip 

qo!'*.r  =  i  Pr  oo 

(L»  I  n  +{UYv:J?=(UY'''p  of 

qn«  r*L  vs  oa 
qou*  r  =  i  Yd  oa 

(U  YMUPN  (Lf  T  )  Y=(L»  I  )  Y  Yd 

q^t  r =l  ay  oo 
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0«FS«(DV 
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pia.  r  =l  r 3  oo 


o  o .  - .  =  {  l  i  q 
0.0*(L)qT 
oo.o=(L)0 
o.o=(i>q 

00.0= (L )U 
0 . 0* ( L ) 1H  [? 

o  p  .  A  P  =  (  Y  )  ^ 

od. Pd*( p  r  )q 

0.  r>  [  =  (  A  H) 

o. py=  ( d  r ) q 

Y  =  L 

y.  y  r  r  =  ( y ) it 

.  Oq.DT,  Af.«-  A4,r  .  v  I  V  jf- ,  L  •  •'P  •  XV  »  ’YY.qXY.qj  )  P.OAHO  J  J  A  ) 


1DFLT  A* VI ♦DDFLMX ) 

CALL  ENTHAL( YYF  ,  XXF  ,  T F  , BBM I  X , ZC , AAOB ♦ A ♦ TC ,PC , FEV *FEL * 
1SW*V1 *DDELMX » DELTA » J * W) 

J  =  1 3 

TF ( 13 ) =268.5 

CALL  CHAOS (TF*XXF*YYF*VK>SW* J * BBM I  X *2C » AAOB * A  * TC *  PC  * W * 
1DELT A  » VI *DDELMX ) 

CALL  ENTHAL ( YYF * XX F * TF * BBM I  X , ZC * AAOB *  A  * TC *PC * FEV » PEL , 
1SW*V1 *DDFLMX» DELTA *J*W) 

HF  ( 7 ) =FFV ( 7 ) 

HF ( 1 3 )  =  FEL ( 1 3 ) 

.  DO  115_J=2 »NPP  _  _ 

XBX= ( EL ( J-l )-EL ( J) )*PROFL( J-l ) 

RXX= ( EL ( J ) -EV ( J ) )* ( V( J)+P ( J ) ) 

BBX= ( HF ( J ) — EL ( J ) ) *  F ( J ) 

BBB*<EL< J)-EV( J+i) ) 

5BV=(X8X+BXX+BBX-Q(J) ) /BBB 

_ BA R  =  A R S. (  (  BRV— V  (  J-f  1  )  )  /V(  J  + 1  )  ) _ _ 

IF  (BBV.LF.0.0.OR.BAB.GF.0.40)  GO  TO  HO 
GO  TO  111 

110  V( J+l )=V( J+l ) 

GO  TO  112 

111  V( J  +  l )  =  BRV 

_ GO _ 1 0_1 1  2 _ 

11?  RBL=F( J)+PROFL ( J-l ) +V ( J+l ) -V ( J ) -P ( J ) -U ( J ) 

ABA= ARS ( (BBL-PROFL ( J) )/PROFL(J) ) 
IF(RBL.LE.0.0.OR.ABA.GE.0.40)  GO  TO  113 
GO  TO  114 

113  PROFL { J ) =PROFL ( J ) 

_ G  0  __LQ__1  JL5_ _ 

114  PROFL ( J ) =BBL 
GO  TO  115 

115  CONTINUE  .  _ _ _ _ _ _ _ 

RETURN 

END 


, 
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Cards  to  be  changed  when  problem  III  is  used  with  a  different 


system : 

Mainline  Program  (Thesis) 

1.  The  number  of  plates  NP  =  21 

2.  The  number  of  components  NCP  =  11 

3 .  The  locations  of  the  feeds  JI  =  7 

NJ  =  13 

4.  The  locations  and  flow  rates  of  the  vapor  and  liquid 
side-streams  respectively  6(16)  =  25.0 

U  ( 4 )  =  15.0 

$  IBFTC  Helen 
Chao-Seader  Subroutine 


1. 

The 

number 

of  components 

NCP  =  11 

2. 

The 

column 

pressure,  atm. 

P  =  250.0 

Temperatures  coming  into  this  subroutine 
have  to  be  in  °F. 

$  IBFTC  Edward 

Subroutine  for  new  temperatures 

1.  Number  of  plates  IF(J.EQ.  21)  GO  to  56. 

$  IBFTC  Coutu 

Subroutine  for  enthalpy  calculations 

1.  Column  pressure,  atms .  P  =  250.0 

2.  Number  of  components  NCP  =  11 
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$  IBFTC  Kathy 
Subroutine  for  new  rates 


1. 

Number  of  plates 

NP 

= 

21 

2. 

Number  of  components 

NCP 

— 

11 

3. 

Liquid  feed  introduction,  moles/U.T. 

J 

= 

13 

4. 

Vapor  feed  introduction,  moles/U.T. 

J 

= 

7 

5. 

Top  vapor  product,  moles/U.T. 

V(l) 

= 

23. 

0 

6. 

Condenser  input,  moles/U.T. 

V  ( 2 ) 

— 

94. 

806 

7. 

Condenser  recycle  stream, 

PROFL(l) 

= 

71. 

806 

8  . 

Bottom  product 

PROFL (21) 

= 

37. 

0 

9. 

Feed  locations  and  rates 

F  ( 7 ) 

= 

34. 

50 

10. 

Sidestream  locations  and  rates 

U  ( 4 ) 

15. 

0 

P(16) 

rr 

25. 

0 

11. 

Feed  temperatures  and  locations 

J 

= 

7 

TF  (7) 

= 

117 

.7 

J 

= 

13 

TF  (13) 

— 

268 

.  5 

12.  Setting  of  feed  enthalpies  to  the  calculated  ones 

HF ( 7 )  =  FEV ( 7 ) 

HF ( 13 )  =  FEL ( 13 ) 


Data  Cards  to  be  Changed 

1 .  Initial  Profile 

X01080. 023. 000071. 806 

Format  (IX,  12,  F5.1,  F6.3,  F7.3) 
Stage  number  =  12  =  01 

Temperature  =  F5.1  =  080. 0°F 
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(f)V3q 
(£X) J33 


(f  )3H 


(£I)3H 


d08.Ifp000.£SO.080X0X 

(£.fq  ,£.d3  yL.tt  %SI  ,XI)  d*flno3  , 

^scbnun  sped 2 


3°0.080  *  1.23 


9‘iud639qm9T 


Vapor  rate  =  F6.3 

Liquid  rate  =  F7.3 


23.000  moles/U.T. 
071.806  moles/U.T. 


2 .  Distribution  of  components  in  the  two  feeds 

X-02. 0-00. 0 
Format ( IX , 2F  5.1) 

The  rate  of  introduction  of  first  component  is: 

2.0  moles/U.T.  in  first  feed 
0.0  moles/U.T.  in  second  feed 
The  next  card  represents  the  rate  of  introduction  of 
second  component. 

Negative  signs  arise  from  the  setting  of  the  matrix 

solutions  (-F.Z. .) 

3  13 

3 .  Equilibrium  data  for  first  iteration 

In  this  case,  relative  volalities  are  calculated  and 
component  six  is  the  reference  component. 

4 .  Chao-Seader  data  for  the  components  in  the  system 

X034 3. 900673. 10. 00005. 680052.0 

Format  (IX,  F7.2,  F6.1,  F6.4,  F5.3,  F5.1) 

Critical  temperature  °R  =  F7.2  =  0343.90 

Critical  pressure,  psi  =  F6.1  =  0673.1 

=  F6.4  =  0.0000 


Acentric  factor 
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08.C&E0  *  S.V1  ■  fl°  9iLr^Bi9qm9^  Xbox^X^D 
X.CVdO  *  X.d'i  *  xaq  ,9^uea9xq  XboX^x^O 


0000.0 
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6  . 
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F5.3 


5.680 


volume  parameter  =  F5.1  =  052.0 

This  information  is  for  component  1  only. 


5 .  Composition  of  the  feeds 


X0. 00000 

Format  (IX,  F7.5) 

This  means  that  in  the  first  feed,  there  is  none  of 
component  1 . 

The  composition  of  the  first  feed  is  read  in  for 
every  component  then  that  of  the  second. 

6 .  The  parameters  for  the  Chao-Seader  component  enthalpies 
Format  (IX,  4E  14.7) 

AH,  BH,  CH,  DH  are  listed  by  component 
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00000. ox 
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TTRFTC  THFSTF  NOnFCK 

COMMON  J,  IW,NP,NCP,KK,  IK  ,FLOD ,  I 

COMMON  V  (P5  )  , PROFL  (?*)*VK(1.*>,?5),X(1S,;>«5) 

COMMON  FtJMX  (  7  6  )  ,RRM  I X ( ?6 )  ,  ZC( 7F )  ,  A  AO^ ( os ) 

COMMON  A(17»S)  , DFL  T  A  (  1  F  )  ,V1  (  1  F  )  ,DDFLMX  (  ?F  )  . F  f  ?F  )  .1  If  OF  ) 

98 

0 1 MFNF I ON  G(?f)  ,R ( ?F  )  , F ( 7P  >  ,C  ( ?F  )  ,  P  (  ?f )  ,Q( ?f  )  ,AT(?f)  , 
1D(?5),T(?5) 

FORMAT ( 1H1 »44X , 18HPR09LFM  NO.  6) 

09 

100 

101 

FORMAT ( 1X,I?,F5.1 , F8  .  ] *  F5  •  1 ) 

FORMAT ( THJ,48X,1fhTNTTIAL  PROFILF) 

FORMAT ( 1 HK,?OX ,FHSTAGE»6X ,4HTFMP , 1  ox , 1 1 H VAPOUR  RATF  .  ov 

107 

100 

1*11 HL IQUIO  RATF ) 

FORMAT ( IX  *o?X» 1 HF, 1 7X  *  1 OHMOLEF/U.T. *  1 ox ,  lOHMOLFF/l  i  ,T . / ) 

FORMAT (1X,70X,I7«FX,fq.o,7(6X,F14.7) ) 

104 

10$ 

FORMAT ( 1 HK ,71 X *?ohFUM  (DELTA  T**2)  AT  I T FR A T I  ON , 1 X , I 0 , 

1 IX , 1H=,F6.3 ) 

FORMAT ( 1 HK ,44X , 1 RHFTNAL  PROF  ILF ) 

1 11 
117 

FORMAT ( 1  Hi  ,01 X , 77HPRODUCT  STRFAM  COMPOSITIONS) 

FORMAT ( 1 HK,77X,9HC0MP0NFNT , OX, 1 6HCQNDFNSFR  STREAM, 6 X ,1 

7  "3  upo  T  TQM  FTRFAM) 

118 

114 

JL1A 

FORMAT ( 1HJ,04X,F14.7,7X,F14.7) 

FORMAT  (  1  HK  ,00X  ,  0  i  hRECOVFRY  OF  METHANOL  ACHIEVED  =  ,  F  8  -  ■?  ) 

FORMAT ( 1  HI  ) 

NCP  =  0 

N  P  =  7  0 

NPP=NP-1 

DO  8  J=1  ,NP 

G  (  J  )  =  0 . 0 

r( J)=G#Q 

C 

C 

8 

I !(  J  )  =0.0 

INITIAL  PROFILF  IF  R  F  A  D  IN 

c 

Rfao(F,oo)  ( J,T ( J)  ,PROFL ( J )  , V ( J ) , J  =  1  , N P ) 

DO  7 00  J  =  1  ,NP 

c 

c 

pnn 

T ( J ) = ( T ( J )*1 • 8 )+o? .O 

KK 

IF  SOLVENT  FEED  LOCATION.  IK  IF  FRESH  FFFD  LOCATION. 

c 

c 

c 

F  Ic  FFED  RATE  IN  MOLES/U.T..  FLOD  IF  L/D  RATIO.  U(l)  IF 
DISTILLATE  rate. 

HR  I T E ( 6  ,  °8 ) 

WRITE  (6  *  l.oo  ) 

!,'R  I  TF  (  6 , 1°1  ) 

OQl 

hRItf(6,10?) 

DO  ooi  J=1 ,NP 

WRITE  (6,1.  OS)  J  ,  T  (  J  )  ,  V  (  J  )  ,PROFL(  J) 

KK  =Q 

IK  =  ?0 

F ( KK  )  =11 8  ,F6 

F (  IK  ) =407,87 

1 1  ( 1  ) =00 . no 

FL00=4. 1 0 

9? 

I  v;  =  o 

IF ( IW.FO.O)  GO  TO  77 

GO  TO  6 E 

6F 

66 

DO  66  J  =  1 ,NP 

CALL  WILFON(T) 

GO  TO  78 
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f  oo 

Os'i'i 

rvs>l 

^  ^  r  f  -  o  'y? 

X  P  .  c  o  \  =  (  N  I  )  t 

r)n,cf  =  i  f  >n 

p  f  #  =  n  o  J  ~ 

n*wT 

on  (o.oq.wn^T 

SO 

P*  OT  oo 

q,‘ .  r  =  '  A  a  on 

(  r  )  :  <  T  !  M'O 

PT  OT  00 

77  no  7  Jal ,NP 
no  7  1=1 *NCP 
CALL  KVAL(T) 

7  COMTINUF 

78  no  1 7  T  =  1  ,NCP 

_ DO  10  J  =  1  »NP _ _ _ 

10  R ( J ) =- { ( V ( J ) +0 ( J ) )  *VK  (  I  *  J  )  +  (  PROFL  (  J  )  +1 1  (  J  )  )  ) 
f ( 1  )=o.nn 

_ DO  .11  _ _ 

JJ=J-1 

U  F( J)=PROFL( JJ) 

_ DO  1?  J=1,NPP _ 

1?  C(J)=V(J+1)*VK(  I  » J  +  l ) 

C ( NP ) =n . on 

_ DO _ 1 =JL„*7iE _ 

17  n ( j ) =o . on 

C 

C  I  RRPRFSFNTS  THF  COMPONFNT.  D  ic  THF  COMPOS  T  T  T  ON  0^  THF 

c  particular  fffd  in  molfs  of  componfnt  i  pfr  unit  timf. 
c  FOR  fxamplf*  the  solvfnt  ffed  ha s  nonf  of  COMPONFNT  1 
C  I ACETONE)  AND  THF  FRFSH  FEED  HAS  0.071 8  MOLES  OF  ACETONE 
C  PrR  rOLF  of  FFpn 
C 

T  F ( I »  QQ  .  1  )  op  TO  ft  1 


I  p  (  T  •  FQ . ?  )  00  TO  A  ? 
IF (  I  .RQ.7 )  00  TO  A  7 

on  jo  m 


61 

D  (  KK ) =-F ( KK ) *0.0 

D( IK) =-F ( TK)*0.07l R 

GO  TO  oi 

67 

D  (  KK  )  =-0 .0 

D (  T  K ) =  — F (  TK)*0.0o7? 

00  TO  91 

60 

n ( KK  )  =-F ( KK ) *  1  .00 
n(  IK ) =-F ( TK)*0.8710 

GO  TO  91 

91 

ROT  =  F ( KK ) +F ( TK)-U(  1  ) 

Pfl)=C(l)/R(l) 

0(1  )=D(1  )/R(l  ) 

14 

DO  14  J=?*NPP 

P{J)=C(J)/(R(J)-(F(J)*P(J-1 ) ) ) 

P(NP)=n.nn 

1  R 

DO  16  J=?*NP 

0(J)=(D(J)-(F(J)*Q(J-1 ) ) ) / ( R ( J ) — ( F ( J ) *P ( J— 1 ) ) ) 

X(  T ,NP) =0 ( N  P ) 

1  A 

DO  1  A  J=1  * N P P 

L*NP— J 

X(I,L)rQ(L)-(P(L)»X(M+D)   .  . .  .   ... 

17 

CONT  INI  IF 

TF(IW.GF.4)  00  TO  707 

GO  to  R07 

c 

c 

c 

THF 

TO 

ACFTONF  CONTFfMT  IN  THF  TWO  ROTTOM-MOST  STREAM  If  FFT 

Z.RRO  to  TNDIICF  CONVFROFNCF. 

c 

300 

V  (  1  ,NPP)=o.oooo 

X  (  1  .NP ) =  0 . 0 0 0 0 

c 

c 

c 

THE 

COMPOSITION  PROFILF  IS  BEING  NORMALIZED. 
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70?  no  j=i ,np 
Qi  imx  ( J  )  =  o.or> 
no  IQ  i  =  i, NCP 
SUMX ( J ) =SUMX ( J )+X ( T  »J) 

1°  CONTINUF 

_ no  Q7  J= 1  <NP _ 

no  Q 7  1=1  ,NCP 
X ( T  ,  J) =X ( I  9 J ) /FUMX ( J ) 

CONTINUE 

IF  c  1W.EQ.01  00  T n 

GO  TO  71 

6Q _ np  7 n  J  =  i ,NP _ 

7n  GALL  WILFON(T) 

T W= I W+l 
JJJ=  T # 

GO  TO  78 
71  IW=IW+1 

J  J  J  =  I  w _ 

IF ( JJJ.F0.7Q )  GO  TO  8Q 
GO  TO  21 1 
211  gall  TNEWf At) 

Cl  lMniT  =  n  #  f) 

no  47  J  = 1 ,  N  P 

47  FUMnT  =  SUMnT  +  ( AT ( J)-T ( J)  )**? _ 

IF  (  FllMOT  .LF.  1 .0  )  GO  TO  78 
GO  TO  48 

40  DO  lx o  j  =  i  ,mp 
AO  T  (  J  )  =  A  T  (  J ) 

CALL  PROFYL(T) 

GO  TO  66 _ 

78  RFCT=X ( ? ,NP )*R0T/47 .007 
V'RITE(6*  116  ) 

WR ITE ( 6  » 106 ) 

WRITE (6 *101 ) 

VR  ITF (6  9 1 Q?  ) 

no  Q 0?  J=1,NP _ 

007  W  R I T  F ( 6  *  1 0  7 )  J  ,  T (J)  ,V(J)  ,  PROFL ( J ) 
'•'R  ITF  (8 , 1  04)  lW*EUMnT 
WRITF(6 , Til) 

WRITE (6  » 1 1 
no  pn 7  1  =  1, NCP 

QOT  VR1TF(F,117)  X (  T  ,  1  )  ,x (  I  ,NP ) 

'»'R  ITF(6»114)  RFCT 
6Q  QT  OP 
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TRFTC  RAV 

oi  IRROI  IT  I NF  KVAL(T) 

COMMON  J  » IW  *NP  *NfP  ,K^  ,  TK  .FLOD,  I 

COMMON  V  (  2  6  )  .PROFL  (7*)*VK(16,26)*X(1S,?6) 

COMMON  SHMX ( 7  5 )  >Y(  1 6 , 2  6  )  *RRMIX ( 26 )  .7C f  7  6  )  *  A  AOR ( 7 6 ) 

c 

COMMON  A  (  1 7 *6  )  *DFLTA (  16 )  *V1  (  1  *  )  *DDF  LMX  (26)  ,F(?5)  ,IJ(  2  6  ) 
DTMFNFION  FN1 ( 10.76 ) *T(?6 ) 

c 

c 

c 

TH  I  c  FURROllTINF  CALojLATFF  COMPOb  t  T  TON-INDEPFNDFNT  VAPO¬ 
RIZATION  CONST  ANT  6  FROM  THE  CHAO-SEADFR  MODFL  RY  ASSUMING 
K(I,J)=F/P  RATIO  for  THF  LIQUID  PHASF. 

c 

p=l 4.7 

A  0  =  6 .76748 

A  1 =-7.01761 

A  2  =— 4 .PSSpo 

A^=2 . 072QO 

A4=o . OO 

A  6  =  n  ,  08  4  7 7 

A6=o. 26667 

A  7  =  — 0 . 7 l  1 78 

AR=-0. 02666 

A9=0 . 028R7 

A  1  0  =  — 4 • 27897 

AT  1 =8 • GK  708 

A  1 2  =  - 1  .77060 

A1 7=-7. 15774 

A1 4=-0.026 

IF(I.EQ.l)  GO  TO  U 

IF (I .EQ.2)  GO  TO  12 

IF ( I .FQ.7  )  GO  TO  17 

11  TC=°1 4.0 

'•'  =  0.616° 

Pr=6°o. o 

GO  TO  16 

12  TC =977,0 
w=n,60Q 

PC  =  1 1 6  6.4 

GO  TO  16 

17  TC=1 164,7 
'•'=0.6676 

PC  =  72or>  .  ? 

GO  TO  1 6 

16  TR=(T( J)+46Q.0)/TC 

PR-P/PC 

FPl =A0+A1/TR+A7*TR+A7*( TR**? )+A4*( TR**7 )+ ( A6+A6*TR+A7* ( TK**2 ) ) *PR- 
1 { A8+AQ*TR)*(PR**7)-ALOG(PR) /7.7026R 

FP2  =  A1 0+ ( A1 1*TR )  +  ( A1 2 /TR )  +  ( A1 7* ( TR**7 )  )  +  ( Al 4* ( PK-0.6  >  ) 

FN 1  (  I  » J  )  =  { F P 1 +  W * F P  2 ) *2.^0268 
'/K  (  I  ,J)=FXP(FN1  (  I  ,  J)  ) 
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MRFTC  LTLI 

SNRROUTINF  WTLSON(T) 

COMMON  J  ,  IW*NP  ,NCP,<K « I<  *  F  LOO  ,  I 

COMMON  V(  ?5  )  *PROFL  (  2*  )  *VK  (1  *>  *?5  )  ,X  (1  s  *  25  ) 

_ COMMON  SUMX  (  ?  5  )  >  V  (  1  5  ,  ?  S  )  »RRMIX  (  ?5  )  »Z  C  (  ?S  )  »  AAO«  (  ?S  ) 

COMMON  A  ( 1  7  *5  )  *  DFLTA  (  1  s  )  *  VI ( 1 s )  >DDFLMX(?S)  ,F  (  25  )  *U ( ?s ) 
DIMENSION  ALAMDA ( ^ , S )  .  7 ( ?5 )  »C ( ?S )  *RHO(15).R(15)  .WILSPR 
1  (1C.1C)»G(1S»?5)  *GAMM  AJ  is  *  V.P  (  1.5.*  2..5J  *.5J  .'MV  ( ?6  )  « 51?  FLY. 
2(15)  *SUM1  ( 15 )  »  SUM?  ( 1  5  )  »Sl )M ^  (  i  c  )  »  T  (  ^  } 

C 

C  THIS  SUBROUTINE  CA  Ln  >[_  a  tf  s  COMPO  s  I  T  T  ON-DFP  FNPFNt  vaPQRT- 

C  ZATION  CONSTANTS  USING  THE  SIMPLIFIED  PRAUSNIT7  FQm  I L I - 
C  BRIUM  MODEL  WITH  LIQUID  ACTIVITY  COEFFICIENTS  CALCULATED 
C  GY  T H s  WILSON  MODEL 
C 
C 

C  COMPONFNT  1  IS  ACETONE  . 7  is  METHANOL  » ^  Ts  WATFR. 

C 

DO  7  JI=1,NCP 
no  7  13*1 »NCP 
7  ALAMDAt JI * J J)*0.0 
N  =  1 

_ M=? _ _ _ __ _ _ _ 

L  =  s 

C 

c  ALAMDA  Ic  THE  BINARY  INTERACTION  ENFRG1FS  IN  CAL/GM-MO 
C  OBTAINED  FROM  ONE  OF  THr  OPTIMIZATION  TECHNIQUES 
C 

_ ALAMDA(N»L)=-164.879 _ 

ALAMDA ( L  » N ) =1679.65? 

ALAMDA ( N  *M ) =-315.893 
ALAMDA (M »N ) =719> 198 
ALAMDA ( M*L ) =  24  1 .0^7 
ALAMDA ( L  *M ) =  S 84 • 8 7 0 
_ c_K.n  =  IIJ  Jj-s?  .oj_/ L«P. 

Zf J)=C( J  J+273. 16 
C 

C_  RHO  is  THE  DENSITY  CALCULATED  BY  EQUATIONS  OFTEN  REPORTED 
C  IN  THF  CRITICAL  TABLES.  B  IS  TH^  MOLAR  VOLUME  OF  THE 
C  COMPONENTS ( MOLE  WT/DENSTTY) 

C _ 

RHO(N)=n#pi?&p-i  ,inn*c(J)*(10.**(-'3)  ]-o,«e;q^(r(J)**7)* 

1 (10. **(-6) ) 

_  RH0(M)sQ»792 _ 

0 ( L )  =  1  .000 
n ( N ) =58 . 046 /RHO ( N ) 

_ b(MI)=32.031/RHO(M) 

R ( L ) = 1 8  •  /RHO ( L ) 

C 

C _ WILSPR  IS  THF  BINARY  W I L SON  CONSTANTS. 

C 

no  i ?  11=1 ,ncp 
DO  1?  J I =1 » NCP 

WILSPR (  I  I  > JI  ) =P ( JI  ) / ( R ( I  I ) *FXP ( ALAMDA (  I  I  * J I  ) / ( i.98  70*z" 

1  ( J  )  )  )  ) 

12  CONTINUF 

Z 

C  TP  TS  TOTAL  PRES  SURF  IN  MM-HG. 
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C 

c  thf  following  is  thf  calculation  of  the  liquid  pha^e 

C  ACTIVITY  COEFFICIENTS  ACCORDING  TO  THF  WILSON  EQUATION. 

c 


11  =  1 


21 

17 

SUM! (111=0.0 

DO  17  J I =1 ,NCP 

SlJMl  (  I  T  )  =  SUM  1  (  I  I  )+X  {  JI  ,J)*WILSPR(  I  I  ,  JT  ) 

DO  19  I  I  I  =  1  , NC P 

SUM? ( I I I ) =0.0 
no  18  J 1=1  *  NCP 

1  P 

19 

SUM? ( I T  I  ) =SUM? (  I I I  )+X ( JT  * J  > *w  TLSPR (  I  I I  ,  JI  ) 

SUM? (  I  I  I  )  =  ( X (  I  I  I  * J )*WILSPR (  I  I  I  ,  I  I  )  ) /SUM? (Ill) 

S|  im  =  o  .  o 

20 

DO  ?o  JJJ=1,NCP 

SI.JM  =  SUM  +  SUM?  (  JJJ  ) 

G  (  I T  »J)=-ALOG(SUMl  (IT)  ) +1. O-SUM 

IF (  I  I .FQ.NCP)  GO  TO  ?? 

11=11+1 

GO  TO  21 

2? 

27 

DO  ??  11=1  » NCP 

GAMMA ( T  T  * J ) =FXP ( G ( 1 1 , J )  ) 

CONTINUE 

C 

C 

c 

VP  IS  THE  VAPOR  PRESSURE  CALCULATED  BY  THE  ANTOINE  EQUA¬ 
TION  AND  REPORTED  IN  STANDARD  HANDBOOKS. 

c 

VP ( N » J) =FXP( ? . ?0?sp* ( 7. 0?447- ( l l .0/ ( ?  ? A . 0+C ( J)  )  )  )  ) 

VP  (M,  J  )  =FXP  (  ?.  ?0?SR*  (  7.P7  8S?-  (  1  A  7?  _  1  1  /  f  ?  ?  o  .  o  4  r  <  .i)  )  )  )  ) 

?A 

VP ( L  * J ) =FXP ( ?. ?0?SR* ( 7.98681  - ( 1  AS 8.? 1 / ( ??P .0+C ( J )  ) 

no  ?f>  i  =  i  » ncp 

Y(  I ,J)=X (  I  * J)*GAMMA (  I  « J)*VP<  I  ,J) /TP 

)  )  ) 

2  5 

DO  25  1=1 »NCP 

VK (  I  ,J)  =Y( I  ,  J) /X (  I , J) 

RETURN 
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SIBFTC  FOWARO 

SURROUTINF  TNFW(AT) 

COMMON  J , IW ,NP ,NCP  *KK  ,  IK  ,FLOD  ,  I 

COMMON  V  (  )  ,PROFL  (?‘i)»VK(l^?5),)((lSl?'i) 

_ COMMON  SIJMX  (  ?F  )  ,Y  (  IF  ,  ?F  )  ,RRMIX  (  2F  )  ,7C  (  ?F  )  ,AAOo(  ?F  ) 

COMMON  A(17,M  *  D  F  L  T  A  (  1  F  )  *  V 1  (IF)  ,ODELMX  (  ?F  )  ,F ( ?F  )  ,i)(  ) 

OIMFNSION  TT ( ?F  )  »AT(?f)*D(?F)  ,S ( ?n , pn )  ,SUMY ( ?F ) * GAMMA ( 

11  F ,?F ) »T(25  ) 

J=1 

87  JJ=1 

16  0  =  J  J _ _ _  _ 

TT(JJ)r(o,c;#o  +  n,^)#ion,o 
T ( J  )  =TT ( JJ) 

_ CALL  WILSON (T) 

SUMY! J) *0*00 
DO  1?  I =1  » NCP 

12  SUMY ( J ) =SUM Y ( J ) +  V< ( I , J)*X(  I  ,J)  _  _ 

S( J,JJ)=SUMY( J)-l .00 
JJ=JJ+1 

IF (JJ.FQ.4)  GO  TO  IQ 
GO  TO  1 4 

IQ  K  =  4 

QQ  HtK-1  ) = ( TT ( K—l  )  — TT(K— 2)  )/(TT (K~2)~ TTfK-8)  ) _ _ 

p=S(J,K-p)*(D(K-i  )  **2  )  -S ( J ,K-? )* ( i  . 0+n ( K-l  )  )**2+S(J,K~l  )*  (  ]  . 0+ (2  . 
ID ( K-l ) ) ) 

C=S(  J,'<-p  )*0(  K-]  )-S(  J,K-2)*(  1  .0  +  0  (  K-l  )  )+S(J,K-l  ) 

00*-2.0*S(  J,K-1  )*(l  ,o  +  D(K-U  ) 

AA=( APS ( R ) **? ) -4.0*S ( J,K-1 ) *0( K-l ) * ( l .o+D( K-l )  )*C 

_ &p.P  =  p  +  (ARS(AA)**0,F) _ _ _ _ _ _ _ _ _ 

APM=R_ ( APS ( A A ) **0. F ) 

PCI =00/ A BP 

BC2gQ0/ABM  _ _ _ 

TA  =  TT(K-1)  +  (TT(K-1  ) - T  T ( K - ? )  ) *  R  C 1 
T°  =  TT (K— 1 )  +  ( TT (K-l  ) —  T  T ( K  — 2 )  )*BC2 

_ IF ( TA.GT.0.0 )  GO  TO  6Q  _  _ _ 

GO  TO  70 
6Q  T ( J  )  =T A 

_ CALL  WI LSON ( T ) 

St ,MY  f  J)  *0*00 
DO  81  1=1, NCP 
81  SUMY ( J) =SUMY( J)+Y( T ,J) 

SA=SUMY( J)-l .00 
IF(TR.LT.O.O)  GO  TO  Q2 
GO  TO  70 
70  T(J)*TB 

CALL  WILSON (T) 

_ SUMY  (  J  )  =0.0 _ _ _ 

DO  88  1=1, NCP 

88  SUMY ( J ) =SUMY( J )+Y ( I , J ) 

SB«SUMY ( J ) -1 .00 

!F(ABS(SR).LT.ABS(SAM  GO  TO  q i 
GO  TO  Q  2 

Q 1  IF ( ARS( SR) .LP. 0.0001 )  GO  TO  Q7 
TT (K-l ) =  T  P 
S ( J  ,K-1  ) =QR 
GO  TO  QQ 

9 2  IF (ABSI SA> .LE. 0.0001 >  GO  TO  QR 
TT (K-l ) =  T  A 

S  (  J  ,  K  - 1  )  =  S  A  
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GO  TO  oo 
Q7  AT ( J ) =  T  R 
GO  TO  1 6 
08  AT(J)=TA 
GO  TO  IS 

IF(J.FQ. NP)  GO  TO  56 

GO  TO  8R 
50  J=J+1 

_ fiO  J  O  _.5J _ _ _ 

56  RETURN 
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$TRFTC  KATHY 

SUBROUTINE  PROFYL(T) 

COMMON  J  *  I W  *NP  *  NfP  *  KK  ,  IK*EL0D,' T 

COMMON  V ( 2* ) »PROFL (  ?s  )  ,VK( is  ,  )  ,X  (1  5  ,  2S  ) 

_ COMMON  SUMX(?5)  ,  Y  (  15,25)  ,PRMIX ( 25 )  .ZC ( 25  )  ,  AAQR ( ?5 ) 

COMMON  A  (  1  7  *  5  )  *  DEL  T  A  (  IF)  ,V1  (if)  ,nPFLMX ( 25 )  , F ( 25  )  ,i|  (  ) 

DIMENSION  SUMY( 2F )  ,SH( 15  )  ,CH( IS )  ,EV( ?F )  ,EL ( ?f )  ,q 

_ 1  (  ?F  )  «HEATDLC15»15)  *FXLl  (15*25)  *EXL7 (15.75)  .FXCSHT  ( 1 «  . 

25 ) *XF ( 1 5  *25 ) ♦ MF ( 75 ) >6 ( 25  >  * T ( 25  > 


C 

C  W,S  ARF  MOLFCMLAR  WEIGHT S 


c 

W 1 =58.046 
’•'2  =  32.031 

W 3=1 8.0 

DO  73  1  =  1, NCP 

DO  78  11=1 , NCP 

73  HF ATDL (  I  *  I  I  )  =  0 • On 
C 

C  HFATDL  ARF  HEATS  OF  SOLUTION  AT  INFINITE  DILUTION  IN _ 

C  PTU/LP-MOLF  OBTAINED  FROM  THE  CRITICAL  TABLES. 

C 

HFATDL (1.3) =4520.0 

HFATDL (3,1) =4520.0 
HFATDL ( 2  *3 ) =?6os.o 
HFATDL  r*  .2 ) =**n3 . o 
C 

C  THE  FOLLOWING  IS  THF  CALCULATION  OF  PARTIAL  MOLAR  HEATS 
C  OF  solution  AT  THE  plate  compositions  USING  EQUATION _ 

C  1 1-22  IN  B.D. SMITH. 

C 

DO  op  J=1 ,NP 
DO  op  1=1 ,NCP 
FXL]  (  I  ,  J  )  =0.0 

_ DO  0  8  IJ  =  1  ,3 _ _ _ 

FXLl ( I  * J )=X ( I J* J )* (HFATDL ( I J, I ) +HF A TDL ( I ♦ I J )  ) +FXL 1 ( I , J 
1  ) 

OR  continue 

DO  00  j= 1 ,NP 
DO  00  1  =  1  ,NCP 

_ EXL2 ( I , J ) =r . 0 _ 

DO  oo  iJ=i , 8 

IF ( I J.EQ. I )  GO  TO  09 

IF ( I J.EQ.3 )  GO  TO  90 

E XL? (  I » J ) =X ( 1  , J ) *X (  I J  * J ) * ( HFATDL (  I J , l  ) +HEATDL (1,IJ))+FXL?(I,J) 
99  CONTINUE 

_ DO  19  J  =  1  ,NP _ _ 

DO  19  1  =  1, NCP 

FXCSHT (  I  ,J) =0.5* ( 1 .o-X ( I , J)  )*FXLl ( I ,J ) -O.S*EXL7  (  I  ,  J  ) 

_ 10  _CONtjnuf _ _ 

C 

C  SH  AND  CH  ARF  LIQUID  AND  VAPOR  ENTHALPIES  RESPECTIVELY. 

C  THE  FIRST  PART  OF  FACH  FQUATI0N  IS  A  LINFAR  REPRESENTATION 

('  OF  THE  IDEAL  ENTHALPIES. 

c 

DO  74  J=1 » NP 
1  =  1 

SH ( I , J ) =0.5279*T ( J ) *W1 +FXCSHT ( I , J ) 

CH( I *J) =SH( T  *J)+2 24.^ *W1 
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74 

SH( I+l • J )*0*5996*T ( J) *W2+EXCSHT ( t + 1 
CH  (  T  +  1  »  J  )  =  .SM(  I  +  • 

SH  (  T+?,J)=0.ooQ6*T(J)  *W**  +  FXCSHT  (  T  +  7 
FH(  I  +  ?  9  J )  =  S  H  (  I  +  ?»J)+O47.06*WR  ’ 

DO  67  J=1 »NP 

8 !  )M  Y  (  J  )  =  n  .  a 

9  J  ) 

♦  J) 

66 

DO  66  ! =1 *NCP 

S IJ M Y  (  J  ) =SUMY( J)+Y<  I  ,J) 

DO  67  1=1 ,NCP 

67 

Y( I , J)=Y ( I , J) /SUMY ( J ) 

DO  76  J=1,NP 

FV ( J ) =0  •  0 

FL  ( J  )  =0.«  00 

DO  76  1=1 ,NCP 

FV( J)=FV( J)+CH(  I . J)*Y(  I  ,J) 

76 

FL ( J ) =FL ( J ) +SH ( I . J ) *X ( I « J ) 

DO  61  J=1,NP 

0( J ) =8. on 

51 

G  (  J  )  =  0 . 0 

HF ( J ) =  9  •  0 

DO  68  J= 1 »NP 

c 

68 

DO  68  1  =  1 » NT P 

XF (  I  » J ) =0 . on 

c 

c 

XF 

IS  FFPP  COMPOSITION* 

XE,(  R  »,KK  )  =1*0  . .  ....  . . 

XF (1  *  IK  )  =0.0718 

XF ( ?  » IK ) =0.8077 

XFJR  9  TK1  =o.8R1  0 

PROFL  (NP  )=F  (  KK  )+F.(  IK)  -U(  1  ) 

PROFL  (  1  )  =FLOD*U(  1  ) 

_ V  ( ? ) = PROFL ( 1  ) + U ( 1 ) 

C 

c 

c 

HF 

IS  FFED  ENTHALPY. 

DO  69  J=1 9NP 

1  =  1 

H  F  (  J  )  =  n  .  6  7  7  9  *  T  (  J  )  *  W 1  *  X  F  (  I  ,  J  ) 

69 

HF ( J) =HF ( J) +( 0.8Q96*T ( J )*W2*XF ( I  +  i • 
HF(J)=HF(J)  +  (0.9OOS*T(J)  *'/'P*XP  (  1+7  , 
NPP=NP-1 

J)  ) 

J)  ) 

DO  79  J  =  2  9  NPP 

XRX= ( FL ( J-l  ) -RL ( J)  ) *PROFL ( J  — 1  ) 

R  X  X.=  .(.  EL  (  J  )  -FV  (  J  )  )  *  (  V  (  J  )  +GJ  J  )  ) . 

R  D  X  =  (HF  (  J)— F(_(  J)  )  *  F  (  J  ) 

Prp  = (FL ( J ) — FV ( J  +  l )  ) 

V ( J+l ] = (  XRX+PXX+QRX-OtJ)  )/P.QR 

7Q 

PROFL  (  J  )  =F  (  J  )  +PROFL  (  J-l  )+V  (  J+l  )-v  (  J  )-0  (  J  )-l)  (  J  ) 

RFTl  IRN 

END 
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Cards  to  be  changed  when  problem  VI  is  used  with  a  different 


system: 


Mainline  program  (Thesis) 


1. 

Number  of  components 

NCP  =  3 

2. 

Number  of  plates 

NP  =  25 

3. 

Feed  introduction  points 

solvent  KK  =  9 

fresh  IK  =  20 

4. 

Feed  flow  rates,  moles/U.T 

.  F  (KK)  =  115.56 

F  ( IK)  =  492.87 

5. 

Distillate  product  rate,  moles/U.T.  U(l)  =  33.09 

6. 

L^/D  ratio 

ELOD  =  4.13 

7. 

Compositions  of  the  feeds 

D  (KK)  =  -F  (KK) 

*0.0 

D ( IK)  =  -F(IK) 

*0.0718 

0.0  is  concentration  of  component  1  in  solvent 
0.0718  is  concentration  of  component  1  in  fresh  feed 
This  is  repeated  for  each  component. 

8.  Recovery  calculation 

RECT  =  X(2,  NP)  *BOT/47 .907 

where  47.907  is  the  moles/U.T.  of  component  2  in 


the  fresh  feed 
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$  IBFTC  Ray 


Ideal  enthalpies  are  calculated  from  Chao-Seader  by  assuming 
=  f/p  ratio  for  liquid  phase. 

1.  Column  pressure,  psi  P  =  14.7 

2.  Critical  temperature  TC 

Acentric  factor  w 

Critical  pressure  PC 

For  every  component 

$  IBFTC  Lili 


(Wilson) 


subroutine 


1.  The  binary  interaction  energies 
AL AMDA ( N , L )  etc. 

where  N  =  component  1 

M  =  component  2 

L  =  component  3 

2.  Parameters  for  the  temperature  dependent  densities 

RHO(N)  etc 

3.  The  molecular  weights 

B (N)  =  5  8 . 04  6/RHO (N) 

where  B(N)  is  molar  volume  for  component  N. 

4.  System  pressure,  mm  Hg.  TP  =  760.0 

5.  Vapor  pressure  equations 

VP  (N , J)  etc. 

The  parameters  are  those  for  the  Antoine  equation  (Obtained 
from  Lange's  Handbook). 
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New  Rates  Subroutine 

1.  Molecular  weights  W1 ,  W2 ,  W3 

2.  Heats  of  solution  at  infinite  dilution 

HEATDL (1,3)  etc. 

These  are  obtained  from  the  critical  tables. 

3.  Ideal  enthalpy  equations 
liquid  enthalpy 

SH(I,J)  =  0.5279  *T(J)  *W1  +  EXCSHT(I,J) 

vapor  enthalpy 

CH  ( I ,  J)  =  SH(I,J)  +  224.3  *W1 

where  224.3  is  the  heat  of  vaporization  in  calories , 

gm 

0.5279  is  heat  capacity  in  cal _ 

gm  -  degree  F 

This  is  repeated  for  each  component. 

4.  Feed  compositions 

XF ( 3 , KK)  etc. 

This  is  composition  of  component  3  in  the  solvent  stage  KK. 

5.  Feed  enthalpies 

The  heat  capacities  of  the  components  should  be  changed 
in  the  calculation  of  HF(J). 

Data  Cards  to  be  Changed 

1.  Initial  profile 

X01050. 0152.0000. 0 
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01 


Stage  number 


Temperature 

=  F5.1  = 

050.0 

°c 

Liquid  rate 

=  F5.1  = 

152.0 

moles/U . T . 

Vapor  rate 

=  F5.1  = 

000.0 

moles/U . T . 

This  is  for  stage  1  (Condenser)  and  is  repeated  for  every 
stage . 
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